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Abstract 


ie  compatibility  of  a representative  500  pound  weight  boost- 
glide  air-to-ground  missile  with  the  trajectory  constraints  imposed 
by  a Synthetic  Aperture  Radar  - Retransmission  guidance  system  was 
investigated  using  a digital  flight  simulation.  A demonstration 
flight  profile  was  assumed,  with  a minimum  of  20  seconds  of  tracking 
on  the  aircraft- to- target  line  of  sight  required.  A guidance 
algorithm  was  developed  which  produced  satisfactory  trajectories. 

A first  order  gradient  technique  was  employed  in  an  unsuccessful 
attempt  to  optimize  the  trajectories  for  maximum  range.  A useable 
launch  envelope  for  this  missile  was  determined.  The  azimuthal 
extent  of  the  envelope  was  limited  by  radar  system  constraints  to 
15*-90  from  the  aircraft  velocity  vector.  A mavtmimi  slant  range  of 
21  nautical  miles  was  obtained  from  a launch  altitude  of  35,000  ft. 
Range  deteriorated  rapidly  with  decreasing  launch  altitude,  with 
5,000  ft  being  the  lowest  altitude  at  which  a useable  launch  envelope 
was  obtained.  Maximum  slant  range  at  5,000  ft  launch  altitude  was  5.5 
nautical  miles. 
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PERFORMANCE  OF  AN  AIR-TO-GROUND  MISSILE 
EMPLOYING  SAR-RETRAN  GUIDANCE 

I.  Introduction 

Background 

Synthetic  Aperture  Radar  (SAR)  is  a coherent,  pulse-doppler  radar 
which  la  capable  o£  finer  azimuth  resolution  than  conventional  non- 
coherent radars.  It  is  also  capable  of  distinguishing  small  surface 
targets  from  background  clutter  and  can  track  such  targets.  With 
these  capabilities,  SAR  is  a prime  candidate  for  an  all-veather  guidance 
system  for  tactical  air-to-ground  missiles.  However,  the  operation  of 
SAR  imposes  some  unique  trajectory  requirements  on  the  missile. 

SAR  attains  fine  azimuth  resolution  by  utilizing  a relative  velo- 
city discrimination  approach  (Ref  1:  4.d.(2)-3).  An  array  of  points 
located  at  the  same  range  R from  an  aircraft-mounted  SAR  is  depicted 
in  Figure  1.  Each  of  these  points  has  a velocity  relative  to  the 
radar  antenna  of 

V - VA  cos  0 (1) 

where  is  the  velocity  of  the  aircraft  and  0 is  the  angular  displace- 
ment of  the  point  from  the  aircraft  velocity  vector  (squint  angle). 

This  relative  velocity  will  produce  a doppler  frequency  shift  in  the 
reflected  radar  signal  of 

„ 2Va  cos  0 

£< X (2) 


1 


2 


One  technique  for  guiding  an  air-to-ground  missile  using  SAR  could 
be  to  locate  and  track  the  target  with  an  aircraft-mounted  SAR  and 
command  guide  the  missile  to  the  target.  Under  this  scheme,  the  terminal 
accuracy  of  the  system  could  be  no  better  .than  the  minimum  resolution 
of  the  SAR,  which  decreases  as  the  aircraft's  range  from  the  target 
increases.  It  would  be  more  desirable  to  have  a convergent  guidance 
system  whereby  resolution,  and  hence  accuracy,  can  be  refined  as  the 
missile  approaches  the  target.  One  way  of  obtaining  such  convergent 
guidance  is  through  a technique  called  retransmission  (SAR-Retran) . 

With  SAR-Retran,  the  missile  is  command  guided  into  the  radar  beam 
while  an  antenna  aboard  the  missile  receives  the  reflected  radar  signal 
and  retransmits  this  signal  to  the  aircraft.  When  the  missile  is 
squarely  on  the  aircraft-to-target  line  of  sight,  the  doppler  phase 


shift  of  the  retransmitted  signal  will  be  identical  to  that  of  the 
direct  reflected  signal  received  by  the  aircraft's  antenna.  Now  the 
SAR-Retran  system  can  "see"  the  target  through  the  missile's  antenna 
as  it  tracks  the  line  of  sight  to  the  target.  As  the  missile  closes 
on  the  target,  the  system's  resolution  increases,  concurrently 
improving  terminal  accuracy.  It  is  important  for  the  missile  to 
precisely  track  the  line  of  sight  in  the  terminal  phase,  for  if  it 
is  displaced  from  the  line  of  sight,  an  additional  phase  shift  will 
result  (more  fully  discussed  in  Section  XI)  causing  inaccuracy  in 
target  resolution. 

The  requirement  for  aircraft  velocity  vector  offset  causes  the 
aircraft-to-target  line  of  sight  to  rotate  about  the  target  as  the 


i 
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aircraft  files  by.  Thus,  the  missile  must  Intercept  this  rotating  line 
of  sight,  then  must  track  the  line  of  sight  for  a minimum  period  of 
time  In  the  terminal  phase. 

Problem 

The  Air  Force  Avionics  Laboratory  is  considering  available  air-to- 
ground  missiles  for  flight  test  and  demonstration  of  the  SAR-Retran 
system.  However,  the  trajectory  constraints  imposed  by  the  SAR-Retran 
technique  far  exceed  the  original  design  criteria  of  the  missiles  in 
the  current  inventory.  Therefore,  the  aerodynamic  compatibility  of 
these  missiles  with  SAR-Retran  must  be  demonstrated  by  computer 
simulation  before  flight  testing  is  undertaken  (gef  2:1). 

Purpose  of  Study 

This  study  investigates,  through  computer  simulation,  the  ability 
of  a representative  boost-glide  air-to-ground  missile  to  fly  the 
trajectories  required  by  SAR-Retran  guidance.  The  constraints  which 
must  be  satisfied  are: 

1.  A minimum  of  20  seconds  tracking  of  the  line  of  sight 
prior  to  target  impact.  The  missile  is  considered  "on" 
the  line  of  sight  when  the  retransmitted  frequency 
error  is  lees  than  or  equal  to  8 hz.  .(Constraints 
imposed  by  radar  and  guidance  system.) 

2,  Terminal  velocity  at  target  impact  must  be  greater  than 
700  ft/sec  (minimum  missile  maneuvering  speed). 

An  attempt  was  made  to  maximize  launch,  range  by  applying  the  results 
of  a first  order  gradient  optimization  scheme  to  the  guidance  algorithm. 


Range-altitude  envelopes  from  which  the  aforementioned  constraints  can 
be  satisfied  were  determined. 


Scope 

1.  A hypothetical  missile  characteristic  of  a boost-glide 
air-to-ground  missile  will  be  considered.  This  missile  will  be 
referred  to  as  the  Retran  Missile. 

2.  This  study  is  limited  to  the  aerodynamic  performance  of 

the  representative  missile.  It  is  assumed  that  all  other  elements  of 
the  guidance  system  operate  perfectly. 

3.  The  study  concentrates  on  determining  the  mid-course 
trajectories  and  required  guidance  algorithms.  Tracking  of  the  line 
of  sight  in  the  terminal  phase  was  accomplished;  however,  terminal 
guidance  and  accuracy  were  not  among  the  goals  of  this  thesis. 

4.  Optimization  of  the  midcourse  trajectory  to  achieve  maximum 
launch  range  was  attempted.  The  results  of  this  optimization  were 
considered  in  formulating  the  guidance  algorithms;  however,  the 
resulting  trajectories  are  not  necessarily  optimal  for  any  launch 
conditions. 

5.  Launch  envelopes  were  determined  for  demonstration  flight 
conditions.  The  target  was  stationary  and  the  aircraft  did  not  make 
tactical  maneuvers. 

6.  The  stability  of  the  guidance  loop  with  the  missile  autopilot 
was  investigated.  Modifications  in  commands  and  biases  provided  to 
the  autopilot  were  made  but  the  circuitry  and  mechanization  remained 
unchanged  from  the  standard  autopilot. 
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II.  Analytic  Formulation  of  the  Missile  Flight  Simulation 


The  Retran  Missile  simulation  program  numerically  Integrates  the 
missile  equations  of  motion  over  the  time  of  flight  to  determine  the 
time  history  of  missile  velocity,  attitude,  and  position.  The  program 
also  updates  the  launching  aircraft's  position,  computes  guidance 
commands,  simulates  the  missile  autopilot  to  determine  control 
deflections,  calculates  the  retransmitted  doppler  frequency  error 
and  determines  when  target  impact  is  accomplished.  The  mathematical 
relationships  used  in  the  program  are  presented  here,  while  the 
program  itself  is  discussed  in  Section  IV. 

Basic  Assumptions 

The  following  assumptions  are  incorporated  in  the  Retran  Simulation: 

1.  The  target  remains  stationary. 

2.  Perfect  target  tracking  is  achieved  by  the  radar  system. 

3.  The  aircraft's  radar  antenna  gimbal  angles  are  physically 
limited  to  ±95  degrees. 

4.  Missile  accelerometers  and  rate  gyros  are  error-free. 

5.  Through  data  link  of  missile  accelerometer  and  rate  gyro 
outputs  back  to  the  guidance  computer  aboard  the  aircraft, 
the  computer  can  determine  accurate  missile  position, 
velocity,  and  attitude. 

6.  Missile  roll  control  operates  perfectly. 

G 
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The  Missile  Model 

The  hypothetical  missile  used  in  this  study  will  be  referred  to  as 
the  Retran  missile.  This  missile  is  pictured  in  Fig.  2.  Its  size, 
aerodynamic  characteristics  and  performance  are  representative  of  a 
boost-glide  air-to-ground  missile  in  the  500  lb.  weight  class. 

The  physical  dimensions  of  the  missile  are  listed  in  Table  I,  the 
mass  properties  are  presented  in  Table  II  and  the  thrust  schedule  of  the 
two-stage  solid  propellant  motor  is  given  in  Table  III. 


Table  I 

Physical  Dimensions  of  the  Retran  Missile 
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Overall  Length  (In.) 

120.0 

Maximum  Diameter  (In.) 

12.3 

Wing  Planform  Area  (In.2) 

475.0 

(2  Wings) 

Control  Surface  Planform  Area  (In.2) 

105.0 

(2  Surfaces) 

Aerodynamic  Reference  Area  (In.2) 

.825 

Table  II 

Mass  Properties  of  the  Retran  Missile 


Time  (sec.) 

0.0 

0.6 

4.5 

Weight  (lb.) 

500.0 

460.0 

400.0 

Center  of  Gravity  (in.) 

55.0 

54.0 

52.0 

Ijot  (slug-ft2) 

2.5 

2.5 

2.5 

Iyy  (slug-ft2) 

67.0 

64.0 

58.0 

Izz  (slug-ft2) 

67.0 

64.0 

58.0 

Table  III 

Thrust  Schedule  of  the  Retran  Missile 

Burn  Time  (sec.)  Thrust  (lb.) 

1st  Stage  0.0-0. 6 9500 

2nd  Stage  0.6-4. 5 2000 

Note:  Thrust  Is  assumed  to  be  constant  throughout 

burn  time. 
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The  missile  is  aerodynamically  symmetrical  about  the  longitudinal 
axis.  Therefore,  the  normal  force  and  pitch  moment  coefficients  are 
used  for  both  pitch  and  yaw.  The  aerodynamic  coefficient  curves  are 
presented  in  Figs.  3 thru  11.  These  coefficients  are: 


% 

ACof 

Cc 

ACC 

C»2 

CNX 

**2 


- Zero  lift  drag  coefficient  at  sea  level 

- Skin  friction  drag  coefficient  (corrects  Cj)q  for  altitude) 

- Chord  force  coefficient  with  no  control 

- Incremental  chord  force  coefficient  due  to  control 
(Cg  + ACc  account  for  induced  drag) 

- Normal  force  coefficient  with  no  control 

- Incremental  normal  force  coefficient  due  to  control 

- Pitch  moment  coefficient  with  no  control 

- Incremental  pitch  moment  coefficient  due  to  control 


- Pitch  damping  coefficient 


' A block  diagram  of  the  missile  control  system  model  incorporated 
in  the  simulation  is  presented  in  Fig.  12.  The  system  includes  two 
accelerometers  aligned  with  the  Y],  and  axes  and  two  rate  gyros 
aligned  in  the  planes  of  the  wings.  A roll  channel  which  generates 
differential  displacements  of  opposite  surfaces  to  maintain  P*0  is 
assumed  to  operate  perfectly  and  is  not  modeled  in  the  simulation.  A 
vertical  g-bias  acceleration  is  normally  added  to  the  vertical  accelera- 
tion command  in  order  to  maintain  missile  trajectory  above  the  line  of 
sight  in  the  early  part  of  the  flight  for  terrain  clearance  (use  of  the 
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g-blas  in  this  study  Is  discussed  in  Section  VI).  A guidance  delay  is 
Incorporated  to  allow  for  safe  separation  from  the  carrier  aircraft. 
Thus,  no  accelerations  are  commanded  for  the  first  .7  seconds  and  while 
longitudinal  acceleration  is  greater  than  9.5  g’ s. 

Effective  pitch  and  yaw  surface  deflections  are  calculated.  These 
are  defined  as 


6C  - 


»(6i  + 63)  + (62  + <5») 

4 


(61  + 63)  + (62  + 64) 
4 


Jl 


JT 


(4) 


(5) 


where  61,  62*  5s  and  64  are  the  actual  control  surface  deflections. 
Reference  Frames 

Four  reference  frames  are  employed  in  the  Retran  Simulation:  an 

earth  fixed  frame,  a missile  body  frame,  an  aircraft  body  frame  and  an 
antenna  gimbal  frame.  The  missile  body  frame  is  fixed  in  the  missile 
airframe  and  is  obtained  from  the  earth  fixed  frame  by  a 3-2-1  rotation 
through  Euler  angles  <|),  6 and  <t> . This  rotation  is  illustrated  in 
Fig.  13.  Similarly,  the  aircraft  body  frame  is  obtained  by  a 3-2-1 
rotation  through  angles  0A  and  $A.  The  antenna  gimbal  frame  is 
obtained  from  the  aircraft  body  frame  by  a 3-2  rotation  through  angles 

and  0qA.  The  gimbal  frane  X axis  is  always  aligned  along  the  aircraft 
to  target  line  of  sight. 

Transformations  between  reference  frames  are  presented  in  Appendix  A. 
The  transformations  are  calculated  on  the  basis  of  instantaneous  angular 
relationships. 
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Fig.  13  3-2-1  Euler  Angle  Transformation 


Equations  of  Motion 

The  Retran  missile  is  symmetrical  about  its  longitudinal  axis 

I 

hence,  the  missile  body  axes  are  principal  axes  and  the  moments  of 
inertia  are  1^,  Iyy  - I2Z  and  I*y  ” 1^2  - Iy2  ■ 0.  Consequently, 
the  missile  equations  of  motion  are  (Ref  3:  2.24) 


m(U  - VR  + WQ)  - 

-mg  sin0  + Ax  + T 

(6) 

■(V  + OR  - WP)  - 

mg  sin<J>  cos0  + Ay 

(7) 

m(W  - UQ  + VP)  - 

mg  cos(j>  cos0  + Az 

(8) 

IxxP  - 

L 

(9) 

IyyQ  + (Ixx  “ IZZ)PR  ■ 

M 

(10) 

+ (Iyy  “ Ixx>PQ  ” 

N 

(ID 

These  equations  neglect  the  rate  of  change  of  mass  and  moments  of 
Inertia  because  these  rates  are  small  after  the  first  stage  burnout, 
while  no  flight  maneuvering  is  commanded  prior  to  first  stage  burnout. 

The  Retran  simulation  assumes  that  the  autopilot  roll  channel 
operates  perfectly  in  nulling  roll  rate.  Therefore,  it  is  assumed  that 
P“0  and  L=*0.  Incorporating  this  assumption,  the  equations  of  motion 


become: 

U - - g sin0  - QW  + RV  (12) 

m 

V * + g cos0  sln<J>  - RU  (13) 

m 

W m ^2.  + g cos0  cos<(>  + QU  (14) 

m 

P - 0 (15) 

Q - IM  + (Izz  - Ixx)pR]/lyy  (16) 

R - IN  + (1^  - Iyy)PQ]/l22  (17) 

The  aerodynamic  forces  and  moments  are  determined  by  the- equations: 

A*  ■ "1%  + ACC  + Cc  + ACof  ] q S (18) 

Ay  m ”1^2  " 8®n  (6C>]  q s (19) 

A*  " -I°N2  ” ^ 8«n  q S (20) 


M ■ iQnq  + Cjj^  - sgn  (6g)]  q Sd  (21) 

N - iCmq  - Cm  + Cj,  sgn  (6C))  q Sd  (22) 

M ^ 1 

where  d is  a reference  length:  d ■ 1 ft. 
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Aircraft  Maneuvering  and  Rotation  of  the  Line  of  Sight 

Following  missile  launch,  the  aircraft  is  assumed  to  fly  a nearly 
constant  arc  about  the  target,  maintaining  a horizontal  radar  gimbal 
angle  of  85°  ±7°.  The  aircraft  initially  turns  away  from  the  target 
until  the  commanded  gimbal  angle  is  reached,  then  turns  on  its  arcing 
maneuver.  The  aircraft  maintains  constant  airspeed  and  altitude.  Bank 
and  pitch  angles  are  neglected,  and  the  aircraft  is  assumed  to 
accomplish  the  turning  maneuver  through  heading  change  only. 

Two  parameters  which  are  utilized  in  the  guidance  algorithms  are 
the  angular  orientation  of  the  line  of  sight  and  its  rotation  rate. 
Angular  orientation  is  determined  by  first  calculating  the  radar  antenna 
gimbal  angles,  which  could  be  measured  directly  in  the  physical  mechani- 
zation of  the  system.  These  angles  are  calculated  by  determining  the 
position  of  the  target  relative  to  the  aircraft  in  the  aircraft  body 
frame.  This  position  is  located  by  the  position  vector 


Rat 


XAT  *A  + yAT  Ja  + ZAT  kA 


(23) 


From  the  geometry  of  this  relationship,  as  shown  in  Fig.  14,  the  gimbal 
angles  can  be  determined: 

ZAT 


0GA  " 8ln"1  [*AtJ 


♦ GA  " 8in 


-l  *AT 


(24) 

(25) 


[Rat]cos0ga 

Since  aircraft  bank  aod  pitch  angles  are  assumed  to  be  zero,  4>qA 
may  be  added  to  the  aircraft  heading  to  obtain  the  heading  of  the 
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Fig.  14  Gimbal  Angles 


line  of  sight  in  the  earth  fixed  frame: 

*LOS  " 4»A  + *GA  (26) 

To  determine  the  rotation  rate  of  the  line  of  sight,  the  antenna  gimbal 
angle  rates  are  determined  first,  as  these  rates  could  be  measured  in 
the  actual  mechanization.  The  velocity  of  the  target  relative  to  the 
aircraft  in  the  aircraft  body  frame  is 

ft  (RAT>  “ f^  ^(rATJ  + UGA  x RAt]  (27> 

Since  the  aircraft  is  flying  a nearly  constant  radius  arc  about  the 
target,  -jjj-  (n  ) can  be  neglected.  Now  the  angular  velocity  of  the 
gimbal  frame  relative  to  the  aircraft  can  be  determined  from  the 
relationship: 


J 


— * 


— — d ^ / p \ 

WGA  X RAT  “ dt  lRAT> 


(28) 


where  0)g  ■ 6qA  Jq  + kA.  T*ie  8imbal  angle  rates  are  therefore: 


*CA- 


®ATc°d® 


GA 


AT 


'GA 


R 


AT 


(29) 


(30) 


Again,  since  the  aircraft  is  in  level  flight,  the  rotation  rate 
of  the  line  of  sight,  defined  as  the  vertical  component  of  the  angular 
velocity  of  the  line  of  sight,  can  be  found  by  adding  the  horizontal 
glmbal  angle  rate  to  the  turning  rate  of  the  aircraft: 


u 


Z " *GA  + *A 


(31) 


Determination  of  Missile  Position  and  Velocity  Errors 

The  distance  of  the  missile  normal  to  the  line  of  sight  and  the 


velocity  component  of  the  missile  normal  to  the  line  of  sight  are 
determined  for  use  in  the  guidance  algorithms. 

The  distance  normal  to  the  line  of  sight  is  defined  by  the  YAK 
and  components  of  missile  position  relative  to  the  aircraft  in  the 
glmbal  frame,  as  shown  in  Fig.  15.  These  are  obtained  by  transformation 
from  the  earth  fixed  frame: 


Likewise,  the  velocity  of  the  missile  normal  to  the  line  of  sight 
Is  determined  by  finding  the  velocity  of  the  missile  relative  to  the 
aircraft  In  the  glmbal  frame: 


fit  G(RAM)  = dt _ (Rj- 


“GI  x RAM  ^3) 


The  vertical  and  horizontal  components  of  this  relative  velocity,  Y 
and  AZ,  are  then  the  velocities  normal  to  the  line  of  sight,  as  shown 


in  Fig.  15: 


XAM~  “ZYAM  + “YZAM 

- [L^]  + 0)zXam  - 0)xZam  (34) 

zam  " ‘Mam  + U)xyam 


Fig.  15  Missile  Position  and  Velocity  Errors 
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The  Guidance  Algorithms 

The  guidance  algorithms  which  were  incorporated  into  the  simulation 
for  the  launch  envelope  studies  are  presented  here.  Other  guidance 
algorithms  which  were  tested  and  rejected  are  discussed  in  Appendix  E. 
How  and  why  these  guidance  algorithms  were  chosen  is  revealed  in 
Section  VI. 

Three  guidance  algorithms  were  utilized:  one  to  determine 

acceleration  commands  in  the  vertical  plane  of  the  glmbal  frame  and 
two  to  determine  commands  in  the  horizontal  plane  of  the  gimbal  frame. 
These  commands  were  then  transformed  to  the  missile  body  frame  before 
being  input  to  the  missile  autopilot. 

I ( 

The  Vertical  Algorithm 

In  the  glmbal  frame  vertical  plane,  XgZg,  the  missile  is  always 
conmanded  to  intercept  and  track  the  line  of  sight.  Acceleration 
commands  are  proportional  to  the  missile's  vertical  distance  normal 

I 

to  the  line  of  sight,  AZ,  and  to  its  vertical  velocity  normal  to 
the  line  of  sight,  AZ.  The  algorithm  is  similar  to  the  algorithm 
for  a beam  riding  missile  discussed  by  J.  Clemow  in  reference  4 
(Ref  4:  48).  A g-bias  was  added  for  trajectory  shaping.  Note  that  the 
vertical  acceleration  ay  Is  commanded  in  the  -kg  direction.  The 
vertical  algorithm  is: 


( 
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(35) 


" 5 (kl  ~ ) AZ 

Ctg 


' 1 t 


tg 


The  term  kj  .HS_  is  analagous  to  *the  natural  frequency  wn  of  a damped 
fctg 

second  order  system,  where  kj  is  an  arbitrary  constant,  ttg  is  the 

tlme-to-go  to  target  Impact  and  tg  Is  an  arbitrary  reference  time  which 

corresponds  very  roughly  to  the  desired  line  of  sight  tracking  time. 

The  factor  in  effect  increases  the  guidance  system  gain  with 

Ctg 

decreasing  tlme-to-go  in  order  to  speed  system  response  as  the  target 
Is  approached.  £ represents  the  damping  ratio  of  the  system.  The 
values  chosen  for  these  parameters  are  discussed  in  sections  V and  VI. 


The  Horizontal  Algorithms 

In  the  early  portion  of  the  flight,  one  of  two  horizontal 
algorithms  is  employed,  the  choice  of  algorithm  depending  upon  missile 
position  relative  to  the  line  of  sight  after  the  first  two  seconds 
of  flight.  If  the  missile  remains  ahead  of  the  line  of-  sight  (4AY), 
then  the  missile  is  commanded  to  turn  to  a heading  parallel  to  the  line 
of  sight.  Since  the  line  of  sight  is  rotating  toward  the  missile, 
it  will  catch  up  to  the  missile,  at  which  point  the  guidance  system 
reverts  to  command-to-track  the  line  of  sight.  ' If  the  missile  falls 
behind  the  line  of  sight  (-AY)  then  the  command-to-track  algorithm 
is  employed  throughout  the  flight. 


The  command- to- heading  algorithm  calculates  a commanded  accelera- 
tion parallel  to  the  earth-fixed  XY  plane  and  perpendicular  to  the 
missile  longitudinal  axis.  The  acceleration  commanded  is  proportional 
to  the  angular  error  between  missile  heading  and  the  heading  of  the 
line  of  sight: 


ahr  - °t 


• TR-p 


Cr  is  a constant  chosen  to  convert  the  reference  time  TR  into  an 
equivalent  reference  range  and  Rjjj  is  the  mlssile-to-target  range, 
so  that  system  gain  in  this  case  increases  with  decreasing  range  to 
the  target.  C?  is  an  arbitrary  proportionality  constant.  Again, 
discussion  of  the  values  chosen  for  the  constants  is  presented  in 
sections  V and  VI.  a^  is  transformed  into  components  axg  and  a^g 
in  the  gimbal  frame. 

The  horizontal  command  to  track  algorithm  is  similar  to  the 
vertical  algorithm  except  that  in  place  of  a g-bias  term  a feed- 
forward rotational  bias  term  is  included.  This  term  provides  the 

r' 

acceleration  required  to  maintain  the  missile  on  the  rotating  line  of 
sight  with  zero  displacement  error.  The  algorithm  is 


«h  - -K  (kj  f*_)  Ay  . (k!  iS-)2  ay  + ( . *MlV>Z  , „„ 

£tg  * TvS — ’ (37) 


The  bias  term  was  determined  as  follows,  assuming  that  the  missile 
is  established  on  the  line  of  sight  in  the  horizontal  (XgYg)  plane  as 
in  Fig.  16.  The  lateral  acceleration  required  to  maintain  the  rotating 
line  of  sight  is  (Ref  5:  4.d.(l)-3): 


I* 
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Fig.  16  Lateral  Acceleration  Required  to  Maintain  Rotating  Line  of 
Sight 

Lo  " (*M]Wz  + 2rM3wz*  cos0  - (®MT  ” 8^n0  (38) 

• • 
substituting  cosip  - ^MT  an<j  sin0  « ^SfT^Z 

VM  VM 

T _■  %rRwra)z  + 2®mt2u>z  ” ^Mr^rr^z  + ^fr^z*  ,„„x 

L0  ■ (39) 

VM 

Since  the  aircraft  is  flying  a nearly  constant  speed  arc  about  the 

• • 

target,  we  can  neglect  If  we  also  assume  <|>  and  4>  are  small,  such 

* • ••  ( • 

that  R^ff  » Vjf  and  R^  * VM,  the  acceleration  required  becomes 

Lo  - -«OT  «z  (*0) 

M 


Transition  from  the  command- to-heading  algorithm  to  the  comma nd- 
to-track  algorithm  is  made  at  the  point  where  the  damping  term  in  the 
command-to-track  algorithm  reverses  the  direction  of  the  acceleration 
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A first  order  gradient  algorithm  employing  penalty  functions  was 
developed  in  an  attempt  to  optimize  the  missile  trajectory  for  maximum 
range.  The  algorithm  is  similar  to  the  gradient  algorithm  employed  by 
Anderson  and  Smith  in  reference  £ (Ref  6:  230),  which  follows  the 

< « 

technique  described  by  Kelley  in  reference  7 (Ref  7:  205-254). 

The  algorithm,  as  applied  to  the  SAR-Retran  problem,  is  designed 
to  determine  the  optimal  set  of  acceleration  conmands  a^  and  ay  which 
produce  the  trajectory  yielding  maximum  range.  The  steps  used  in  the 
algorithm  are  as  follows: 

1.  The  state  equations  X(t)  are  integrated  forward,  using 
commands  calculated  by  a guidance  algorithm  similar  to  the  one  employed 
in  the  Retran  simulation.  These  commands  are  stored.  Forward 
integration  is  terminated  when  the  missile  intercepts  the  line  of 
sight  (within  50  feet).  This  yields  X(t)  and  an  estimate  of  the 
final  time  tf. 

2.  The  payoff  function  J,  including  penalties,  is  calculated. 

3.  The  adjoint  vector  X(tf)  is  found  from 


I 


5*  The  following  differential  equation 


X1  - - dn/ix. 


(46) 


la  integrated  backward  to  tQ,  using  the  states  obtained  from  the  forward 
Integration.  The  controls  determined  from  step  1 or  the  previous  step  6 


are  used . Simultaneously,  the  gradients 
calculated  and  stored. 


J*{ 


t)  and  (t)  are 

O 3v 


6.  The  incremental  corrections  to  the  acceleration  commands  are 


6«h(t)  - -a  (t)  and  6ay(t)  - -a  ^~~(t).  With 


( 


•h  ■ + 

(47) 

\ ' *v(c)old  + 

(48) 

an  alpha  search  is  performed  to  minimize  J with  respect  to  a.  The  new 
a^Ct)  and  ay(t)  are  stored. 

7.  Using  the  new  controls  resulting  from  step  6,  the  state  equations 
are  again  Integrated  forward  to  the  stopping  condition  and  steps  2 through 
6 are  repeated.  This  process  is  carried  on  until  J is  minimized  within 
the  desired  tolerance. 

1 

. 

The  Missile  Model 

To  facilitate  the  use  of  the  gradient  technique,  a simplified 
missile  model  was  developed.  Rotational  dynamics  of  the  missile  were 
Ignored.  The  missile  was  assumed  to  maintain  zero  bank  angle  throughout 
the  flight.  Aerodynamic  coefficient  curves  were  approximated  by  continuous 
functions.  Small  angles  were  assumed  for  a and  $ and  control  deflections 
were  neglected,  except  in  the  determination  of  ACC< 

•j 

J 

29  . . 
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Equations  Of  Motion. 

Incorporating  the  above  simplifications,  the  vector  equation  of 
the  missile  Is 


mVM  “ AxB  + + AzB  + T + ®8 


which  yields  the  scalar  equations 


Vm-£  I(Ax+T)  cosacosS  + Ay  sing  + ^ sinotcosB]  - g siny  (50) 
]J>VjjCosa  = ^ [-(Ajj+T)  cosasinB  + Ay  cos 3 - sinctsinB]  (51) 

-yVM  - ^ [-(Ajj+T)  sino  + cosa]  + g cosy  (52) 


Aerodynamic  Forces 

The  aerodynamic  forces  are  determined  from  the  relations 


Ajj  - -Sq  [Cqq  + Cg  + ACc  + ACDf  J (53) 

Ay  - -Sq  [Ch2  “ sgn  (5g)]  (54) 

A*  - “Sq  [Cu2  “ CNi  sgn  (6*)]  (55) 

where  S is  missile  cross-sectional  area. 

Continuous  functions  approximating  the  aerodynamic  coefficients  were 
obtained  by  a least-squares  fit  to  the  coefficient  curves  (using 
AFITSUBROUTINE  PLSQ) . These  curve  fits  are  illustrated  in  Figs.  18 
through  21.  Note  that  Cg  is  assumed  to  be  zero  (valid  for  M - .8), 

Cg  is  assumed  to  be  zero  (valid  for  6*  0)  and  the  ACg  curve  selected 


3< 
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is  the  curve  for  6 ■ 4°,  which  is  approximately  the  average  value  of 
control  deflection  observed  in  simulated  flights.  The  following  functions 
were  obtained: 

CDq  - -51.658M7  + 380.974M6  - 1151. 244M5  + 1844. 332M4  (56) 

-1689. 781M3  + 885.882M2  - 246.298M  + 28.250 

ACC  - -(2.851  x 10~5)  Oj.2  - (1.829  x lO”3)  Oj.  + 8.571  x 10~3  (57) 

where  Oj.  ■ /a2  +0*  (deg) 

AcDf  * ~ (7.406  x 10”12)(-Z)2  - (1.079  x 10~7)(-Z)  - 8.878  x 10“s  (58) 
Cjj^  " 20.1  a (59) 

where  a is  expressed  in  radians. 

Thrust  and  Mass. 

Thrust  is  obtained  directly  from  the  thrust  schedule  in  section  II. 

t 

Time  of  Burn  (sec)  Thrust  (lb.) 

0.0-0. 6 9500. 

0.6-4. 5 2000. 

Mass  is  represented  as  a linear  function  between  the  stage  burnout 
masses. 

M - Mo  - Cxt  + Ci  <t-tx>  - C2<t-t1>  + C2<t-t2>  (60) 


where  the  stage  burnout  masses  are 


time  (sec) 
mass  (slugs) 


0.0 


0.6 


4.5 


(Note  that  <t-tQ>  indicates  a singularity  function  which  does  not  become 
effective' until  t>tQ) . 


Control. 


The  Independent  controls  in  this  algorithm  are  the  acceleration 
connands  in  the  missile  body  axes,  aft  and  a^.  In  the  simplified  missile 
model,  it  is  assumed  that  the  autopilot  delivers  sufficient  control 
deflection  to  obtain  the  comnanded  accelerations,  within  the  autopilot 
limit  of  ±20  g's.  The  lateral  equations  of  motion  can  then  be  written  as 


ah  “ ~ I-(Ajt+T)cosotsin6  + Ay  cos8  - Az  slnasinS] 


Sv  = ^ I (Ajj+T)  sina  - cosa]  -g  cosy 
Assuming  small  angles  for  a and  8 


sin  a - a 


sin  8 * 8 


cos  a ■ 1 


cos  8 m 1 


and  neglecting  the  terms  containing  a and  8 as  compared  to  the  other 
terms,  the  equations' of  motion  can  be  approximated  as 


ah  - — i 
m 


: Fz  _ 


g cosY  - 


-g  cosY 


Assuming  that  the  a and  8 necessary  to  attain  the  commanded  accelera- 
tions may  be  obtained  instantaneously  (since  the  change  at  each  time 

step  is  normally  very  small),  a and  6 nay  be  approximated  by  the  following 
relationship: 
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(32.2  ay  + g cosy)  m 
20.1  Sq 

.(32.2  %)  m 
20.1  Sa 


a and  $ are  restricted  to  a maximum  of  ±20  degrees,  which  Is  consistent 
with  the  maximum  values  observed  in  simulated  flights. 


Atmosphere 

The  atmosphere  is  modeled  according  to  the  1959  ARDC  standard 
atmosphere. 


Tenperature:  t ” -aH  + 


Density: 


: Mach  No: 


- .003566Z  + 518.69  °R 
P " C2T  €«</ 

- 6.3806  x 10"15  T4,26 

ft3 


.0204  -M- 

ST 


The  State  Equations. 


The  state  equations  consist  of  the  three  equations  of  motion 

" m I (A^-fT)  cosacosB  + Ay  sinB  + A®  sinOtcosB]  -g  sinY  (70) 

■ m^cosY — I-(Ax+T)cosasin8  + Ay  cos  8 - A2sinasin8]  (71) 
Y " ' I-(Ax+T)sina  + Az  cosa  ] - g—  cosY  (72) 


and  the  three  kinematic  relationships 
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X • Vjj  CO8YCO8* 


(73) 


Y • VM  cosYeini|/ 
Z - -VM  slnY 


(74) 

(75) 


The  forward  integration  of  the  state  equations  follows  this 
sequence : 

1.  Flight  is  initialized  at  t * 1.0  seconds  to  avoid  complications 
introduced  by  the  safe  separation  delay.  Initial  parameter  values  were 
obtained  from  the  Retran  simulation  at  the  end  of  the  first  second  of 
flight  for  each  of  the  launch  positions  investigated.  % and  ay  are 
initialized  at  zero. 

2.  Atmospheric,  mass  and  thrust  quantities  are  determined. 

3.  a and  8 are  determined. 

4.  Aerodynamic  forces  are  calculated. 

5.  Missile  equations  of  motion  are  integrated. 

6.  Kinematic  equations  are  integrated. 

7.  The  state  values  are  stored 

8.  A test  is  made  to  determine  if  the  stopping  condition  has  been 
reached.  The  stopping  condition  is  defined  as  (X-XL0S)2  + (Y-Y^g) 2 + 
(Z-Zlos) 2 < (50) 2 where  XL0S,  YL0S  and  ZjXjg  are  the  coordinates  of  the 
point  on  the  line  of  sight  at  the  same  range  from  the  target  as  the 
missile. 

9.  If  the  stopping  condition  has  not  been  reached,  then  for  the 
initial  forward  integration,  new  guidance  commands  are  calculated  using 
the  guidance  algorithms.  For  subsequent  forward  integrations,  the  stored 
values  of  at,  and  ay  are  used* 


! 
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10.  Steps  2 through  9 are  repeated  until  the  stopping  condition  Is 
achieved . 

Payoff  Function  and  Penalties 

Missile  range  Is  primarily  limited  by  the  necessity  of  achieving  a 
minimum  value  of  terminal  airspeed  at  target  Impact.  Thus,  any 
trajectory  that  maximizes  range  must  maximize  this  terminal  airspeed. 
Consequently,  the  negative  terminal  airspeed,  -V^,,  was  chosen  as  the  basic 
payoff  function;  the  algorithm  will  attempt  to  minimize  -VT,  thereby 
maximizing  VT.  However,  in  order  to  avoid  complications  imposed  by  line 
of  sight  tracking  in  this  program,  the  flight  is  terminated  when  the 
line  of  sight  is  intercepted  and  tracking  of  the  line  of  sight  to 
the  target  is  assumed.  Prom  the  results  of  the  Retran  simulation,  it 
appears  safe  to  assume  that  the  further  decrease  in  airspeed  while 
tracking  the  line  of  sight  can  be  approximated  as  a function  of  range  from 
the  target;  i.e.,  the  terminal  velocity  can  be  represented  as 

VT  " VMf  “ (k01  Rf2  + *02  Rf>  (76) 

where  is  the  missile  airspeed  at  line  of  sight  intercept  and  Rf  is 
the  range  to  the  target  at  the  point  of  interception.  The  basic 
payoff  function  is  therefore 

J - -VMf  + koX  Rf2  + koj  Rf  (77) 

The  values  of  and  kQ2  were  obtained  by  a least  squares  fit  of  the 
velocity  versus  range  profile  obtained  from  the  Retran  simulation  for  a 
missile  tracking  the  line  of  sight.  Profiles  were  obtained  for  launch 
altitudes  at  35,000  feet  and  20,000  feet  to  account  for  the  change  in 
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Inclination  of  the  line  of  sight  with  varying  altitude.  The  curve  fits 
to  these  profiles  are  illustrated  in  Figs.  22  and  23.  The  values  of 
and  ligj  obtained  were,  for  the  35,000  foot  launch  altitude: 

kO!  - -2,700  x 10"8 

koj  “ 3.395  x 10"3 

and  for  the  20,000  foot  launch  altitude: 

kgj  - 2.815  x l<f8 

1^02  - 3.031  x 1(T3 

Note  that  the  function  representing  the  decrease  of  airspeed  on  the  line 
of  sight  is  analagous  to  a penalty  function  which  penalizes  the  missile 
for  intercepting  the  line  of  sight  too  early  if  lower  drag  could  other- 
wise be  achieved. 

In  order  to  track  the  line  of  sight  following  intercept,  the 
missile  must  be  generally  aligned  with  the  line  of  sight;  otherwise, 
high  g maneuvers  are  required  to  align  the  missile,  resulting  in 
increased  drag.  To  account  for  this  extra  drag,  penalty  functions  were 
Included  to  penalize  the  missile  for  heading  errors  at  line  of  sight 
intercept.  The  heading  tolerances  were  established  as 

0 i ^ 

0 i It  * SL0S  ± l0‘ 

Likewise,  a penalty  function  was  added  to  penalize  the  missile  for 
intercepting  the  line  of  sight  at  too  short  a range  to  provide  adequate 


tracking  time.  The  range  tolerance  was  established  as 
Rf  > 16,000  ft 

to  allow  a minimum  of  20  seconds  tracking  time  at  an  average  missile 
airspeed  of  700  to  800  ft/sec. 

The  resulting  payoff  function.  Including  penalties,  is 

2 I 

J - -VMf  + kojRf  + ko2Rf  + 2 k1(i|>f-ipL0S)2 

+ I k2  (yreL0S)2  + I k3  (Rf-16,000.)2  (78) 


The  values  of  kj,  k2,  and  k3  were  chosen  through  a sensitivity 
analysis  which  is  discussed  in  section  V. 

Adjoint  Vector.  Hamiltonian  and  Backward  Integration 
The  adjoint  vector  is  defined  as 

V^tf)  - dJ/dX  tf  (79) 

where  dj/dx  - [dj/dV  dJ/dy  dj/dx  2>J/dY  dj/dz] 

This  vector  is  calculated  using  the  final  conditions  obtained  from 
forward  integration  of  the  state  equations. 

The  Hamiltonian  is  defined  as 

H - X1  • X (80) 

The  differential  equation 

t1  - -dH/dX  (81) 


where 


dH/dX  - [dH/dV  dH/Jt|>  dH$  y dH/dX  dH/dY  dH/dZ] 


is  Integrated  backward,  along  with  the  state  equations,  to  obtain  A(t) 
and  X(t).  The  X(t)  and  the  controls  ah(t)  and  ay(t)  which  were  calculated 
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and  scored  during  the  forward  Integration  pf  the  state  equations  are 
usedln.thls  calculation. 

Having  found  X(t)  and  x(t) , the  gradients  (t)'  and  (t) 
can  be  determined  and  stored. 


Alpha  Search 

The  adjustments  to  be  made  in  the  controls  in  order  to  improve  the 


solution  are  $%(t)  - 


(t)  and  day  (t) 


(t).  In  order  to 


determine  the  value  of  a which  will  minimize  the  solution  at  this 
iteration,  a one-dimensional  alpha  search  is  performed  to  minimize  J 
with  respect  to  a.  Having  determined  a,  the  new  acceleration  commands 
are  calculated  and  stored: 


«h(t)  - 


6ah(t) 


- «y(t)old  + 6ay(t) 

Using  these  new  camnands,  the  state  equations  are  again  integrated 
forward  to  obtain  a new  solution  and  the  iteration  is  repeated. 


Algorithm  Termination 

The  algorithm  terminates  when  further  iterations  will  make 
negligible  reductions  in  the  payoff.  This  point  is  determined  when 
successive  payoffs  are  within  a chosen  tolerance  of  one  another;  a 
tolerance  of  .01Z  was  employed. 


( 
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Program  Retran  Is  the  Ret ran  missile  simulation  program.  It  is  a 

! modification  of  a Maverick  missile  five  degree  of  freedom  simulation 

program  which  was  developed  by  Mr.  Merril  Habbe  of  the  Air  Force  j 

Avionics  Laboratory.  The  Maverick  missile  data  and  parameters  were 

] 

replaced  by  the  Retran  missile  parameters,  and  the  algorithms  necessary 
to  represent  the  launching  aircraft  and  to  incorporate  SAR-Retran 

• ' ' * “ i 

guidance  were  Included.  A Calcomp  plotting  routine  was  added.  The 

. 

program  is  written  in  FORTRAN  Extended,  Version  4,  for  use  on  the  CDC 

i 

6600-CYBER  74  computer  system.  A source  program  is  attached  as 
• - i 

Appendix  B. 


The  program  employs  simple  Euler  numerical  integration.  Time 
steps  of  less  than  .05  seconds  were  found  to  yield  reasonable  accuracy 
(within  5Z  of  results  using  a very  small  time  step  of  .005  seconds); 
the  time  step  is  reduced  by  a factor  of  10  in  the  final  10  seconds  of 
flight  for  increased  accuracy  in  the  critical  terminal  phase.  The 
program  contains  an  abbreviated  ARDC  1959  atmospheric  table  and 
tables  of  aerodynamic  coefficients  for  the  missile. 

The  coefficient  tables  were  obtained  from  the  coefficient  curves 
shorn  in  Figs.  3 through  11.  Thrust  and  mass  parameters  are  also 
tabulsrlzed.  A look-up  routine  employing  linear  interpolation  is  used 
to  extract  values  from  these  tables. 

Flight  parameters  are  initialized  at  t0  - 0.  Data  is  input  on 
three  cards  (required  input  is  explained  at  the  beginning  of  the 
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program  listing  in  Appendix  B)  and  Includes  aircraft  position  coordi- 
nates, Euler  angles  and  velocity  components  at  launch,  together  with 
target  position  and  velocity  (however,  the  program  as  listed  assumes  a 
stationary  target  and  Ignores  input  target  velocity).  The  missile 
parameters  are  initialized  to  the  aircraft  parameters  at  launch. 

The  program  terminates  when  the  missile  passes  through  target 
altitude  or  when  a specified  maximum  flight  time  is  exceeded.  Output 
data  is  printed  during  the  flight;  these  data  are  explained  at  the 
beginning  of  the  output  listing,  a sample  of  which  is  contained  at  the 
end  of  the  program  listing  in  Appendix  B.  By  setting  the  proper  integer 
flag,  aircraft  and  missile  position  coordinates  are  punched  on  data 
cards.  These  can  be  utilized  with  a separate  plotting  routine  to 
produce  multiple  trajectory  plots  such  as  those  presented  in  section  V. 
Setting  another  integer  flag  will  produce  the  following  Calcomp  plots: 

1.  Plan  view  of  trajectory 

2.  Vertical  profile  of  trajectory 

3.  Acceleration  commands  vs.  time 

4.  a and  8 vs.  time 

5.  Deceleration  due  to  aerodynamic  drag  vs.  time 

6.  Retransmitted  frequency  error  vs.  time 

7.  Missile  distance  normal  to  the  line  of  sight  for  the  final 
20  seconds  of  flight 

Examples  of  these  plots  may  be  seen  in  Section  V and  Appendix  D. 

Multiple  flights  may  be  simulated  in  one  input  run  by  including 
the  appropriate  number  of  data  cards. 
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Program  Optral 


C 

Program  Optraj  is  the  first  order  gradient  optimisation  program. 

This  program  is  also  written  in  FORTRAN  Extended,  Version  4,  for  use 
on  the  CDC  6600— CYBER  74  computer  system.  A source  program  listing 
is  attached  as  Appendix  C. 

This  program  utilizes  simple  Euler  numerical  integration.  Time 
steps  of  .5  seconds  were  found  to  yield  satisfactory  results.  Flight 
parameters  were  initialized  at  tQ  ■ 1.0  seconds.  Data  is  input  on  6 
cards  (required  Input  is  explained  at  the  beginning  of  the  program 
listing  in  Appendix  C)  and  includes  aircraft  and  missile  position  coordinates, 
Euler  angles  and  velocity  components  at  t0»  together  with  initial  values  of  a 
and  The  penalty  function  coefficients  are  also  input. 

The  program  terminates  when  successive  payoff  values  lie  within 
a designated  tolerance  of  each  other.  The  program  may  also  terminate 
if  an  initial  solution  is  not  obtained  or  if  improvement  of  the  Initial 
solution  cannot  be  made.  A message  is  printed  for  these  abnormal 
terminations.  The  following  output  data  is  printed  during  the  optimiza- 
tion process: 

1.  The  missile  velocity  and  position,  velocity  and  distance  normal 
to  the  line  of  sight  and  acceleration  conmands  (in  missile  body  axes) 
are  printed  for  the  initial  trajectory  and  again  for  the  optimized 
trajectory.  Only  final  values  are  printed  for  the  other  forward 
integrations  of  the  state  equations. 

2,  The  value  of  the  payoff  function  and  the  calculated  terminal 
velocity  are  printed  for  each  forward  integration.  The  angular  errors 
and  range  error  at  line  of  sight  interception  are  also  printed. 
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^ H ^ H 

3.  The  values  of^jjj-  and^£-  are  Printed  flDr  each  backward 


Integration. 


4.  The  alpha  values  and  payoff  values  are  printed  for  each  alpha 


search. 


Printed  output  for  the  Optraj  program  can  be  seen  in  Appendix  F. 

A nunber  of  Calcomp  plots  are  drawn  at  the  completion  of  the 
optimization  process.  These  plots  present  the  initial  curves  and 
optimized  curves  for  the  following  parameters: 


1.  Trajectory  plan  view 

2.  Trajectory  vertical  profile 

3.  *h<j 


5.  a 


6.  $ 


7.  V* 


8.  AT 

9.  AZ 


Examples  of  these  plots  may  be  seen  in  section  V. 

Only  one  trajectory  may  be  optimized  on  an  input  run. 


Validity  of  Program  Results 

Ret ran  Simulation.  Because  the  Retran  missile  is  a hypothetical 
vehicle,  no  flight  data  is  available  to  validate  the  results  of  the 
simulation.  However,  the  Maverick  missile  5 degree  of  freedom  simulation 
upon  which  this  program  is  based  has  been  successfully  validated  against 
flight  test  data.  While  the  performance  of  the  Retran  missile  differs 


I 

t 


45 


n ■ 

from  that  of  the  Maverick,  the  simulation  results  appear  to  be  reasonable 
In  the  light  of  the  Maverick  simulation  results.  Unfortunately,  the 
Maverick  simulation  results  are  classified;  therefore,  a direct 
comparison  cannot  be  presented  in  this  report. 

Optra!  Optimization.  An  attempt  was  made  to  validate  the  performance 
of  the  Optraj  program  by  obtaining  the  same  optimal  solution  from  two 
videly  differing  initial  solutions.  However,  it  became  apparent  from 
these  attonpts  that  numerous,  relatively  closely  spaced  local  minimums 
were  encountered  in  the  payoff  function.  Since  the  first  order 
gradient  technique  will  converge  on  the  local  minimum  within  the  same 
range  as  the  initial  solution,  the  program  converged  on  different 
optimal  solutions  for  widely  differing  initial  solutions.  In  order 
to  obtain  convergence  on  the  same  optimal  solution,  the  initial 
solutions  could  differ  only  slightly. 

Because  the  trajectory  changes  made  by  the  optimization  algorithm 
were  too  slight  to  duplicate  in  the  Retran  simulation,  the  effects 
of  the  optimization  could  not  be  validated  in  the  simulation. 


V.  Results  . 


Guidance  Algorithm  Selection 

Several  guidance  algorithms  were  tested  in  the  Retran  simulation 
In  order  to  select  the  algorithms  most  suitable  for  employment  of  the 
missile  with  SAR-Retran  guidance.  The  algorithms  which  were  rejected 
are  discussed  in  Appendix  E,  while  discussion  of  the  algorithms 
presented  here  may  be  found  in  section  II.  The  algorithms  selected 
were : 

vertical  command-to- track: 
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horizontal  conmand- to-heading: 


ahT  - &r  Wco8  - i|>) 


CiC  tR 


L % J 


horizontal  command-to-track: 
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The  coefficient  values  selected  for  these  algorithms  were 
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Horizontal  and  vertical  trajectories  produced  by  these  algorithms 
are  presented  in  Figs.  24  and  25  for  a long  range  medium  squint  angle 
launch  (squint  angle  is  the  angle  between  the  aircraft  velocity 
vector  and  the  line  of  sight). 

Performance  parameters  which  were  considered  important  in  the 
selection  of  the  guidance  algorithms  were  the  schedule  of  accelera- 
tions commanded  by  the  algorithms,  the  time  histories  of  angle  of 
attack  a and  sideslip  angle  3,  the  profile  of  deceleration  due  to 
aerodynamic  drag  and  the  missile  distance  errors  in  the  terminal  phase. 
The  retransmitted  frequency  error  was  also  considered,  as  maintenance 
of  this  error  below  8 hz  was  the  criteria  for  line  of  sight  tracking. 
These  performance  parameters  are  presented  in  Figs.  26  thru  30  for 
the  long  range  medium  squint  angle  launch.  Parameter  plots  for 
launches  from  different  altitudes  and  ranges  may  be  seen  in  Appendix  D. 

Trajectory  Optimization 

An  attempt  was  made  to  optimize  the  acceleration  command  schedule 
to  produce  a trajectory  yielding  maximum  range.  An  initial  squint 
angle  of  45  degrees  was  chosen  for  the  optimization,  with  launches 
made  from  altitudes  of  35,000  ft.  and  20,000  ft.  Initial  trajectories 
were  obtained  using  essentially  the  same  guidance  algorithms  employed 
in  the  Retran  simulation,  although  the  values  of  Cj,  tR,  C,  and  gjj^as 
were  varied  slightly  in  order  to  achieve  trajectory  shapes  similar 
to  those  obtained  from  the  simulation. 
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A sensitivity  analysis  determined  that  penalty  function  values 
of  ■ 1000.,  V.2  = 1000.,  and  k3  = 1.0  x 10-6  would  yield  satisfactory 
solutions. 

The  first  order  gradient  algorithm  directed  only  small  changes 
in  the  acceleration  command  schedule  before  reaching  a local  minimum 
of  the  payoff  function.  The  resulting  trajectory  changes  were 
negligible,  as  were  changes  in  the  performance  parameters.  Plots  of 
the  initial  and  optimized  trajectories  and  performance  parameters 
for  the  35,000  ft.  launch  are  presented  in  Figs.  31-39.  For  this 
launch,  the  payoff  function  was  improved  from  J = -731.88  to 
J *=  -731.89,  with  terminal  velocities  of  771.78  for  the  initial 
trajectory  and  771.79  for  the  optimized  trajectory.  The  trajectories 
and  performance  parameters  for  the  20,000  ft.  launch  are  presented  in 
Figs.  40-48;  here,  the  payoff  function  went  from  J * -307.49  to 
J ■ -307.57,  with  terminal  velocities  of  452.89  and  452.96. 

Because  the  trajectory  and  performance  parameter  changes  were 
too  small  to  be  distinguishable  on  the  Calcomp  plots,  the  output 
listings  for  these  two  program  runs  are  included  in  Appendix  F. 

Launch  Envelope  Determination 

Using  the  algorithms  selected  through  the  simulation  Investigation, 
the  launch  envelope  for  the  Retran  missile  was  determined.  The  criteria 
used  in  determining  a successful  launch  were: 

1.  Minimum  terminal  airspeed  of  700  ft/sec. 

2.  Minimum  tracking  time  on  the  line  of  sight  (Af<j  <_  8 hz)  of 
20  seconds. 
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Simulated  flights  were  made  to  determine  the  maximum  and  minimum 
ranges  obtainable  at  a given  initial  squint  angle.  Squint  angles  of 


1 


15°,  30°,  45°,  60°,  75°,  and  90°  were  Investigated,  with  launches 
made  from  altitudes  of  5,000  ft,  20,000  ft  and  35,000  ft.  Thus, 
three  sets  of  planview  launch  envelopes  were  obtained  - one  for  each 
launch  altitude.  The  launch  envelopes  are  shown  in  Figs.  49-51  as  they 
would  appear  relative  to  the  aircraft;  these  figures  indicate  that 
any  target  lying  within  the  envelope  corresponding  to  aircraft  altitude 
could  be  successfully  struck  by  a Retran  missile.  A different  view 
of  the  launch  envelope  is  presented  through  Figs.  52-57.  These  show 
the  planview  trajectories  for  maximum  and  minimum  range  launches  at 
each  of  the  three  launch  altitudes,  thus  together  they  indicate  the 
size  and  shape  of  one-half  of  the  launch  envelope  relative  to  the 


target . 

Note  that  the  minimum  squint  angle  of  15  degrees  was  chosen  as 
a practical  minimum  value  at  which  target  acquisition  and  tracking 
would  be  accomplished  by  the  SAR-Retran  system.  This  is  not 
necessarily  a missile  limitation,  although  no  launches  were  attempted 
from  lower  squint  angles.  Likewise,  90  degrees  was  the  maximum 
practical  squint  angle  obtainable  with  the  aircraft  gimbal  system. 
35,000  ft  was  taken  as  the  maximum  operational  altitude  for  a loaded 
fighter  aircraft,  while  5000  ft  appears  to  be  the  minimum  altitude 
from  which  a useful  launch  envelope  can  be  obtained  due  to  missile 
range  limitations. 
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Retran  Missile  Launch  Envelopes  Relative  to  Aircraft  - 20,000  ft  Launch  Altitude 
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Fig.  53  Trajactory  Plaa  Flaws,  35,000  ft 
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Trajectory  Plan  Views , 5,000  ft  Minimum  Range  Launches 


VI.  Discussion  of  Results 


Guidance  Algorithm  Selection 

The  guidance  algorithms  were  selected  on  the  basis  of  missile 
performance  obtained  in  the  Retran  simulation.  Several  algorithms 
were  studied,  with  primary  emphasis  placed  on  obtaining  maximum  range. 
Preliminary  investigations  had  been  made  by  the  Air  Force  Avionics 
Laboratory  utilizing  algorithms  that  were  similar  to  the  comraand-to- 
track  algorithms.  At  lew  initial  squint  angles  (below  30  degrees), 
where  the  missile  initially  falls  behind  the  line  of  sight,  these 
algorithms  produced  a long  curving  flight  path,  as  can  be  seen  in 
Fig.  97  . The  first  algorithms  tested  in  this  study  produced  a straight 
line  path  to  the  line  of  sight  intercept,  then  a relatively  high-g 
(4  g's)  turn  onto  the  line  of  sight  (these  algorithms  are  discussed 
in  Appendix  E).  Flight  time  and  flight  distance  were  reduced;  however, 
the  drag  penalty  incurred  during  the  intercept  turn  resulted  in  lower 
terminal  velocities  than  had  been  obtained  with  the  cammand-to-track 
algorithms.  The  algorithms  were  modified  to  reduce  the  turning  angle 
required  at  line  of  sight  intercept,  but  terminal  velocities  were 
still  low. 

The  desirable  feature  of  the  command-to-track  algorithm  for  low 
to  medium  squint  angle  launches  is  that,  with  proper  selection  of 
frequency  coefficient  and  damping  ratio,  a fairly  uniform  schedule  of 
l(W-g  commands  can  be  achieved,  keeping  angle  of  attack  and  sideslip 
angle  l<w  throughout  the  flight.  Consequently,  a smooth  deceleration 
due  to  drag  was  maintained,  with  no  sharp  airspeed  losses  (see  the 
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parameter  plots  In  Appendix  D).  However,  at  high  initial  squint  angles 
(greater  than  45  degrees)  where  the  missile  accelerates  out  ahead  of 
the  line  of  sight,  the  conmand- to- track  algorithm  commands  large 
accelerations  (6-9  g's)  in  the  initial  turn,  then  directs  a fairly 
large  intercept  angle  back  to  the  line  of  sight  which  necessitates  a 
high-g  (10-15  g's)  turn  onto  the  line  of  sight.  At  very  high  squint 
angles,  the  accelerations  required  can  exceed  missile  maneuvering 
capability. 

To  achieve  satisfactory  trajectories  from  high  squint  angles, 
the  command- to-heading  algorithm  was  added.  With  proper  choice  of  the 
turning  coefficient  Gp,  the  schedule  of  acceleration  commands  can  be 
kept  uniformly  low.  The  missile  is  conmanded  to  turn  to  parallel  the 
line  of  sight,  with  line  of  sight  rotation  toward  the  missile  providing 
the  closure  rate.  Since  this  closure  rate  is  relatively  low,  a smooth 
transition  to  the  command-to-track  algorithm  can  be  made,  with  no  large 
accelerations  required. 

Launches  from  minimum  ranges  present  a different  problem.  Here, 
large  accelerations  are  required  to  achieve  early  Interception  to  the 
line  of  sight  so  that  minimum  tracking  time  can  be  attained.  These 
large  accelerations  are  obtained  by  scheduling  the  gain  of  the  conand- 
to-headlng  algorithm  proportional  to  the  inverse  square  of  misslle-to- 
target  range.  Likewise,  the  frequency  of  the  conmand- to- track 
algorithm  Is  calculated  proportional  to  the  inverse  of  time-to-go  to 
target  Impact,  thus  speeding  response  at  shorter  ranges.  Range  was 
chosen  for  gain  scheduling  of  the  command-to-heading  algorithm  because 
the  large  angle  between  missile  velocity  and  the  line  of  sight  during 


75 


V / 


initial  employment  of  this  algorithm  makes  tlme-to-go  an  unreliable 
indicator. 

The  coefficient  values  t^,  K,  and  Cj  were  chosen  to  yeild  a 
fairly  uniform  schedule  of  acceleration  comnands  throughout  the 
flight.  The  value  of  the  range  conversion  coefficient,  C^,  was  chosen 
as  1000  ft/sec,  which  corresponds  roughly  to  the  average  missile 
airspeed  during  the  initial  turning  portion  of  the  flight. 

The  damping  ratio  c;  = .5  was  the  highest  damping  ratio  that  did 
not  produce  large  acceleration  commands  during  the  transition  to  the 
horizontal  command -to- track  algorithm.  This  same  damping  ratio  was 
utilized  in  the  vertical  algorithm  and  produced  smooth  convergence 
to  the  line  of  sight  after  only  one  overshoot  in  both  the  vertical  and 
horizontal  planes  (this  can  be  seen  from  the  distance  error  plots  in 
Appendix  D).  This  smooth  convergence  was  obtained  by  not  allowing 
the  algorithm  frequency  to  increase  beyond  5 seconds  time-to-go;  if 
frequency  were  allowed  to  Increase  without  bound,  the  missile 
trajectory  became  erratic  in  the  final  seconds  of  flight.  In  the 
actual  guidance  system,  a larger  damping  ratio  would  probably  be 
desired  in  the  terminal  phase,  with  no  limitation  on  maximum  frequency 
in  the  algorithm.  However,  it  was  felt  that  the  terminal  accuracy 
obtained  with  £ ■ .5  was  sufficient  for  the  purpose  of  this  study. 

A g-bias  of  lg  produced  the  desired  trajectory  shaping  for  terrain 
clearance,  as  can  be  seen  in  the  vertical  trajectory  profiles. 

Autopilot  Compatibility 

One  of  the  goals  of  this  study  was  to  determine  the  stability  of 
the  guidance  loop  with  the  roll  stabilized  autopilot  of  the  Retran 
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missile.  Initially,  the  autopilot  g-blas  was  set  at  lg  and  guidance 
commands  were  not  transformed  to  the  missile  body  frame;  that  Is,  It 
was  assumed  that  the  missile  body  frame  was  sufficiently  aligned  with 
the  gimbal  frame.  This  Is  the  configuration  In  which  the  Retran 
missile  autopilot  would  be  employed  with  conventional  sensors  and 
guidance  (laser,  electro-optical,  etc.).  However,  In  the  turning 
trajectories  required  by  SAR-Retran,  missile  bank  angles  greater  than 
30  degrees  were  frequently  encountered.  These  large  bank  angles 
caused  sizeable  errors  In  the  orientation  of  the  resulting  missile 
acceleration.  Thus,  missile  response  to  guidance  commands  was  erratic 
and  sometimes  uncontrollable. 

To  eliminate  this  erratic  behavior,  the  autopilot  g-bias  was  set 

*■ 

to  zero  and  a g-bias  was  included  in  the  vertical  guidance  algorithm. 
Guidance  conmands  calculated  by  the  algorithms  in  the  gimbal  frame 
were  transformed  to  the  missile  body  frame,  so  that  no  matter  what 
missile  bank  angle  was,  the  resulting  accelerations  were  properly 
oriented  in  the  gimbal  frame.  With  these  modifications,  no  further 

i ' 

problems  were  encountered  with  the  missile  autopilot. 

Trajectory  Optimization 

The  attempt  at  trajectory  optimization  using  the  first  order 
gradient  algorithm  was  largely  inconclusive.  The  problems  which  were 
encountered  are  enumerated  here. 

Obtaining  satisfactory  initial  trajectories  was  difficult  with 
the  simplified  missile  model.  Since  instantaneous  response  was  assumed, 
missile  performance  and,  hence,  accelerations  commanded  by  the  guidance 
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algorithm  fluctuated  extensively.  Trajectories  achieved  in  the  simula- 
tion could  not  be  duplicated  in  the  Initial  solution,  but  reasonably 
shaped  trajectories  were  obtained. 

A sensitivity  analysis  was  carried  out  to  determine  what  values 
should  be  assigned  for  the  penalty  function  coefficients.  It  was 

“ *)h 

found  that  a value  of  K-»  > 1.0  x 10  6 caused  the  values  of  r — (t)  and 

J <>ah 

(t)  to  become  very  large,  while  values  smaller  than  1.0  x 10  6 had 

O Sy 

little  effect  on  moving  the  solution  toward  longer  range  intercept 
points.  Thus,  K3  » 1.0  x 10~6  was  chosen.  Similarly,  it  was  determined 
that  values  of  and  K2  larger  than  1000.  had  little  increasing 
effect  on  the  heading  error;  in  fact,  an  Improvement  of  10Z  in  heading 
error  was  the  maximum  which  could  be  achieved  and  this  was  obtained 
using  and  equal  to  1000. 

In  the  first  order  gradient  algorithm,  the  solution  converges 
to  a local  minimum.  In  the  Retran  problem,  local  minima  of  the 
payoff  function  were  encountered  with  only  minor  trajectory  changes. 
Widely  differing  initial  solutions  converged  to  different  optimal 
solutions,  with  similarly  small  trajectory  modifications.  Because  of 
the  occurrence  of  numerous  local  minima  in  the  function,  it  could  not 
be  concluded  that  any  of  the  resulting  trajectories  were  optimum. 

Indeed,  the  optimal  trajectory  could  only  have  been  obtained  by  inputting 
an  initial  solution  already  very  close  to  the  optimum  solution. 

Since  the  guidance  algorithms  and  coefficient  values  previously 
•elected  yielded  highly  satisfactory  results,  these  algorithms  were 
used  unchanged  for  the  launch  envelope  determination. 
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Launch  Envelope  Determination 

The  extent  of  the  usable  launch  envelope  was  determined  by  finding 
the  maximum  and  minimum  range  launch  points  from  which  the  missile 
could  satisfy  the  constraints  of  minimum  terminal  airspeed  and 
minimum  line  of  sight  tracking  time.  Maximum  launch  range  was 
limited  by  the  airspeed  constraint;  a terminal  airspeed  between  700 
ft/sec  and  750  ft/sec  was  chosen  as  criteria  for  a maximum  range 
launch.  The  tracking  time  constraint  was  always  satisfied  for  maximum 
range  launches.  Minimum  launch  range  was  limited  by  the  tracking 
time  constraint.  Terminal  airspeed  was  sufficient  In  all  minimum 
range  launches. 

It  was  felt  that  launch  points  can  safely  be  Interpolated  between 
the  test  points  used  to  determine  the  launch  envelope  boundaries, 
although  no  thorough  investigations  of  launches  from  points  within 
the  envelope  was  undertaken.  This  hypothesis  was  never  contradicted 
by  any  of  the  intermediate  range  launches  that  were  attempted  In 
the  process  of  determining  the  launch  envelope  boundaries. 

One  common  feature  among  the  20,000  ft  and  35,000  ft  launch 
envelopes  is  that  missile  flight  path  angles  were  such  that  the 
missile  maintained  an  almost  constant  airspeed  of  750  ft/sec  - 
775  ft/sec  during  the  latter  tracking  phase  of  the  trajectory. 
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VII.  Conclusions  and  Recommendations 


f U 

Conclusions 

Based  on  Che  performance  of  Che  hypoCheCical  ReCran  missile, 

1C  appears  chac  Che  SAR-Retran  guidance  syscem  could  be  demonscraced 
using  an  existing  boosc-gllde  air-Co-ground  missile.  A large  launch 
envelope  is  achievable  aC  high  alClCude,  Chus  providing  good  sCand- 
off  capabiliCy.  However,  Che  launch  envelope  becomes  very  narrow 
aC  a launch  alciCude  of  5000  fc;  Chus,  Chis  appears  Co  be  Che  minimum 
alClCude  for  successful  operacion. 

The  missile  aucopiloc  can  be  made  compacible  wich  SAR-ReCran 
guidance  by  eliminacing  Che  g-bias  from  Che  aucopiloc  and,  inscead. 

Including  Chis  bias  in  Che  guidance  algorichm.  The  guidance  commands 
muse  be  resolved  inco  Che  missile  body  frame  before  being  inpuC  Co 
Che  aucopiloc. 

The  firsC  order  gradienc  technique,  as  applied  Co  chis  problem, 
produced  inconclusive  resulcs.  Although  ic  indicaced  chac  only 
minuCe  CrajecCory  changes  were  required  Co  achieve  a local  minimum 

9 

in  Che  payoff  funcCion,  Che  soluCion  Co  which  ic  converged  was  noC 

necessarily  an  absoluCe  minimum.  I 

Recommendations 

A 

The  effect  of  aircraft  post-launch  tactical  maneuvering  on 
missile  performance  and  launch  envelope  should  be  studied.  A 

I maneuvering  aircraft  could  be  included  in  Che  Retran  simulation. 

The  effect  of  maneuvering  targets,  such  as  a tank  or  truck, 
should  be  investigated. 

( 
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A higher  order  optimization  technique  might  be  employed  to 
optimize  missile  trajectories  for  maximum  range* 

In  view  of  the  range  limitation  of  present  air-to-ground  missile, 
as  represented  by  the  Retran  missile,  particularly  at  low  altitude, 
performance  specifications  should  be  developed  for  new  tactical 
missiles  to  be  employed  with  SAR-Retran. 
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Appendix  A 

Coordinate  Transformations 

The  following  coordinate  transformations  are  utilized  in  both 
the  Retran  simulation  and  the  Optraj  optimization  program.  The  order 
of  rotation  from  the  earth  fixed  frame  to  the  missile  and  aircraft 
body  frames  is  (Ref  1:  2.16) 

1.  Yaw  (ip  >4>A) 

2.  Pitch  (0,0A) 

3.  Roll  (<J>,<j>A) 

The  order  of  rotation  from  the  aircraft  body  frame  to  the  glmbal  frame 
is 

1.  Yaw  (i|>GA) 

2.  Pitch  (0GA) 
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where 
[LiaI  - 
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[L^]  is  identical  to  [LM]  except  for  deletion  of  the  subscripts 
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Missile  Body  Frame  to  Earth  Frame 


Y 

% 

Y 

- t^r1  [er1  ttp 

ym 

Z 

— m 

ZM 

(93) 


■ (limJ 


r^M 

yM 

ZM 


(94) 


Ilim1  is  identical  to  [Lj^]  except  for  deletion  of  the  subscripts. 
Aircraft  Body  Frame  to  Gimbal  Frame 

(95) 


V 

V 

yg 

ya 

_ZG 

_ZA_ 

" tLGAJ 


where 


[LgaI  " 


C0GAC<^GA 

*S^GA 

S9GA%A 


C9GAS*GA  -S0 


GA 


Op, 


GA 


s9gas^ga  c9ga 


(96) 


(97) 
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«*ROGRiH  1- T?ftM(irpjr,0JTPUT,rAPE5=IN9JT,rftP;9  = 0UTPJT,>L0T,PjNCHI 
C 



C* 

C* 

C • RE  TRAN  HISSIE  SIMD.ATION  . 

C* 

C* 

c 

c . . 

c» 

C*  ' DEFINITION  OF  INPUT  VARIABLES 

C* 

C*  XEA,  YEA,  ZEA  - INITIA.  A I P.CRAFT  POSITION  IN  EARTH-FIXED 
C*  COORDINATES  « - T> . 

C*  . 

C*  XEAO,  YE  AO , Z E AO  - INITIAL  AIRCRAFT  VELOCITY  IN  EARTH-FIXED 
C*  COORDINATES  (rT/SEC). 

C* 

C*  THAOG,  3SA0G»  =HADG  - INTTIAL  AIRCRAFT  PITCH,  HEATING,  ANO  ROLL 
C*  ANGLES  IN  DEGREES  WHICH  ARE  CONVERTED  TO  RADTANS  (TH,PS,»H) 

c* 

C»  XET,  YET,  ZET  - TARGET  POSITION  IN  EARTH-FIXEO  COOROINATES  (FT) 

C* 

C*  H - INTEGRATION  STEP  SIZE  IN  SECONDS. 

C» 

C*  G8IAS  - BIAS  CDMH AND  IN  VERTICAL  CHANNEL  IN  G*5 
C* 

C*  GHHIAS  - BIAS  COMMAND  TN  HORIZONTAL  CHANNEL  IN  G'S 
C» 

C*  TINE  - NAXINJN  FLIGHT  TINE  (SEC).  ( PROS  RAH  TERMINATES  HHEN  MtSSLE 
C*  GOES  THROUGH  TARGET  ALTITUDE  OR  HHEN  THE  IS  EXCEEDED. ) 

C* 

C*  ITABL  - INTEGER  SWITCH  CONTROLLING  PRINTING  OF  AERODYNAHIC  TABLES. 
C*  PRINTING  TAKES  PLACE  ONLY  IF  ITABL  = 1 . 

C* 

C*  IPR  - THE  NJH3ER  OF  SECONDS  BETWEEN  PRINT  CYCLES  (INTEGER) 

C* 

C»  PBIAS  - 3IA5  ON  AUTOPI.OT  ROLL  RATE  GYRO  IN  RADIANS  t SECONO. 

C* 

C*  TR  - ACCELERATION  COMMAND  COEFFICIENT  REFERENCE  TIME  (SEC). 

C* 

C*  ZI  - OAM^TNG  RATTO  IN  COHMANO  TO  LOS  GUIDANCE  ALGGRIf HNS 
C*  • 

C*  CT  - GAIN  COEFFICIENT  ?0R  COMMAND  TO  HEADING  GJIOANCE  A.GORITHM 
C* 

C*  F - PLOTTING  SIZE  FACTOR  F = 1.0  YIELDS  b.IN  3Y  BIN  PLOTS. 

C* 

C*  IPU  - INTEGER  SWITCH  CONTROLLING  PJNCH  OUTPUT,  PUNCHING  TAKES 

c*  place  only  if  i»j  = 1, 

c* 

C*  IPT  - INTEGER  SWITCH  CONTOLLING  PLOTTING.  PLOTTING  TAKES  PLACE 
C*  ON. Y I?  IPT  » l, 

C* 

C*  ILEG  -INTEGER  SNITCH  CONTROLLING  PRINTING  OF  LEGINO  OF  OUTPUT 


YARIA3LiS.  LEGEND  IS  POINTED  ON.Y  IF  ILEG  =’  1. 


u 


C* 

C*  • 

•'***•*» 

c 

COMMON 
COHO'! 

COMMON 

1 

COMMON 

1 

COMHO'! 

COM MOM 
COHMO'l 
COMMOH 
COMMOM 
COMMOH 
COMMON 
COMMON 

1 
1 
1 

COMMON 

1 

COMMON 

1 

COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
C 

c 

C* 

C*  FLIGHT  INITIALIZATION  * 

C*  * 

C 

10  CONTINJE 
C 

INITIALIZE  INTEGRATION  STEP  COJNTER 
C 

IJK  » -1 
C 

C»*****»  INITIALIZE  PLOT  COJNTER 
C 

NPTS=0 
NPT2  * 0 
MPT  3 =0 
C 

C****»*»  REAO  INITIAL  FLIGHT  CONDITIONS 
C 

REA9(5 , 3013 1 XEA,  YEA*  ZEA,  XEAO,  YSAO , ZEAO,  TMA03»  PSAOG,  PHADG 
1 * X£T , YET,  ZET,  H,  G9IAS , GHOIAS,  TIME,  IIA3L,  I PR,  PBIAS.TR, 

1?I,C7,c»I3J,IPT,ILEj 
C 

9010  FORMAT!  3-10  2,  / 5F3.?,  2F&.2,  F3.2,  12,  14,  e5.  3 , *4.1  /JF10.0, 


/ACS/  axb,avq,azo,ax£, aye, aze.poot, oodt, root 
/AEJ.ER/  PHAODT, TMAOOT , PSAOOT 

/AERY  OEL,OOC,3N1 , CM  1 , CC ,CM2 , 0N2 , CM, 300 , OEL SO,  AM ACH, A LPHA 
,alt,9Eta,::y,:m2y,on2y,cmy,cniy,:miy,occy,:o,:n,cny 
yam;/  PS,  TH,  ’■«,  PSOG,  TM  OS , »M3S  ,P3A,  T M A , aH  A, 

PSADG,  TM1DG,  PHAOG,  =»S3ADG,  THGA05,  HGAOG 
/ 3.0  SKI/  XXI, YYt ,7ZI, AMASS, CG, HG,PBI AS, << 

/•T.D;K2/  AHS,  4Y3,TTG,DELF,ISUID,TR,UA,  RANSE.ZliCT.YOOT 
/CMC/  47,  AM,  MY,  rf? 

/CMr.i/  no, oo, ot, acz,acy,ac<,s3Ias,ghbi as,i jk 

/CNT.2/  R,a,3, J,0,31,P1,P2,<,» 

OER/  O,RL,CMA0,CMCG,QQ,S, YEL»THR,CMADY,CM03Y,RY 
/INI TL/  H, TIME, ITAOL,IOR, ILEG 

/MATRIX/  CETl,0;T?,CET3,CEr‘+,0ET5,0ETs,0ErZ,CEr8,CET3, 

:iSl,0IS2,:iS3,:iG4,CIG5,0ISo,:iGJ',CIS3,CI53, 

CM31,CMG2,CMG3,:hG4,CMG5,0MG6,0MGT,CMG3,CMG3, 

SPSS  A,  O^SGAtSH  5 4, CTHGA 

/°L04TA/  PXET<  5 2 ) , PYE T ( 32J,9XE(  52),3YE(  52),®XE4(  52), 

DYE4 ( 52), PZE<  r 2 ) » PZEA ( 52),- 

/P.3AT2/  PAH3 ( r2),®AYG<  52),>T(  52),=3e(  5 »), 

“R(  52)  , PAL  ( 52),PBE(  52),P70K  5?),«>0Y(  32),»0Z<  52) 
/°L3AT3/  PAYS  C 72>,P3ES<  32>,STS<  32),  PR3(  32)  , IPRO 
/POS/  XET, YET, ZET, XE A, YE A, ZEA, R AN AT, XE, YE, ZE 
/®RNTl/  DY,OZ,07Y,OVZ, AYP, AZP 
/7TARGET/  XET 3 , YETO, ZET0.X3T 0, YBT3,Z3T3 
/THRST/  TIME3(2l) ,THSTLO  (21) 

/YE. 0/  XEO, YEO, ? EO, XEA O.YEAO, ZEAO, JB,Y3,N9 


j 

I 
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I 


j 


jt 


i 


i:2,*X,T’fMfI2> 

c 

c** ** * •**  t :gt  end  oc  jdb 

c 

IF(E0F(5n  75*15 

c 

c** *****  S£T  oommanoeo  azimjth  gimoal  angle  of  aircraft  ipsgc> 

c 

15  OSGC  = -55. 

C 

o:rIS"  AIRCRAFT  VELDCITY  (UA> 

c 

IJA  = X:  AO 
C 

C*-****‘*  PRINT  INITIAL  CONDITIONS  AND  AERODYNAMIC  TA3LE3 
C 

CALL  ODUMT. 

IF(I°U» "0. 1)  PUNCH  12, XE  A > YEA  , ZCA 
12  cOPHAT ( 3P13 . 2> 

C 

c*******  initialize  flight  parameters 

c 

CALL  INTLF.T 
IPRO  = I PR 

IPR  = T3R/M  * 

HO  = H. 

C 

Cm««**m*******»»*»*m»m*»***m*mm*m*m»m«om*m*m*»m*ommi> 

C* 

c*  THIS  SE3UEN3E  15  FOLLOWED  FDR  EACH  INTEGRATION  STEP 

C* 



C 

20  CONTINUE 
C 

C*»«***  INCREMENT  INTEGRATIDN  STEP  COUNTER 

c ; 

IJK  = IJK  F 1 
C 

C*******  03T1IM  TARGET  VELDCITY 
C 

>'  CALL  TARGET  (X) 

C. 

C**“*  + M+  COHPJTE  GDOROlNATi  TRANSFORM AT I DNS 
C 

CALL  EULER! 1) 

C • 

C*;****»*  COM0 JTE  GUIDANCE  3 OH  HANDS  ANO  CONTROL  DIS3L  ACEHENTS 
C 

CALL  COMMAND 
. CALL  CDNTRDL 
C 

C*»‘**‘*  COMPUTE  AIR  DENSITY  ANO  SPEED  Dr  SOUND 
C 

ALT  * -7E 

CALL  ATMOS  (ALT ,RHO, SD5) 

C 

£**»»» **  GOMPJTE  ANGLE  DF  ATTACK, SIOS-SLIP  air  SPEED,  MACH  HUMMER,  AND 


( 
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c oynh!:  pressure 

c 

ALPHA  s ArAN(H3/U9)  *57  3 

beta  = &T>^(va/un)  •;?,  3 

VEL  = SORU  XEO’XED  f YiD’YEO  ♦ ZcO'ZEO) 

AMACH  = VE;./ SOS 
QO  = 5*VE.*VEL*RHO 
C 

C»*«**»*  COMMUTE  THRUST,  MOUNTS  OF  INERTIA,  AND  CENTER  OF  GRAVITY 
C 

IF(X- LT.  4.  5)CALL  THRJST (X) 

C 

£«,«»•«*  pg^  AE RD3Y NAMIS  FORCES  AND  MOMENTS  ON  MISSILE  AND 
C COMMUTE  THE  ACCEL ERDMETER  OUTPUTS  IN  BODY  r RAME  AXES. 

C 

CALL  AERO 
C 

c*«.*»«*  COMMUTE  MISSILE  AIRSPEED  ACCELERATION/DECElERATION 
C 

VOOT  = (AX3) *COS( ALPHA/57.  235) *C0S(3ETA/>7.  295) *AY3*SIN( 

1 BETA/57  235)  fAZB*SIN(  ALPHA/57.  29  5)  *00S  ( 3£T  A/5  7.  23  5> -32. 2*  SI  N ( TH- 
iALPHA/57. 235)  • 

C 

C*»»****  TRANSFORM  ACCELERATION  COMPONENTS  OF  MISSILE  FROM  300Y  AXES  TO 
C EARTH  SYSTEM.  <A<3, AYB.AZB)  TO  ( AXE, AYE, AZE) 

C 

CALL  E JIER( 2) 

C 

c. ««»,,,  If|T  E SRA  T E ACCELERATION  COMPONENTS  OF  MISSILE  ( AXE, A YE, AZE) 

C - GIVTNG  VELOCITIES  (XEO ,YEO, ZEO) 

C 

YED  = YEO  v AYE*H 
ZEO  = ZED  ♦ AZE*H 
XEO  = XEO  ► AXE*H 
C 

C*******  INTEGRATE  ANGULAR  ACCELERATIONS  0-  MISSILE  (500 T , ROOT) 

C GIVING  300Y  ANGULAR  RATES  (Q,  R) 

C 

0 x Q i-  QOOT*H 
R « R ► RDOT*H 
C 

COM»JTE  EULER  AN3.E  RATES  (PSOOT,  HOOT, PHOOT) 

C 

PSOOT  = (3* SIN (PH)  ♦ R*COS(PH))  / COS(TH) 

THDOT  = R*:03(°H)  - R*SIN(PH) 

PHDOT  « PSOOT*SIN(TH)  ^ PBIAS 

C * * 

C*****»»  INTESRATE  EULER  ANG.S  RATES  (PSOOT,  THOOT,  ®HDO  T) 

C GIVtNG  MISSILE  ATTIT JOES  • (PS,  TH,  P«) 

C 

PS  « PS  ♦ >SOOT*H 
TH  * H ♦ f HQOT*H 
®H  » PH  ♦ *HOOT *H 


C 

C**‘**‘*  COHMANO  AIRCRAFT  TJRNING  RATE 
C 

C PS5ERR  IS  THE  ERROR  BETWEEN  THE  COMMANDED  AIRCRAFT  SIM9A.  ANGLE ' 

C AMP  THE  ACTUAL  SIM3AL  ANGLE.  THIS  ERROR  IS  LIMITEO  TO  PLUS  OR 


” 


c 


INTEGRATE  AIRCRAFT  NELOCITIES(XEAO,YEAO>  GIVING 
C AIRCRAFT  POSITION  ( X£A»YEA) 

C 

30  XEAQ  = JA*003<THA> *C3S(°SA) 

YEAD  = JA*C03  (THA)*SIN<PSA> 

ZEAO  = JA*3IN( THAI 
XEA  = XEA  f XEAO*M 
YEA  = YEA  v YEAOM 
ZEA  = ?E  A ► ZEAO*M 
C 

INTEGRATE  HISSILr  VELOCITIES  < XEO, YEO, ZED) 

C GIVING  MISSILE  POSITION  <XE,YE,Z£) 

C 

XE  = X:  ♦ XE9*H 
YE  = YE  ♦ YE3»M 
TE. = ZE  ♦ ZE3*H 
C 

C»*4***»  INTEGRATE  TARGET  /E-OCITIES  <XSTO, YETO,ZETO> 

C GIVTNG  TARGET  POSITION  (XET »YET , ZET) 

C 

XET  = XET  f XETO*H 
YET  = YET  ¥ YETD*H 
ZET  = ZET  f ZcTD*M 
C 

C**‘**«  TEST  for  impact 

c 

IF(ZE,3T.Z£T)  GO  TO  VO 
C 

<;♦*.**♦*♦  INCREME XT  TIME 
C 

IFTTTG. UT. 10. ) H=M0/10 
IFtTTG.'.T.lO  ) IPR  = 1 /H 
X * X ♦ H 
C 

C**‘**‘*  TEST  INTEGRATION  STEP  COUNTER  FOR  PRINTING 


! 


«»vrt,i«»rs>=f£ 

»7EtUPT«?>  * -ZE 
»T(NPT5>  * * 

PP.CIPTD  * rance 
"AHGtNPTS)  s amg 

•TF<3Ah;<N3TS).GE.  20.)  3 AUG  (NPTS)  = 20. 

ir<Pt-r,(N3r3> . LE.-20, ) 3 AHG  < M P T S>  = -20. 
PAVG  ( M3T  S)  = AVG 

IFCPAVSOPTS) « GE.  20.)  PAVG(NPTS)  = 20. 

IFIPAVK’PTS)  .LE.-20.I  PAVGtNPTS)  = *20. 
PAL(NPrS)  * ALPHA 
»Df(NPTS>  s BETA 
PDF(NPTS>  = OELF 

IF(3r>P(N°TS  >.GT.  15.)  3DF( NPTS)  = 15. 
PVOTUIPTS)  = VQOT 

TFCPvnrrprs) . ce.  io),»  pvoknpts)  * 100. 

IFtPVDT (M3rS> . LE. -10). ) PV0TCNPT5)  * -100. 
IF (X- or. 12. ) GO  TO  41 
HPT?  = NPT2H 
PAVS(N3T?)  = AVG 

IF( PA VS (N3T  2) .LE.-20.)  3AVS(NPT2>  * -20. 
TF(PAV?(N3r?) . GE.  20.)  PAVS(NPT2>  a 20. 
P3FS(M3r?>  = BETA 
f*TS(M°T2)  s X 
60  TO  43 


41  HPT 3 = .NPT3H 
®0Y(HPT3)  s DY 
PD? ( HPT 3)  * OZ 
PR3  (NPT.3)  a RANGE 

43  IF( IDU.  ER.  1 ) PUNCH  42,<E,YE,ZE,X£A,YEA,ZEA 

42  FORMAT (3-10.2/3F10.2* 

C 

PRIST  PLIGHT  VARIABLES 
C 

' CALL  OJTFlT 
C • 


50  CONTINUE 
C 

TrSr  r3?  CONTINUATION  OF  FLIGHT 
C 

IFCX  LT.TISE  ANO.ZE..T  ZED  GO  TO  20 

C >• 

C**-*»‘»  ENO  Oe  FLIGHT 
C 

C*******  PLOT  TRAJECTORY 
C 

IF(IDT. EQ. 1)  CALL  3I<(SPTS,NPT2,NPr3) 
C . . 

GO  TO  10 
C 

75  CONTINUE 
C 

c»*«»*«*  ENO  OF  J03 
C 

STOP  1 

c 

ENO 
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I 


subroutine  aero 


c 

r.*< 

c* 

c* 

c* 

c* 

c* 

c*> 

c 


thts  routine  computes  forces  and  homents  given  angles  d-  attack, 

CONTROL  DEFLECTIONS,  AMO  MACH  NUMBER.  DUT°'JT  IS  TRANSLATIONAL  AND 
ROTATIONAL  ACCELERATION  COMPONENTS  IN  BODY  AX-J. 


C 

c*> 

c* 

c* 

c* 

C*‘ 

c 


COMMON  /ADC/  AXP, AY3.A73, AXE, AYE, A2S,POOT,aOOT, ROOT 
COMMON  /AS*/  DEL » OGC » CNl'.CM  1,  CC »CM2 » CM2 » CM,  COO » DEL  CO.  AMACH,  ALPHA 
1 ,ALT,QETA,CCY,CM2Y,CM2Y,CMY,CM1Y,CM1Y,  3CCY,  CO,  CN.CMY 

COMMON  /AERT3L/  0DCT<9,5,6),  CMT( 8 , 5 » 5) , CMT (3,5,5)  , <1(19>, 

1 CCTliO.S) ,5M2T( 10,5) ,CN2T (10,6) ,CMTT(iO,S),<2<15», 

2 CDDT ( 10) , MCMCOT (10)  ,H(5)  , Oi.COT (61 
COMMON  /AMS/  RH02(12> ,?0S2( 12) , ALTT (12) 

COMMON  /OlDCKI/  XXI, YYI, 221, AMASS, CG,HG,PBIAS, XX 
COMMON  /CMT.!/  DB,05,DT,ACZ,ACY,A5<,S3IA3,SM3IAS,IJK 
COMMON  /CMT  L2/  R,  0,  G , J,  0 , B1 , »1,  P2 , X , » 

COMMON  /OER/  D.RL.CMAD.CMCG.QQ.S, VEL.THR.CMADf.CMCSY.RY 
COMMON  /TMRST/  TImE5(21) ,TMSTL0(21> 

COMMON  /VE.O/  XED»YE)»?£0,X£A0,YEA),ZEA0,U8,V3»W3 

DIMENSION  MA(2) ,XA(2t , V A 1 ( 3) , XA1( 31 


4ER00YNAMIC  TAB.ES 


C 

c 


DATA  {((0051(1, J,<),t *1,8), J= 1,5)  ,<=t,  3) 


i 

0.00 

9 

.010 

9 

.0  45 

, .110 

t 

. 230 

9 

.333 

9 

.300 

9 

.300 

9 

2 

-.01 

9 

0.  000 

9 

.033 

, . 090 

9 

. ISO 

9 

.253 

9 

,253 

9 

.250 

9 

3 

-.02 

9 

-.010 

9 

.020 

, .075 

9 

. 140 

9 

.203 

9 

,203 

9 

.200 

9 

4 

03 

9 

-.025 

9 

0.  30 

t .040 

9 

.033 

• 

.130 

, 130 

9 

.130 

9 

5 

-.04 

-.040 

9 

-.32 

, .005 

.030 

9 

.050 

9 

.050 

9 

.050 

9 

6 

0.00 

9 

.010 

f 

.04? 

, .110 

.200 

9 

.390 

9 

.300 

9 

.300 

9 

7 

-.01 

9 

0.000 

9 

.033 

, .090 

9 

. ISO 

9 

.253 

9 

.253 

9 

.250 

9 

8 

-.02 

9 

-.010 

9 

.020 

, .075 

9 

. 140 

9 

.203 

9 

.200 

9 

.200 

9 

9 

-.03 

« 

-.025 

t 

0.  30 

, .040 

9 

.030 

9 

.139 

9 

.130 

9 

.130 

9 

A 

-.04 

9 

-.040 

9 

-.32 

• .005 

9 

.030 

9 

.050 

9 

.053 

9 

.350 

9 

B 

0.00 

9 

. 010 

« 

.0,5 

* .110 

9 

. 200 

9 

.300 

9 

,300 

t 

.300 

9 

C 

-.01 

9 

0.000 

t 

.033 

, .090 

9 

.ISO 

9 

.2,0 

9 

.253 

9 

.250 

9 

0 

-.02 

t 

-.010 

9 

.020 

, .075 

9 

.140 

9 

.209 

9 

,200 

9 

.200 

9 

c 

-.03 

9 

-.025 

» 

0.00 

* .040 

9 

.030 

9 

.130 

9 

. 130 

9 

.130 

9 

F 

-.04 

9 

-.040 

9 

-.32 

, .005 

9 

.030 

9 

.053 

9 

, 05) 

9 

.050 

/ 

DATA  (((0 

CCTd.J.O ,1*1, 

6)  , 0=1,5)  ,<=♦,&> 

- 

/ 

1 

0.00 

9 

.010 

• 

.0,5 

, .110 

9 

.200 

9 

.300 

9 

.300 

9 

.300 

9 

2 

-.01 

t 

0.000 

V 

.033 

, .090 

9 

. ISO 

9 

.250 

9 

, 250 

9 

.250 

t 

3 

-.02 

9 

-.010 

9 

.020 

, .075 

9 

. 140 

9 

.200 

9 

,20) 

9 

.200 

• 

4 

.-,03 

9 

025 

9 

0.30 

, .040 

9 

.030 

9 

.130 

9 

,130 

9 

.130 

t 
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5 

-.04 

9 

-.040 

9 

-.3  2 

9 

.005 

9 

.030 

9 

.030 

9 

.053 

9 

.050 

. 9 

6 

0.00 

9 

. 010 

9 

« 0 45 

9 

. 110 

9 

. 200 

9 

.300 

9 

.300 

9 

. 300 

9 

y 

•.01 

• 

o.noo 

9 

.033 

9 

.090 

9 

.1>0 

f 

.230 

9 

.253 

9 

.250 

9 

8 

-.02 

• 

-.010 

9 

.020 

9 

. 07 5 

9 

.1.0 

9 

.200 

9 

.200 

9 

.200 

9 

<» 

-.03 

• 

-.025 

9 

0.00 

9 

. 040 

9 

. 030 

9 

.133 

9 

.131 

9 

.130 
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OAT  A (C(Cir  (I.J.O ,I=t ,6>,J=1,5>,<=1,3>  / 
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F 
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9 
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9 
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DATA 

« C C CMT  (I*. 

J,<>  .1  = 1,8) 

t Js If  5) 

/ 

• 

c 

; 

1 

.0 

» 1.0 

* 2.3 

9 

3.0 

t 4.9 

9 

5.0 

9 

5.0 

9 

5.0 

• 

2 

.0 

, 1.0 

• 2.3 

9 

3.0 

* 4.3 

9 
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9 

5.0 

9 

5.0 

3 

.0 

» 1.0 

» 2.0 

9 

3.0 

* 4.3 

9 

5.3 

9 

5.0 

9 

5.0 

4 

.0 

» 1.0 

t 2.0 

9 

3.0 

t 4.0 

9 

5.0 

9 

5.0 

9 

5.0 

5 

.0 

* 1.0 

t 2*1 

9 

3.0 

» 4.3 

9 

5.0 

9 

5.0 

9 

5.0 

w 
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.0 

, 1.0 

» 2.3 

9 

3.0 

* 4.3 

9 

5.0 

9 

5.0 

9 

5.0 

1 

7 

.0 

» 1.0 

> 2.0 

9 

3.0 

t 4.0 

9 

5.0 

9 

5.0 

9 

5.0 

8 

.0 

» 1.0 

. 2.1 

9 

3.0 

> 4.3 

9 

5.0 

9 

5.0 

9 

5,0 

9 

.0 

» 1.0 

. 2.9 

9 

3.0 

» 4.0 

9 

5.0 

9 

5.0 

9 

5.0 

. 

A 

.0 

* 1.0 

• 2.3 

9 

3.0 

» 4.0 
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9 
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t 4.0 
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.0 

1 1.0 
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9 

3.0 
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9 
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9 

5.0 

• 

5.0 

/ 
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OAT  A (( C O'*  T <I,J,<>,t*l,R>,J=l,5>  ,<  = 1,3>  / 
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OATA  C < CC'IT  ( I , J , K)  ,1  = 1,8),  J=  1,5)  » < = V, 6)  t 
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OATA  AITT 
C 

1 0 

2 19890. 

C 

C*»***‘» 

C 

DATA  S3S2 
C 

1 Ills. A 

2 1037.7 
C 

C*****»» 

c 

OATA  RH02 
C 

1 .23770 

2 .12762 


C 

C 

OATA  COOT 
C 


/ 
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, 24000. 
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. 1110.  5 

, 1020.3 
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, .22753 
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C’ 

c 
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oata  h:hcot 
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/ 


• 3000. 

, 28000. 
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C' 

c 

c 
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C' 
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c 

c 

c** 

c 

c 

c 
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c* 

c* 
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1 . 4 , • 6 , 75  , . 35  » .9-  , 


OATA  H / 
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OATA  O'LCDT  t 

100  « *.0008  , -.0014 


OATA  N»  / 10,  6 / 


OATA  Nil  / 8,  5,  5 / 


SOIPUTATIOMS  BEGIN 


1.0  , 1,1 


, 12000. 


, -.  0024 
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, 1.3,  1.5,  2.0 


, 20000.  , 30000. 


, -.0052  , -.0100 
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c 

Finn  AXIAL  FORCE  EFFICIENTS 
C 

CALL  imLJI’.  6,H,0E.CDT  , AL  T » OELCD,  IE » 0) 

CALL  NDTLJ(2,10,MCHCDT,CDOT,AMACH,  COO, IE, 1) 

c 

£«*«*»««  FIND  iONENT  COEFFICIENTS 

c 

XA  ( 1)  = SETA 
XAC2)  * ANACH 

CALL  NDTL'J{3»NA,X2» 3 M2T  , XA,CM2Y»IE»  t) 

CALL  NDTLJ( 3,NA  ,X?,CN2f ,XA,CN2Y,IE,  1) 

CALL  N3TLJ<  3,NA,X2,Cm,  XA,CMY,IE,l> 

C 

C*******  FIND  “ITCH  MOMENT  CD  EFFICIENTS 
C 

XA(  1)  = AL=MA 

CALL  NDTLJ(  3,.NA,X2,Ci2r,XA,CM2,IE,l> 

CALL  NDTLJ(3,NA,X2,CN2r, XA,CN2,IE,l) 

CALL  NOTLU(3,NA,X2,CUT,XA,CM  ,IE,1) 

C 

C*******  FIND  IND’JCEO  DRAG  COEFFICIENT 
C 

ALPHAT  = AT  AN ( SORT ( W?*WB  ♦ VB*VB  > ( U3  >*57.3 
XA ( 1)  * AL3HAT 

CALL  NDrLJ(3,NA,X2,CDT  ,XA,CC  ,IE,1) 

C 

C** *****  FIND  Y AX  -ORCE  AND  MOMENT  COEFFICIENTS 
C 

XAiCi)  = A33 (DC) 

XA1 (2)  = A3S<  BETA  > 

XA1 (3)  = AMACH 

CALL  N^TLJf  4,NAl,Xl,CMr ,XA1,CM1Y,IE,1) 

CALL  NDTLJt  !*,NAl,Xl,DNTtXAl,CNlY,IE,l» 

C 

C*******  FIND  PITCH  FORCE  AND  MOMENT  COEFFICIENTS 
C 

XA1 (1)  = A3S(OB) 

XA1 (2)  * A 3 S ( ALPHA! 

CALL  NOrL‘J(V,NAl,Xl»CMr,XAl,CMl,It,ll 
CALL  Un.TLJ('4tNAl,Xl,CNr,XAl,CNl,IE,l» 

C 

C*»»***»  FIND  INDJCED  ORAG  CDEFFICIENT 
C 

XA1 (1)  = A3SCOB) 

XA1(2)  = A3S(ALPHA> 

CALL  NDTLJI  ^»,NA1,X1,  DCCT  ,XAi,  OCC,  I E,  1) 

XAi(l)  = A3S(OC) 

XA1 (2)  = A 3 S C BETA) 

CALL  NDrLJU,NAl,Xl,3CCT,XAl,OCl,IE,l> 

OCC  = DCC  -»  OC  1 
C 

C**»**,,  COMBINE  COEFFICIENT? 

C 

CO  * CDOOCC  4CC40SLC0 
SINOEL  = SI GN < 1 . ,03) 

SINOER  » SIGN ( 1. ,0C> 


99 


0 


ON  = CN2-0N1*  SINOEL 

omao  = om?-cmi.*sinoe. 

0MC0=C-.T  7-0*,>  *CN/12. 

CNY  = CN.RY-CNIY'SIND:  R 
CMAflY  = CM2Y-CM1Y*SIM0ER 
CMC3Y  S U3.  7-CG)*0Nr/l2. 


C 

C« 

C 

C 


»*•  LINEAR  ACCELERATION  COMPONENTS  OF  MISSILE  IN  BODY  AXES  OJE  TO 
AERODYNAMIC  FORCE  5 AND  THRUST 


AXR  = CTU-S*Q0*CD> /MASS 
AYB  = -3*00*DNY/AMA3S 

AZ3  = -s‘od*cn  /amass 

ACX  = A XI/ 32. 2 
4CY  = AY3/32- 2 
ACZ  = AZ9/32. 2 


C*»**»»»  ANGULAR  ACCELERATION  COMPONENTS  Or  MISSILE  IN  9D0Y  AXES  DUE  TD 


C 

c 


AERODYNAMIC  MOMENTS 


°OOT  =0. 

OOOT  = 5*QDM  <CM*Q/(2  *VEL>>  + (CMAD-CMCG)  > / YYt  ♦ XX*R 
ROOT  = 3»3DM  <CMY*R  / ( 2 . * VEL ) ) - (C1ADY-CMC3Y)  ) / ZZI  - XX*Q 


C 

C 


RETURN 


END 


C 

C*a 

C* 


SUBROUTINE  ATMOS  ( AL T , RHO, SOS ) 


C*  THIS  ROUTINE  COMPUTES  AIR  OENSITY  AMO  SPEED  OF  SOUND  AS 
C*  A FUNCTION  OF  ALTITJOE. 

C* 


c 

c 

c 


COMMON  /AMS/  RM02( 12> , S 0S2 ( 12 ) , A. TT ( 121 
DIMENSION  N A C 1 ) ,XA(1) 


NA  ( i ) = 1 2 
XA(1) =ALT 


C*“*»*»  AIR  OENSITY 
C 

CALL  N0TLJ(2,NA, ALTT , RM02 ,X A , RHO, IE. 0J 
RHO= . 01  * RHO 
C 

£•»*»***  SPEED  0?  SOUNO 
C 

CALL  NOTLJ(2»NA,ALTTfSOS2*XA»SOS»IE»D) 
C 

RETURN 

c 

END 

SUBROUTINE  COMMANO 
C 


- 


* " 

t a; 


res’*  I 

fft 


r 

i*n. 


L ^ 

fc—  v* . 

F r.'v 

Her m 


100 


c*  t hts  routine  computes  th:  position  ano  velocity  op  the  missle  with 

C*  PESPECT  n THE  AIRCRAFT-TO-TARGET  .INE-DP  -SIGHT  AMD  IS3JES 
C*  ACCELERATION  COMMANDS  (AV,  AH)  TO  3UIOE  THE  MISSLE  DOWN  THE 
C*  L I ME-OF-FIGHT  • 

c* 

c 

COMMON  /a:;/  AXQ.AYO, A73.AXE, AYE, AZE, aOOT,OODT, ROOT 
COMMON  /AEJ.ER/  PHA03T, THAOOT .PSAD3T 

COMMON  /AMS/  PS,  TH,  PH,. PSDG,  THDG,  PHOS  ,P3A,  THA,  PHA, 

1 PSAOG,  THA  OG,  PH A DG,  =S3AOG,  THGADG,  HGADG 

COMMOM  /3..D0K2/  AHG , A V 3 , TTG  , OELF,  I 3 UI  D , TR,  UA,  * AN3E  , ZI  , C T,  V ODT 
COMMON  /CMC/  AV, AH, HT, HZ 

COMMOM  / C N T _ 1 / 03, DC , DT , ACZ , AC Y ,ACX , 331 AS, GHBI AS , I JK 
COMMOM  /CMT-2/  R , 0, 3 , J , C , 81 , PI, P2 , < , ? 

COMMON  /OER/  0,RL,CMAD,CMCG,QQ, S,  VEL  , THR, CM  AD Y , CMC  3 Y, RY 
COMMON  /INITL/  H, TIME , I T A3L  , I PR ,1 . EG 

COMMON  /MATRIX/  CET l , C : T 2 , CET 3 ,CE T ♦ , CE T5, CE To , CET 7 ,CET 8 ,CE T 3, 

1 CI3l,CIG2,CIG3,CIG4,CIG5,CI3o,CIG7,CI3J,CI33» 

1 CM31,CMG2,CMS3, CMG4, CMG5 ,CM36, CMS 7 , CM3  3 , CM3  3 , 

i S°SSA,CPSGA,SrHr'A,CTHGA 

COMMON  /®D3/  XET,YET,Z‘T,XEA,YEA,ZEA,RANAT,XE,/£,ZE 
COMMON  /°RNT 1/  DY,DZ,OVY , OVZ » AYP, AZP 
COMMOM  /VE-D/  XED, Y£D,7ED,XEAD, YEA 3 ,ZEAQ,JB,/3»N3 
COMMON  /VTAR3ET/  XETD, i l TO, ZETD,X3T 0, Y BTO, ZBTO 


C*******  FI  NO  POSITION  DP  TARGET  RELATIVE  TO  AIRCRAFT  (EARTH  AXES) 

C 

XAT  = XET  - XEA 
YAT  = YET  - Y£A 
ZAT  = ZET  - ZEA 
C 

c, ♦***.,  find  VE. D3ITY  op  r A R3ET  RELATIVE  TO  AIRCRA  rT  IM  EARTH  FRAME 
C (EARTH-  AXES) 

C 

XTAOI  = XETD-XEAO 
YTAOI  = YETD-YEAO 
ZTAOI  = ZETD-ZE AO 
C 

£»«•«**»  fin:j  am3 JLAR  VELOCITY  OF  AIRCRAFT  FRAME  RELATIVE  TD  EARTH  FRAME 
C (HA)  (EARTH  AXES) 

C’ 

WAX  = PHADDT*COS(THA> * COS (PSA ) -THA DOT* SIN (PSA ) 

WAY  = »HAODT*DOS(THA) *SIN(PSA) +THAD0r*C0S(P3A) 

HAZ  = -®HADDT* SIM( TH A)  fPSAOOT 
C 

FIND  VE.DCITY  DF  TARGET  RELATIVE  TO  AIRCRAFT  IN  AIRCRAFT  FRAME 
C (EARTH  AXES) 

C 

XTADA  = XrADI-WAY*ZATMAZ*YAT 
YTADA  s YTADI-WAZ*XAT*WAX*ZAT 
ZTAOA  = ZTADI-WAX*YATMAY*XAT 


TRANSFORM  VELOCITY  OF  TARGET  RELATIVE  TO  AIRCRAFT  TO 
3IM3AL  AXES 


YTAOAG  * CET,*XTADAOST5*YTADA+CST9*ZTADA 


TTAOAC 


EK’XTAOAOET&’YTAOAfCEn’ZTAOA 


C*******  C0MOJT;  SHOAL  ANSl*  RATES  (GIM3AL  AXES) 
°SGAOOT  = YTAOAG/(RANAr»CTHGA) 

THGAOOT  = -ZTAOAG/RANAT 


C 

%*******  TRANSFORM  ANGULAR  velocity  OF  GIM3AL  frame  relative  to  aircraft 
C FRAME  (MSA)  TO  EARTH  AXES 

c 

WG  AX  = -CISl*PSGA3DT*STHGAfCIG2*rHSA3DT+’0IG3**  SGA33T*  CTHGA 
HGAY  = -CrS  + »PSGAOOT*STHGA*-CIG5*THSAD3T+SrS5*s3SAO)T*3THGA 
WGAZ  = -CI3Z*PSGA03T»5rHGA*CIG8*rH5A30TOIG9»33CA03T*CTHGA 


c* 

******  FI 

NO  ANSJLAR  VELOCITY  OF 

GIMBAL  FRAME  RELATIVE  TO  EARTH  PRAMS 

c 

(HGI ) (EARTH  AXES) 

1/ 

HGI  X 

= MGAX+MAX 

WGIY 

= HGAYfWAY 

WGIZ 

= HGAZfHAZ 

• 

c 

C* 

FIND  ANGULAR  VE.OCITY 

OF  LINE  OF  SIGHT  HLOS)  (EARTH  AXES) 

t 

NX  = 

HSIX 

NY  * 

iHSIY 

HZ  = 

HSIZ 

• 

\0 

c*******  FINJ  o3sitION  OF  IISSLE 

RELATIVE  TO  ATRCRA rT  (EARTM  AXES) 

RAMX 

s XE  - XEA 

RAMY 

s YE  - YEA 

RAMZ 

= ZE  - ZEA 

C 

C****»*»  FIND  VE.OCITY  OF  HI3SLE  RELATIVE  TO  AIRCRAFT  IN  EARTH  FRAME 
C (EARTH  AXES) 

C 

VAMX  = XE3  - XEAO 
VAMY  = YEO  - YEAO 
VAMZ  = ZEO  - ZE  AO 
C 

c* ******  pIN0  DISTANCE  of  MI3SLE  NORMAL  TO  -INS  OF  SISHT 
C 

OY=Cm*RAM<  f CET5*RAMY  + CETG’RAMZ 
OZ*CETZ»RAMX  «•  CET3*RAMY  ♦ CET9*RAMZ 
C' 

c*******  FIN5  VE-OCITY  OF  MISSLE  NORMAL  TO  LINE  OF  SIGHT  (EARTH  AXES) 
C 

OELXO  = VAMX*WZ*RAMY-WY*RAHZ 
OELYO=VSMY-HZ*RAMXfM<*RAMZ 
OELZO  = VAMZ*-MY*RAMX-W<*RAMY 
C 

c*******  TRANSFORM  VELOCI T t of  MISSLE  NORMAL  TO  LINE  OF  SIGHT 
C TO  AIRCRAFT  GIM3A.  AXES 

C 

9VY*crrA*OELxo  ♦ :et»*oelyo  *■  cets*oelzo 

OVZ=CETZ*OE LXO  CST)*3SLYD  ♦ CET9*OELZO 

c 

C*»***,*  FIN0  FROM  NCS3LE  TO  TARSEr 

C 


( 


RANG:  = S3RT((xeT-XF|**2MYeT-YE>**2MZET-7E)**2) 

C 

C*******  PINT  V E.OCITY  3F  HISSLE  RELATIVE  TO  TARGET  (EARTH  AXES) 

C 

VTMX  = X'O-XETO 
VTMY  = YE3-YETO 
VTHZ  = ZED-ZETO 
C 

COM3 JTE  RETRAM  FRiOJENCY  ERROR 
C 

IF( X LT • 2. ) 30  TO  2 
ALAM  = .1 

SIGI1  s ATAN2<SIGN(53RT(YAT**24ZAT**2> ,YAT)  ,<AD 
3IGI2  * AnN3(SIGN(S0RT(RAMY**2mHZ»»2>  ,RAMYI  ,RAM<) 

SIGIT  = 4TAM2(3IGS(S3Rr(YEAO**2tZEAO»*2>,YEAD) ,XEAD> 

EPS II  = ATAM2(SIGN(S3Rr ( ( YET- YE)* *2  + ( ZET-ZE) * * ?) , Y ET- YE) ,X£T-XE> 
EPSI2  * ATAN2<SIGM(5)Rr ( YED** 2+ZED* *2) ,YE0) ,XE3> 

SIG1  = A9S(SIGIl-SIGt3> 

SIG2  = A3S( SISI2-SI3I3* 

EPS1  = ATS(E3SIl-E3St2» 

•PS2  = A9S(E3SI2-3IGI2» 

OELF  = 49S(JA*(C05(St31> -COS ( SIG2) ) ♦VEL*(C0S(E352) *COS ( £PS1) ) ) / 

1 ALAM 
C 

C**** * **  COM3 JTE  HEADING  3-  3IM9AL  X AXIS  (EARTH  AXES) 

C 

2 e>SGAP  s A3IM(SPSG4*C3S(PHA)/COS(THA)) 

MGA  * 3SG43«-3SA 
HGA9G  s H3A  *57 . 295 
C 

C*******  TRANSFORM  VELOCITY  OF  MISSLE  RELATIVE  TO  TARGET  TO  GIM8AL  AXES 


C TO  ?IN3  A3P ROXI M A T E CLOSURE  RATE 

C 

VC  = C:Ti*VTNX*CET2*/TMV4-CET3*VTMZ 
C 

C**‘**‘*  FIN3  APPROXIMATE  HISSLE  TIME  TD  G3 


TTG  * R4M3E/49S (VC) 

C 

C*******  COM3JTE  ACCELERATION  COMMANO  COEFFICIENTS 
C 

IF(TTG.LT.  TR/4.  ) TTG  = TR/4. 

RR  = 1000. 

CH  = CT* (TR*RR/ RANGE) **2 
Cl  = • 01r»5*(TR/TT3)*»2 
C?  = ZT* • 125* ( TR/TTG) 

GGBIAS  = 1. 

CB  =(2.*MZ*VEL/32.2-RANCE*V00T*WZ/(32.  2*VED) 

C 

C*******  ACCOUNT  FOR  GUIDANCE  DELAY 
C 

IF ( (X  5T..7)  OR. (ACX.LT. 9.5))  GO  TO  5 
AH  * 0. 

AV  = 0. 

GO  TO  100 
C 

c*******  IF  MISS.E  is  AMEA3  OF  LOS,  COMMANO  TURN  TO  PARALLEL  LOS 
C 

r>  IF(  IGUIO  Fj,  x>  GO  TO  15 
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AHT  = CMMMSA-PS) 

AXI  = -AHT'SKKPS) 

AYI  = AHT*CDS(PS> 

• axs  = ceti* axi+cet2*ayi 

AHG  = CETKAXI+CET5MYI 
AVG  = CC1*DZ»C2*0VZ) ►G‘9IAS 
C 

C*******  TRANSFDRN  ACCELERATION  COMMANDS  TD  MISSILE  30DY  AXES 
C 

AH  = CM34*AXG+CMG5*AM3*CMGS*AVG 
AV  = -CMGZ*  AX5-CMSB*  AH3+CMG9* AVG 
C 

TRANSITION  TO  COMMAND  TO  LOS  . 

C 

AHTEST  = (-Ct*0Y-C2*DVY) 

IF(HGA*AHTE3T. 3T.0.)  3D  TO  100 

IF (X. GT. 2. ) IGUIO  si  * 

C 

• c*******  COMMAND  ACCELERATION  TO  LINE  OF  SIGHT 
C 

10  AHG  = (-C1»DY-C2*DVY*C3) 

AVG  = (C1*DZ+C2*DVZ) +G39IAS 
C 

C*******  TRANSFORM  ACCELERATION  COMMANDS  TO  MISSLE  300Y  AXES 
C 

92  AH  = CMS5*AH3-CMG3*AV3 
AV  = -CMGS'  AHG  + CM39MV3 
C 

100  RETURN 
C ' 

END 

SUBROUTINE  CONTROL 
C 

C*’ 

C*  THIS  ROUTINE  SIMULATES  THE  ACCELERATION  AUTOPILOT  AND  CONTROL 
C*  SURFACE  DYNAMICS  AND  DETERMINES  THE  EF-qtfTIVE  CDNTRO.  SJRFACE 
C*  OEeLECTIDNS  (D3,  OC>  3IVEN  THE  ACCELERATION  COMMANDS  (AV,  AH). 
C* 



C 

COMMON  / CMD/  AV,AH,WY»<(Z 

COMMON  /CNT.12  OB, DC . DT , ACZ . ACY ,ACX , 3BI AS » GH3I AS, I JK 
COMMON  2CNTL2/  R, 3, G , U, C , B1 , PI ,P2 , X , 3 
C 

C* 

C*  VERTICAL  CHANNEL  AM3LIrIER 

C* 

C*  C = B(  3. 11(.  03  S +1)2(0  1 S + 1 ) ),  -HERE  S > JH 
C* 

C 

C»***»”  TOTAL  CDMMANOEO  VERT  ICAL  ACCELERATION  (LIMITCO  TO  20  G»S) 
C 

AsAV+GRtAS 
IF(A  GT. 20. ) A=20. 

IF(A.LE.-20.)A=-20. 
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c 

AOO  VERTICAL  ACCELEROMETER  DAMPING  COMMAND 
C 

9 * A f ACZ 
C 

C*******  ACCOUNT  FOR  GUIOANCE  OELAY 
C 

IF(X.LT.  Z.  DR. ACX.CT.3  5>B=0.0 
C 

C*******  COM*UTE  VERTICAL  CHANNEL  OUTPJT  (.IMITED  TD  3 DECREES) 

C 

91M31  1*3“10«  *C)  *OT  ► B t 
C=B1*  033*3 
IF(C  Sf. 11. 1>  C = U.l 
IFCC.LE .-11. 1) C=-ll. t 
C 

C»»***«*  AOO  »ITCH  BODY  RATE  DAMPING  COMMAND 
C 

0 * C -ll.A5*Q 
C 



C* 

C*  HORIZONTAL  channel  amp.ifisr 

c* 

C*  P2  = 0(  3. 11(. 03  S »l)'(D.l  S ♦ 1 ) ),  WHERE  S = JW 
C* 

C 

C*******  TOTAL  CDMHANOED  HORIZONTAL  ACCELERATION  (LIMITED  TD  20  G'S) 

C 

H=AH 

IF(H.GT.?0.)H=20* 

IF(H  LE.-20, )H=-20. - 
C 

C***»*»*  AOO  HORIZONTAL  ACCELEROMETER  OAMPING  COMMAND 
C 

0=H-ACr 

c 

C‘****»*  ACCOUNT  ?0R  GUIDANCE  OELAY 

c 

IF(X.LT..T.DR.  ACX.GT.3  5)0=0. 0. 

C 

C*»****»  COM’JTE  MORIZONTA.  CHANNEL  OUTPJT  (LIMITED  TO  3 DECREES) 

C 

Pis (31. 1 * D* 1 0. *P2)*0Tf»l 
P2*P1». 933*0 
IF(P2.3T.ll.l)P2=ll.l 
IF(»2,.E  -11.  l)P2=-ll.t 
C 

C**»»***  AOO  YAW  9D0Y  RATE  DAMPING  COMMANO 
C 

RP=P2-ll. Aa*R 

c- 

(!»**♦•»*  AOO  HORIZONTAL  CROSSOVER  ANO  FEE03ACK  TO  VERTICAL  CHANNEL 
C 

£a-G+R3»0 

C 

c*******  COMMUTE  VERTICAL  CONTROL  SURFACE  RATE  (LIMITED  TO  300  DECREES  ' 


I 

I l. 


i 

i 

I 

I 
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C 


c PER  SECOND) 

C 

FF=66.**E 

IF(FF. “T.  300. > FF=300. 

IF(FF.,E.-300. >FF=-300 
C 

C*******  INTEGRA r E VERTICA.  CONTROL  SURFACE  RATE 
C 

G=FF*OTfG 

C 

C*******  AOO  VERTICAL  CROSSOVER  AND  FEED3AC<  TO  HORIZONTAL  CHANNEL 
C 

Sl=-0+RP-'J 

C 

C** *****  COH’JTE  1DRIZOMT  A.  CONTROL  SURFACE  RATE  (.IMITEO  TO  300  DEGREES 
C PER  SECONO) 

C 

7=66.7*51 

IF(T.GT.300,)T=30Q. 

IF(T  LT.-300. > T=-300. 

C 

C*»*»*»*  INTEGRATE  HORIZONTAL  CONTROL  SJRFAGE  RATE 
C 

U=T*OT*U 

C 

C»**»***  CALC JL A T E EFFECTIVE  VERTICAL  CONTROL  SURFACE  OErLECTION  (03) 

C (LIMITED  TO  15  DECREES) 

C 

OB=(-G*U)/Z. 

IF (03. 3T  15. >03=15. 

IF(03..T.-15>0B=-15, 

C 

C*“***‘*  CALCULATE  EFFECTIVE  HORIZONTAL  CONTROL  SURFACE  DEFLECTION  (DO 
C (LINITEO  TO  15  DECREES) 

C 

oc*(u*;>/2. 

IF(OC.CT-i;.)OC=15. 

IF(OC  LT.-15. > OC=-15. 

c 

RETURN 

C 

’ END 

• SU3R0UTINE  E'JLER(I"LA5> 

C 

c*  * 

C*.  THIS  ROUTINE  COMPUTES  THE  AIRCRAFT  GIM3AL  ANGLES  PSGA  AND  THGA  * 

C*  (ASSUMING  ®ERFECT  TARGET  TRACKING)  AND  COMMUTES  THE  TRANSFORMATION  * 
C*  MATRICES  HHICM  RELATE  THE  EARTH  AXES,  AIRCRAFT  AXES,  31 19 AL  AXES  * 
C*  AND  MISS. E AXES.  * 

C*  * 

c 

COMMON  /ACC/  AX9tAr8,A73, AXE, AYE, «ZE,POOT,aDOT, ROOT 
COMMON  /AN3/  PS,  TH,  »H,  PSDG,  THOG,  °HDC  ,PSA,  T1 A,  >1A, 

1 PSAOG,  TH9G,  PHAOG,  ®SCAOG,  THGADG,  1GADG 

COMHON  /MATRIX/  CET 1 , C • T2, CET3, CE \<* , CE T5,  GETS,  CET7, CE T 8 ,CE T9 , 

1 CICl,CI3Z,CICJ,CIGi»,CIG5,CIC6,0IG7,0ICi,CtC3, 


:.jt6 


f 


I ■ 


i :h31,2.i,;2»:i:;i,:>ig<*,cmgs  15*06, 

i s®s;A,CPsr,A,sn;A,cTM3A 

COMMON  /®0  3/  XET,YET,7-T,xeA,YEA,ZEA,RANAT,X£ ,YF,ZS 
COMMON  /VELO/  X£0,ri),?i0,XEA0,r£A0,Z£40,J0,Y3,IO 
COMMON  /VTARSET/  XETO.f  ETO,ZETD,OTD,VDTO,Z3T3 
REAL  LIP,L?9,L3P,N1P,M?P,M3P,N1P,N2P,S3P 
REAL  L1,L2,L3,M1,M2, 13, N1,N2,N3 


C 

C* 

C 


IF(IFLAG.E5.2)GO  TO  20 

>♦*  TRANSFORMATION  FROM  EARTH  TO  AIRCRAFT  AXES 


LIP* 

L2P* 

L3P= 

Ml9* 

M2P* 

M3P* 

NIP 

N2?  = 

N3P 


COS(THA) *COS(PSA) 

SIM ( PSA ) *COS(TMA* 

-SIN(MA) 

COS (PSA)  * SIN (THAI  * 3 I M( PHA) -SI N( PS A ) *OOS(9HA) 
SIM ( PSA  > •SIN(TMA)  *?IN(PHA)  ♦ COS(  9SA)  *COS  ( 9HA) 
CO?(THA)‘SIN(®HA) 

COS(»SA) *SIN(THA) *ODS(PHA> ♦SIN(P3A)*SIN(PMA) 
SIM  (PS  A)  • SIN  (MAI  •:  OS  (PHA) -COS  (PS  A)*  SIM  (PHA) 
COS(THA)*OOS(PHA) 


c 

TRANSFORMATION  FROM  EARTH  TO  MISSILE  AXES 
C 


LI 

L2 

L3 

Ml 

M2 

M3 

N1 

N2 

N3 


COS(H)*COS(PS) 

SIN(PS) *COS(TH> 

-SIN(T-t) 

COS (PS) *SIN(TH) * SIN (PH) 
SIN(9S)  *S  IN  (T-0  *5  r M (PH) 
COS  f TM) ‘SIN(PH) 

COS(®S) *SIN(TH> *00S (PH) 
SI M ( 95 ) * SIN (TH) *003  (PH) 
COS(T-I)  *COS(P-() 


SIN(PS) *C9S (PH) 
COS(PS) *C05(PH) 


SIN (PS) * SIN (PH) 
COSCPS) *SIN(PH) 


C 

c*»** 

c 


c 

c**»* 

c 


***  find  RANGE  FROM  AIRCRAFT  TO  TARGET 
RANAT=SORT( (XET-XEA) **2  ♦ (YET-YEA) **2  ♦ (ZET-ZEA) **2> 

***  COM9 JTE  AIRCRAFT  3IM9AL  ANGLES  (LIMITEO  TO  30  DEGREES) 


C 

C 

c 


THE  COMPUTATION  0"  CTHSA  AND  CPS3A  IS  VALID  ONLY  Ir  TNSA  IS 
BETWEEN  -33  DEG  AND  *■  90  OEG 


STHGA*-(N19*(XFT-XEA)*N2PMYET-Y£A)+M3P*(ZET-ZEA> ) /RANAT 
IF(STH3 A ■ 3T . 1. ) STH5A*.  3999 
IF(STHSA  LT.-l. ) STMS \=- . 9999 
CTHGA=SORT ( 1 -STHSA* 3TM3A) 

THGADG  * ASIN(STHSA) *5’. 295 


SPSGA*1. /( RANAT’CTHGA) * ( (XET-XEA) *Ml9f ( YET- YEA) *M29f ( ZET-ZEA) *M3P> 
IF(SPSSA.3T.l. ) SPSS A *.9999 
TF(SPSSA  LT.-l. )S933A=-, 9999 

CPSGA  * (L  I'P*  ( XET-XE  A)  *L2P*  (YET-YEA)  t«_3P*  (ZET-ZEA)  ) /(  OTHGA*  R AN  AT ) 
PSGAOG  = ASIN(SPSSA) *57, 295 


C 

c**** 

c 


***  TRANSFORMATION  FROM  AIRCRAFT  AXES  TO  3IM9AL  AXES 
CAT 1*CTHGA*39SGA 
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r 


» -i 


m 


' r*sj* 

• 

!.v» . 


,-3< 

H 


c 

c**“* 

c 


c 

c 


c 

c**»* 

c 


CATI?  = Crtr.ft»5:,SCA 
CAT3*-STHGA 
C AT4=-SPS3A 
CATS^SGA 
CAT6=0.0 

CAT7=srH3A»C®SGA 
CAT3=STMG4»  SPSGA 
CAT9=CTHGA 

**«  TRANSFORMATION  FROM  EARTH  AXES  TO  GIM3AL  A XES 

CETl=CeU*„lP*CAT2*MlPvCAT3*NlP 
CET2=C»n*.?3fCAT2»H>3f:AT3*N2P 
CET3=CAri*.3P*CAT2*M3?*CAT3*N3P 
CET4=OAn.*.ie*CAT5*MlPfCAT5*NlP 
CET5  = CAr"4*„2P«-CAT;*H?3t-CAT6*N2P 
CETB  = CAT4».3?»CAT5*M3Pt-CAT6*N3P 
C£T7=CAT7*.  1P*CAT3MIP4-CAT9*N1P 
CET3  = CSrr»„?:>*CAT8,M?PfCAT9*N2P 
CET9=CAT7*^33*CAT8*M3?fCAT9*N3P 

♦♦♦  TRANSFORMATION  FROM  AIRCRAFT  AXIS  TO  EARTH  AXES 

CIA1  = COS  ( ®3A)  *COS  ( Ml ) 

CIA2  = -C3S(3SA) *31 N ( T < A) *SIN (PHA) -STM (PSA) *33S(®HA) 
Cl  A3  = -C3S(PSA)*StN(TMA)*COS  <PHA ) frSIN  (PSA)  *3IN(PHA) 
CIA4  = SIN ( ®3A) *C3S(rHA) 

CIA5  *=  -SIN(PSA)*SIN(TNA)*SIN(PHAH-C3S(P3A)*:3S(PHA) 
CIA6  = “SIN (°SA)*SIN(TNA )*COS (PHA) -CDS (P3A) *SIN(PHAl 
CIA7  = SIN(TMA) 

CIA9  = C03(THA)*SIN{5Hn 
CIA9  = COS ( TH A) *C3S ( >4A ) 

**»  TRANSFORMATION  FROM  SIM3AL  AXES  TO  AIRCRAFT  AXES 


CAS1 

S 

CPSSA*CTHGA 

CAG2 

= 

-SPSGA 

CAG3 

r 

CPS3A*STHGA 

CAS  4 

s 

S°3  3 A *3  THG A 

CAG5 

= 

C°33A 

CAG* 

r 

S®S3A*STH3A 

CAS7 

= 

-ST  M3  A 

CAG3 

s 

0.0 

CAG9 

= 

CTH3A 

C 

c*«»» 

c 


c 

c»»»» 


*♦»  TRANSFORMATION  FROM  GIM3AL  AXES  TO  EARTH  AXES 

CIG1  = DI41*C431*3I4>»;AG4»CIA3*0A37 
CIG2  = CIAl*CAG2*CIA’*CAG5+CIA3*CA33 
CIG3  * 3IAl*0AS3*0IA?*:AGG«-CIA3*3A39 
CI34  = CIA4*CA31*CIAJ*;aS4»CIA6*3A 37 
CIG5  = CIA«»CAG2+CIA»»CAG5fCIA6*CA38 
CIG6  = 3IA4»CAG3'»CI4>*CAG6+’CIA6*3A39 
CIST  = CIA7*0AG1  + CIAJ»08G4«-CI  A9*CA37 
CIG8  = CIA7*CAG2*CIA3*:AG5»CIA9*0A38 
CIS9  = CIA7*CAS3«OIAJ*CAG6»CIA9*CA39 

♦*»  TRANSFORMATION  FROM  GIMBAL  AXES  TO  HISSLE  BOOT  AXES 
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c 

omgi  = .t*:i;i»L2»:x;v»L3«szc7 
CMG2  = L1*CIS2+L2»CI35*L3*CIG8 

cmst  = .i*cis3*l2»ci;3*l3*cis9 

CMG<*  = Ml*CI31*M2»CIS4*M3*CIG7 
cmg5  = 'U*;ir,2tM2*:r;;f  H3*ciG8 
•CMG&  = mi*ci;3*mz*ci;5*m3*cig9 
CMG7  = Nl*CIGt*N2»CI3i,fN3*CIG7 
CHG9  = Ni*ci32tN2*oi;>*M3*ciG8 
CHG9  = N1*CIS3*N2*CIS3*N3*CIG9 

c 

C**.**»*  TRANSFORM  VELOCITY  COMPONENTS  IN  EARTH  AXES  (XE 3 , Y ED , ZED) 

C TO  MIS3.E  AXES  <U3,V3,W9) 

c 

UB  = Ll*XE3  * L2*YED  L3*ZED 
VB  = M1*XE3  v H2*YE3  ♦ M3*ZEO 
HB  = Nl*XE3  * N2*YE3  * N3*ZEO 
C 

RETURN 

C 

C»***»«*  TRANSFORM  ACCELERATION  COMPONENTS  IN  MISSILE  AXES  ( AXB,  AY9,  AZ3) 
C TO  EARTH  AXES  (AXE, AYE, AZE) 

C 


20  CONTINUE 

AXE  = ,1*A<3  M1*AY3  * Nt’AZB 
AYE  = .?*AX3  f M2*  A Y 3 * N2*AZB 
AZE  = . 3* A X3  + M3*AY3  f N3*AZ8 


C 

C**”**»  ADO  ACCELERATION  3JE  TO  GRAVITY 
C ' 

AZE  = AZE  ^ 32. 2 
C 

RETURN 

C 

END 

SUBROUTINE  INTLFLT 
C • 

c* 

C*  THIS  ROUTINE  SETS  FLISHT  PARAMETERS  TO  THEIR  INITIAL  VALUES. 

C* 

c 

COMMON  7AEJLER/  PHAD3T, THADOT ,PSAODT 

COMMON  /AMS/  PS,  TH,  “H,  PSOG,  THOS,  »HOS  ,«>SA,  THA,  “HA, 

1 PSAOG,  THAOG,  PHAOG,  =>S3AOG,  TMGACG,  HGA3G 

COMMON  /BL3C<1/  XXI, YYT ,7ZI , AMASS, CG, MG, “BIAS, XX 
COMMON  /3L3C<2/  AHS,  A VS , TT5 ,0ELF, ISUID , TR.UA,  RANGE,  ZI , CT, VD3T 
COMMON  /CH3/  AV,AH,HY,HZ 

COMMON  /CNTLl/  08, OC , DT, ACZ , ACY,ACX , S3 IAS,GM3I AS, I JK 
COMMON  /CNT.2/  R,0,  S , J,  C , 31 ,“  1 , P2,X  , “ * 

COMMON  /OER/  0,RL,CMA3,CMCG,QQ,S, VEL , T HR, CMA3 Y ,CMC SY,  RY 
C01M0N  /PIN/  0000,  O-VM 
COMMON  / 1 N [ T L / H , T I M i,  IT  A8L,IPR,ILES 
COMMON  /P33/  XETfVET,ZET , XE A, YEA, ZEA,RANAT ,XE,YE,ZE 
COMMON  /VE.3/  XEO,  YE  3,  l EO, XEAO, YE A3, 1 1 AO,  JO,  V3 *W8 
C 

c 
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C*******  CONVERT  AIRCRAFT  ATTITUDE  FROM  DECREES  TO  RADIANS 

c 

THA  = THAOS  f 57.235 
»HA  = BRADS  / 57.235 
PSA  = PSADS  / 57.235 
C 

C”V*‘»  SET  INITIAL  GUIDANCE  FLAG 
C 

IGUIO  * 0 

c 

C* *»*♦**  INITIALIZE  MISSILE  PARAMETERS  TO  AIRCRAFT  PARAMETERS  AT  .AUNCH 
C 

TH=THA 
PS=PSA 
PH=PH A 
XE=XE A 
YE=YI A 
ZE=ZEA 
XEO=XEAO 
YED=YEAD 
ZEO=ZEAD 
VEL  = DA 
C 

c*.4*»*»  INTIALIZE  ACCELERATION  COMMANDS 
C 

AH  = 0. 

AV  = 0. 

C 

c. SET  INITIAL  VALDES  DM  INTEGRATORS  IN  CONTROL  3U3R0JTINE 
C - 

OT  * H 

Ml  * 0. 

B2  = 0. 

PI  * 0. 

P2  * 0. 

C * 0. 

G a 0. 

U « 0. 


SET  INITIAL  VALUES  OF  MISSILE  AVGJLAR  VELOCITY 


P » 0. 
Q a 0. 
R » 0. 


SET  INITI4L  VALUES  OF  AIRCRAFT  EU.ER  RATES 


PHAOOTaO. 

TMAOOTxO. 

PSAOOT«0. 

C 

SET  INITIAL  MASS  PARAMETERS  ANO  REFERENCE  AREA  OF  MISSLE 
C 

AMASS  * 300. /32. 2 
XXI  a 2.5 
VYI  a 57.0 
ZZI  a 57.0 
S » .525 
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[ j! 

f 


c 

C** *** •*  SET  INITIAL  VALUES  3N  ACCELEROMETER  OJTPJT5 
C 

ACX  3 0. 

ACY  = D. 

AC Z = 0. 

’ VOOT  * 0. 

C 

C»****«  INITIALIZE  TIHE 
C 

X = 0. 

c 


c 


c 

c** 
c* 
c* 
c* 
c* 
c* 
c * 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c * 
c* 
c* 
ft* 
0* 
c* 


RETURN 

END 

SUBROUTINE  NDTLU (ND, NA , X , Z, XA , ZRt IE » NEXTR) 


N-DIMENSIONAL  TABLE  .OOK-UP 

GIVEN  THE  ARGUMENTS  < C 1 > , X( 2) X < N*l>  , 

30NPJTES  Y = F(X<1),X(2),...X('I-1>)  8Y 
.INEAR  INTERPOLATION  FROM  A TA3LE. 

ACCURACY  REAL  SINGLE  PRECISION 


CALLING  SEQUENCE 

CALL  LOOKJ»<ND,NA,X,  Z,  < A ,ZR,IE,NEXTR> 

WHERE 

NO  * DIMEN3I ON  OF  LOOK-UP 
HHEN  f a c (X)  NO  = 2 

MA  = AN  ARRAY  OF  LENGTH  N3-1 

NUMBERS  OF  VALUES  OF  EACH  TA3.E  OF  X 

THE  TABLES  ARE  LISTE3  BY  SIZE,  THE  LARGEST  rIRST • 


C* 

C* 

C* 

ft* 


X a TABLES  OF  EACH  X IN  ORDER 
Z » FJNCTI3N  VALUES 

IF  A a F (X, Y, Z)  THE  3EPEN3ENT  VARIABLE  ARRAY 
MUST  BE  IN  THE  FOLLOWING  ORDER, 


c* 

ASSUME  X = 4»  4=  3, 

Z=2 

c* 

ZCl)a 

F (XI , T 1 , Z1 ) 

Z ( 1 3)  a 

F(X1,Y1,Z2» 

ft* 

Z ( 2)  a 

F(XZ,T1,Z1) 

Z(i*)  = 

F(X2,Y1,Z2> 

c* 

Z(3>* 

F(X3, TL, Zl) 

Z C 1 5 > a 

-CX3,Yt,Z2) 

c* 

Z<0» 

F(X3, f 2, Z 1) 

Z(13) = 

F(X3, Y2,?2I 

c* 

Z ( 9 1 a 

F (X4, T2, Z 1) ’ 

Z<  201  a 

F(X4,Y2,Z2i 

c* 

ZC9I* 

F (XI , T3 , Z 1) 

Z(21)= 

F (XI » Y3  »Z2) 

c* 

Z(10)= 

F(X2,Y3,Z1> 

Z ( 22) = 

-(X2, Y3,Z2> 

c* 

mu- 

F(X3, Y3.Z1) 

Z(23)  = 

F(X3,T3,Z2i 

c* 

Z(12)a 

F(X4,Y3,Z1) 

Z(24)a 

F (X4, Y 3 ,Z2> 

c* 

c* 

ZR  « 

RSSU.T 

c* 

c* 

NEXTR  a 

1 SKTRAPOLATE 

c* 

■ 

0 N3  EXTRAPOLATION 

J 


! 


< 
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c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c 


6 

a 


3’ 

10 

u 

12 

1- 


13 

11 

3 


22 

23 


IS  * EWO)  :ooe 

0 N3  ERROR 

-1  X ARRAY  TOO  SiflLL 

1 X ARRAY  TOO  LAR3E 

2 ARRAY  NOT  IN  ASCENDING  ORDER 


OIMINSIOM  X (1) ,2(1) , NA  (1) ,XA (1) ,NS(5>, WJ(32I 
EOUIVA.SNDE  ( XMM, NM)  , (XMH.MH) 

LI =2 
LFsMO-l 
00  3 1*1, L- 
L2SL1*'JA  (I)  -2 
F0UN0=3, 

00  4 J*L1, L? 

IF<X(J>.GT.X(J-i>>  GD  TO  6 

IE=2 

RETURN 

IFCFOUNO-  NS.O.  ) GO  TD  4 
IF<XA(I)-X<  J-Dl  3 , 4 , V 
IF<  J 3T.L1>  30  TO  10 
IF(NEXTR.ia.O)  GO  TD  32 
FOUND*!. 

NS( I ) =.1-1 
GO  TO  4 
IE=-1 
RETURN 
FOUNOsl. 

NS<I»*J-2 

C0NTINJ2 

IF 1F0UND)  11,12,11 

tF(XA(I»-X(l.?n  13,13,14 

IFINIXTR.NE.O)  GD  T)  13 

IE*1 

RETURN 

NS  II) *.2-1 

L1*L2*> 

CONTINUE 

KF*2**LF 

HH«-2 

00  13  1*1, <F, 2 

Lt*0 

IZ*0 

NPT*1 

00  19  U*l,.s 
HM*2»*CJ-l» 

IF( AND (XMM,  XT<<)  • ED.  0>  30  TO  22 
N*NS(JI»1 

GO  TO  23  - 

N=NSCJ> 

U*N-L1 

L1*L1+NA(J) 

IZ*NPT*fN-miZ 

NPT*NPT»NA(J) 
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WJ(I>=MIZf  L) 

10  WJ(IM)  :Z(IZ»2) 

00  32  T=t,LF 
H = MS  Cl) 

KF=KF/? 

RATI0  = (XA<I)-X<M))/<<<>1H>-X<M>> 

00  3?  J*l,<r 

32  MJ(  J)  =WJ(2» J-l) ♦<WJ(2*J) -HJ(Z*J-i) ) 'RATIO 
ZR=WJ  (t) 

C 

RETURN 
C 

ENO 

SUBROUTINE  OUTFLT 
C 

C* 

C'  THIS  ROUTINE  PRINTS  HE  OUTPUT  DURING  HE  "LIST  T 
C*  ‘ 



c 

COMMON  / AER/  OEL»DSC»S'll»CllllCC»CK12i3N2»C‘1,S03,OELSOfft[1ACHf  ALPHA 

l , alt,  beta,  3CY,3M2Y,3N2Y,CMY,C'iiY,:MiY,oocY,;o,:N,CNY 

COMMON  /an;/  PS,  H,  PT,  P5DS,  THD3,  PHOS  ,*SJ,  HA,  A, 

1 PSAOS,  HA3G,  PHAOG,  3S3A0G,  THSA3G  , HGAOG 

COMMON  /0.03<2/  AHG, AYS, TTG , DELF, ISUI3 , TR, U A, RANGE, ZI,CT,VD3T 
COMMON  /CMO/  AV,AHHY,TZ 

COMMON  /CNTL1/  OB, 03 , 3 T , ACZ , ACY ,ACX , 331 AS ,GM3I AS ,1 JK 
COMMON  /CMTL2/  R, 3,  G , U, 3 , 01 ,P1 , P2 , X , * 

OOMMON  /OER/  3 , RL,  CM  A 3,  3MCG  ,QQ  ,S,  VEL  , HR,  CMA3Y  , CMO  3 Y,  RY 
COMMON  /°3S/  XET,YET,ZiT,XEA,YEA,ZEA,RANAT,<E, YE,ZE 
COMMON  /PRMT 1/  DY,DZ,OYY ,DVZ 
C 

CONCERT  MISSLE  ATTITJOE  FROM  RA3IAMS  TO  3E3REES 
C 

PS0G=5T. 235»PS  . 

TH0G=57.  235*TH 
PH0Gs5T.235'?H 

C ^ 

C*»4»»M  OUT* JT 
C 

WRITE (9,000  ) 

WRITE  (S,0 1)5!  X,TTG 

write <9,<uo>  vel,ama:h,oelf,iguio 

WRITE (5, 015)  OY.OVY, AT3, AH,3ETA  ,DC, A3Y 
WRITE (5 , 92D ) OZ.OYZ, 1Y3, AV, ALPHA, D3* AOZ 
WRITE (3,023)  XE,PS3G»X*A,PSGA0G»HGAD3 
WRITE(9,O30)  YE,THDG,YEA,THGAOG 
WRITE (5,0351  Z E.PM33 , ZE A ,WZ 
C 

RETURN 
C 

900  FORMA T( /) 

105  FORMAT ( IX,  5 H TIME  ,F0. 2, 3X,<»HTTG  ,F9,2> 

910  FORMATdX,  HYEL  ,F9,  2,  JX ,5HMACH  , F3 . 2, 3X, 5M3ELF  ,F3.2,  3X,SHI3UI3  , 
115) 

915  FORMA T(/1X,2H0Y»F11»  l,3X,4H0VY  ,"9. 2, 3X,3HAH3, fll, 2,3<, 3HAM  ,F8.2, 
• 13X,SH0E T A , F7 . 2, 3X, JM DC  ,F9. 2,3X, AHACY  ,F8,2» 

[ 

! • • • 


i 


! 
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920  PORIATt  1X.2  W.Fll.  l,3X,4M0VZ  ,r  9.  2 , 3 X , 3HA  V3 , e-11 . 2 ,3  < , 3HAV  ,FB.Z, 
13X,5HAL3HA,?  ).  ?,3X,.HDD  ,F1.2,3X,4HA3Z  ,F9.2I 
325  FORMA T( /IX,  MXE.Fll.  ?,3X,3HPS  , FI  0.  2 , 3X , 3MX  EA , Fll.  > ,3< , 4HPS3A , F7. 2 
1.3X,^HHr.A  ,F3.2) 

930.  FORMAT  C IX,  2MY  S , FI  t . ? , TX  , 3H  TH  , F10,  2»  3X , 3MYE  A , FI  1,  2 , 3<  , 5HTM3  A , 

1 F6i  2) 

935  FORMA T ( IX , 2HZ E.F11, 2, 3X , 3HPH  , FI 0.  2, 3X, 3HZEA , rll.  2 , 3< , 3HHZ  ,F8,2) 
C 

END 

SUBROUTINE  DJTINTL 
C 

C*  . * 

C*  THIS  ROUTINE  PRINTS  THE  INITIAL  INPUT  VALUES.  THE  A EROD Y NAMIC  * 

C*  TABLES  ARE  A. SO  PRINTED  Ir  ITABL  =1.  * 

c*  * » 

C 

COMMON  /AN3/  PS,  TH,  PM,  PSOG,  TMOG , PHDS  ,P3A,  THA,  3HA, 

1 PSA OG,  TMAOG,  PHADG,  3S3ADG,  TH3AOG  , H3ADG 

COMMON  /3„D0Kt/  XXI,  f f T , ZZI  , AMASS,  CG,  W3, 33I  AS,  XX 
COMMON  /9_0C<2/  AHG,  A V3 , TTG,  OELF,  1 3UI0  , TR,  J A,  R A N3E  , ZI,  CT,  V DDT 
COMMON  /CNTll/  OB, OC , OT, ACZ , ACY ,ACX , 33 1 AS.GH3I AS, I JK 
COMMON  /INI TL/  H,TIME , I TA9L,IPR,ILEG 
COMMON  /PDS/  XET, YET, ZE T , XE A , YEA, ZE A, RANAT, XE , YE, ZE 
COMMON  /VE. 0/  XEO, Y ED , ZED ,XE AO ,YEA3 , ZEA0»U8,V3,W3 
C 

C*******  PRINT  INITIAL  INPJT  VALUES 
C 

WRI TECS, 9000) 

8000  FORMAT(*’1‘*//40X,  "RET  RAN  MISSILE  SIMULATION'*) 

HRITE (3,3000)  XEA,  YEA,  ZEA,  XEAO,  Y'EAD,  ZEAO,  TMADG,  PSAOG,  PHAOS 
1 ,X*T,  YET,  ZET,  H,  33IA3,  GH8IA3,  TIME,  P3IAS  , CT  , TR  , ZI 

C 

c.»******  OPTIONAL  PRINTING  0-  AERODYNAMIC  TABLES 
C 

IF( IT A3L  ED,  1)  CALL  D'JTTBL 
C ' 

C****  * **  OPTIONAL  PRINTING  Or  VARIABLE  LEGENO 
C 

/'  IF(ILE3.  EQ,  1)  WRITE (3,7  0 00) 

IFdLES.ED.  1)  MRITE(  3, 7500) 

C 

RETURN 

C 

7000  FORMAT ( //L3X, "OUTPJT  L: GENO"/ 

i "TTMr  - r t so;rn  ci  f.jr  IT u:  i«i<i  « 


i ins.  - nun  1 1 ns  tbsoi  , 

1*'TT3  * TIME  REMAINING  TD  TGT  IM^AST  (SEC)  “/ 

1**VEL  - MISSILE  AIRS*EEO  (FT/SEC)  ", 

l'*MACM  - MISSILE  MAH  NUMBER  "/ 

1"PELF  - RETRANSMITTED  FREQUENCY  ERROR  (HZ)  "/ 

l**OY  - MDRIZONT AL  DISTANCE  ERROR  (*T)  ", 

f'DVY  - HORIZONTAL  VELOCITY  ERROR  (FT/SEC)  "/ 

l“OZ  - VERT ICA.  DISTANCE  ERROR  (FT)  ", 

l"OVZ  - VERTICA.  VELOCITY  ERROR  (FT/SEC)  "/ 

1 "AM 3 * HOR.  ACCELERATION  COMMAND  - 3IM1AL  r RAME  (G*S)  ", 

1“AH  - MOR.  ACCEL,  3 OMMANO  - MISSILE  CRAME  (5»S>  "/ 

1"AVC  - VERT,  ACCE.,  COMMANO  - 3IM3AL  FRAME  (3»S»  ", 
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IMI  W9$ 


n 


1"A  V 
1 "BETA 
1 "ALPHA 
1"0C 
1"0B 
l-ACY 
1"ACZ 
FORMAT  (" 
1 "XcA 
i "YE 
1"YEA 
1 ”ZE 
1”Z£A 
1 -ps 

1"PS3A  - 
i”TH 

1"TH3A  - 


JERT  . ACCE.,  COMMAND  - MISSILE  FRAME  (3*3) 
SIDESLIP  AN3.E  (OEG) 

ANGLE  Or  ATTACK  (DEG) 

EFFECTIVE  f A <1  CONTROL  DE-LECTION  (OEG) 
ErrECTIYE  * I T CM  CONTRO.  OE-'.ECTION  (DEC) 
MISSILE  ACCELEROMETER  OUTPJT  (3'S) 

MISSILE  ACCE.EROMETER  OUTPJT  (3'S") 

XE“,53X, 

MISSILE  POSITION  COORDINATES  (FT) 

AIRCRAFT  POSITION  COORDINATES  C - T» 


i"PS  - MISSILE  HEADING  (DEG) 

1**PS3A  - HORIZONTAL  3TMDAL  ANGLE  (OE3) 

i’*TH  - MISSILE  PITCH  ANGLE  (OEG) 

1-TH3A  - /ERTICAL  3 I M ? A L ANGLE  (OEG) 

1 "PM  - MISSILE  BANK  AN3LE  (DEG) 

1"HZ  - ROTATION  RUE  OF  LINE  Or  SI3HT  (RAO/SEC) 

1"HGA  - HEADING  OF  LINE  OF  SIGHT  (OEG) 

1"I3UIQ  - HORIZONTAL  A.GORITHM  MODE 

lft9X,“T3UID  s 0 S COMMAND  TO  HEADING"' 

169X,"I3JI3  = 1 J COMMAND  TO  TRACK"//) 

9000  FORMAT ( * INITIAL  F.I3HT  CONDITIONS*  / 

12X,4HXEA  ,F13. 2,3X, 4HYEA  , F 10. 2,3X, 4 MZ E A , F10  . 2, 3X , 3HXE AO  i 
i 3X, 5HYE AO  , "10 . 2,3X, JHZE AO  .F10.2/2X, 4HTHA  , -i 3. 2 , 3X, 4HPSA 
1 3X,  4H°HA  ,fio.  2,3X,!*M<rT  ,F  14.  2 , 3 X, 4 HY ET  ,F11.  2,  3X,  4HZET  ,f 
12X,2HH  , f+. 2,3X,6M33IAS  , F2. 0 , 3X, 7MG H9I AS  » FT . 0 , 3X  » 5M  T I ME  , 
1 3X, 6HPAI AS  ,F12. 2,3X, 3HCT  , F5. 2,3X, 3MTR  , F6. 2 , 4X, 3 HZI  ,F6,< 

% 

END 

SUBROUTINE  OJTTDL 


THIS  ROUTINE  PRINTS  THE  AERODYNAMIC  TA3.E  GENERATED  IN  SUBROUTINE 
AERO. 


COMMON  /AERTRL/  OCCT(5,3,6>,  CMT( 8, 5 ,3 ) ,CNT (3 , 5 i S) , Xl(19), 

1 CCT(10,3) , CMET (10,5) ,CN2T (10,6) ,CMTT (10, 3), <2(18) , 

2 COOT  (10)  ,MCHCOT(10)  ,H(a)  .OE-COKS) 

COMMON  /AMS/  RH02 ( 12 ) , S0S2 < 12 ) , ALTT ( 12 ) 

COMMON  /TMRST/  TIMES ( 2 l ) , THSTLB (21) 


HRITE (6.50) 

FORMAT ( IX/ / 44 X, "AERO DYNAMIC  TABLE"/) 
WRITECi,  Si) 

FORMAT(10X,*0CCT*) 

HRITE (3. 100) DC CT 
HRITE  (3,52) 

FORMAT ( 10X,*CMT*) 

HRITE (3, 103) CM T 
HRITE (3, S3) 

FORMAT ( 11X,*CNT*) 

HRITE (5, 10) ) CMT 
HRITE (ft, 63) 


66  FORMATdOX,  »Xl*> 

WRITE (>, 100) XI 
WRITE (St  5V) 

54  FORMAT (10X»*0CT*) 
WRITE (6, 110) OCT 
WRIT" { 6| 55) 

55  FORMATdOX,  »3M2T*> 
WRITE (5, 110) CM2T 
WRIT" (6,551 

56  FORMATdOX,  »:W2T*» 
WRITE(6,110)CH2T 
WRITE  <6, 52) 

57  FORMATdOX,  *3MTT*) 

writer, uo>:mtt 

WRITE (6*53) 

58  FORMATflOX, »X2*) 
WRITE(5,  110)  X2 
WRITE (5t53) 

59  FORMATdOX,  *TIMES*) 
WRITE(5,123) TIMES 
WRITE (8,60) 

60  FORMATflOX, *THSTL9*( 
WRITE (3,120) THSTLB 
WRITE (8,51) 

61  FORMAT (10X, *4LTT*) 
WRITE (5, 133 ) ALTT 
WRITE (5, 52) 

62  FORMATflOX, *S0S2*) 
WRITE (6, 130)SOS2 
WRITE  (5,53) 

63  FQRMATdOX,  *RH02*> 
WRITE (5,140) RH02 
WRITE (5,54) 

64  FORMATflOX, *300T*) 
WRITE (5,105) COOT 
WRITE  (6,55) 

65  FORMAT (10X,*MCHCOT*) 
WRIT* (5,135)MCHCOT 
WRITE  (5,57) 

67  FORMA T(  10X,  *'H*) 

WRITE (5,100)H 

• WRITE (5, 59) 

WRITE (6,55) 

68  FORMATflOX, ►OSLCOTM 
WRITE (5, 140)  OELCOT 

C 

RETURN 

C 

110  FORMAT  1 10e5 , 4) 

100  FORMA  T( JF10 , 4) 

105  FORMAT! 10E9  « 3) 

120  FORMAT(7E10.2) 

130  FORMAT( 5E10 , 2) 

140  FORMAT! 5F10, 5) 

C 

END 

SU9R0UTTNE  TARGET  (X) 
C 
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C*  THIS  ROUTPIE  SUPPLIES  TARGET  VELOCITY  IN  EARTH-FIXEO  INERTIAL 
C*  COORDINATES  AS  A FUNCTION  OF  TIME  ( X>  • 

c*  • 


COMMON  /VTARGET/  XETD,  YETO,  ZETO, XBTO,YDTO,ZBTO 

c 

c** ***■»•  target  is  stationart  in  this  simulation 
c 

XETO  =0. 

YETO  = 0. 

ZETO  s 0. 


RETURN 

END 

SUBROUTINE  THRUST(T) 


C*  THIS  ROUTINE  COMPUTES  THE  MASS  PARAMETERS  ANO  THRUST  DURING  BOOST. 
C* 


COMMON  /S_0C<i/  XXI, YYI ,7ZI, AMASS, CG,WG,P9I4S, XX 
COMMON  / oer/  o,rl,cmao,cmcg,qq,s,  vel,thr,cmaoy,cmcgy,ry 

COMMON  /THRST/  TIMES<21> >THSTLB<21» 

c 

DIMENSION  NA  ( 1 ) , XA ( 1 ) 

DIMENSION  S T t 3 ) , WT(  31 , XX  IT (3)  , ZZI T ( 3)  ,C3T(3> 

C 

OATA  CGT  f 55.0,  54.1,  52.0  / 

DATA  ST  / 0.,  .60,  4,5/ 

OATA  NT  / 500.,  450.  , 400.  / 

OATA  XXIT  / 3*2.5  / 

OATA  ZZIT  / 57.0,  54.0,  58.0  / 

C 

* NA (1)  =3 
XA(1)=T 
C 

C**** ***  CENTER  OF  GRAVITY  Or  THE  MISSILE 
C 

CALL  N0TLJ(2»NA,ST»C3T,XA»CG»I£,0) 

C 

C****»‘*  HEIGHT  0r  THE  MISSIVE 
C 

CALL  N0TLi(2,NA,ST,HT,<A,HG,IE,0) 

C 

C*******  MOMENT  0*  INERTIA  ABOUT  THE  LONGITUOINAL  AXIS 
C 

CALL  NOTLJI2, NA, ST, X<IT,XA, XXI, IE, 0> 

C 

C*******  .MOMENTS  OF  INERTIA  ABOUT  THE  TRANSVERSE  AXES 
C 

CALL  NDTLJ(2,NA,ST,ZZir,XA,ZZI,I£,0) 

YYI*ZZI 


j 


XX  S (771  . XXI)**niA3/YYI 
C . 

C”****»  HASS  0?  MISSILE 
C 

AMASS  = wr,/37.2 
C 

£*••••*•  COM* JTAf  ION  OF  ME  AXIAL  F0RC2  DE-IVEREO  3Y  TME  ENJINE  AS 
C A FJMCriOM  OF  TIME 

C 

NA(1)=?1 

XA(l>=r 

CALL  NDrLU(2,MA,TIMES.rHSTLO,XA,MR,tE,0> 

C 

RETURN 

C 

£N0 

SUBROUTINE  ®IX(N,M2,M3i 
C 

c*  • * 

C*  THIS  ROUTIN'  DRAWS  A CALCOMP  PLOT  0s  HE  MISSlE  AND  AIRCRAFT 


C*  TRAJECTORIES  AMO  PLOTS  ®; RTINENT  PERFORMANCE  PARAMETERS.  * 


C 

COMMON  /PL3ATA/  PXET(  52>,PYET<  52),’XE(  52),»YE<  >2),PXEA<  52), 
1 PYEA ( 52),PZE(  5 2) »PZEA  ( 52), F 

COMMON  /®„DAT2/  PAH5(  52),PAVG(  52),»T<  52),PDF(  52), 

. 1 PR<  52) , PAL ( 52) , PBE ( 52),PVDT(  52),PDY(  52>,PDZ<  52) 

COMMON  /P.DAT3/  PAVS(  32),P8ES(  32I,»TS<  32),  »R3 ( 32)  , IPRO 
DIMENSION  ROU),ZOU) 

C 

ORAM  HORIZONTAL  TRAJECTORY  PLOT 
C 

IF(PZEMl)  .GT.5000.)  CD  TO  5 
XM  * *32005, 

XO  « 0000. 

YN  * 32000. 

YO  * -5000. 

CO  TO  15 

5 XF(PZEA(1).CT. 20000.)  CO  TO  10 
. XM  « -50000. 

XO  « 20000. 
yn  > aoono. 

YO  ■ -20000. 

CO  TO  15 

10  XN  = -120000. 

XO  * 30000. 

YH  * 120003. 

YO  * -30000, 

15  °XET (NH)  *<N 
®XET (N*2) -XD 
®YET<NH)sYM 
PYET (N*2> *YD 
PXEA(NH)=XM 

®xea(n«-2)=xd  :: 

PYEA  (NH)  = YN 
PYEA(N*2)=YD 


PXE  (‘J*l)=XM 
PXE  (*(♦?)  -O 
. pve(*m-i)=ym 
PYE  (!!♦?)  = YD 
CALL  P.0T3 ( 30 ) 

•CALL  P'.0T(3.  ,-5.  ,-3) 

CALL  P.0T( 1, 5,2. 0,-3* 

CALL  FACTDR(r) 

CALL  S»AXIS  (0.  ,0.  , 1 3 MC RDSSRANGE  (FT  > , 1 5 , 4.  0, 3 0 . , YM , YD,  - . 65 , 1. 20, 
1 1130, 'N*  , 1, ,-l,.0305> 

CALL  S»AXI>(0.  ,4.0,  1 ♦HDD  WNRANGE  (FT)  ,14,7.0,0.  0 , X M , XO  , 2*60,  *4  ,50, 
1 1190,0. ,1. , -l, .0305* 

CALL  LT  ME  (3XEA,PYEA,'J, 1,0,0) 

CALL  LT ME( 3XE » PYE, M » l * 0 , 0) 

CALL  LIHE(?XET,PYET,M, 1,0,0) 

CALL  SYMDD.(4. 00,4. 03,  126,2, 0.,-l) 

CALL  P.DT(-t  ,-.8,3) 

CALL  PLOT(3.0,-.8,2> 

CALL  P.OT(3.0,4.8,2> 

CALL  PLOT (- 1 . ,4.8,2) 

CALL  P.OT ( - 1 ,-.8,2) 

C 

C 

c*-******  ORAW  VERTICAL  TRAJECTORY  PLOT 
C 

IFCPZEA (l> . 3T. 5000.)  33  TO  25 
RM  = 3S000. 

RO  = -5030* 

Z H = 0. 

ZO  = 1500. 

CO  TO  35 

23  IF(PZEMl).  3T. 20000, 1 30  TO  30 
RH  = 90000. 

RO  = -15003. 

7M  * 0. 

ZO  = 5000. 

SO  TO  35 

30  RM  = 132003. 

RO  = -22000. 

ZH  * 0. 

ZO  - 10000, 

35  =R<m-l)  = RM 

PR(M«-2)  = RD 


pZE (M*t ) = ZM 
PZE(N+?)  = ZD 
R0(  1)  = °R( 1) 

Z0(1)  * °ZE A (1) 

R0(2)  =0. 

Z0(  2)  =0. 

R0(3>  * RM 
<R0(4)  = R9 
Z0C3)  = ZM 
Z0 (4)  = ZO 

CALL  PLOT (14.*  0. ,-3) 

CALL  SMXISIO.  ,0.,13«_TITUOE  (FT)  , 13,  4.  0 ,90,  , ZM,  Z Df-.  65 , 1.  30 , 
1 1190, 90. ,l.,-l,. 0305) 

CALL  S:>AXI5{0.  , 0.  , 10MR4MCE  (F  T)  ,-10, 5,  0, 0. 0 , RM,  RO,  2. 50  , 5290  , 


1 
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CALL  LT»ir<P9,PZE,M,  1,0,0) 

CALL  0\SH.M(*0,70,2,  l) 

CALL  STMIO. <3. ,0.,. 123,2,0. ,-l> 

CALL  P.DT(-l.  3,3> 

CALL  PLOT(7.0,-.8,2) 

CALL  PIDT(7.0,4. 8,2) 

CALL  PLOT (* 1.  ,4.8,2) 

CALL  PLOT(-l. ,-.8,2) 

C 

0*AW  AHS  AND  Ai/G  i/S  TIME  PLOT 
C 

XFtPzrmi.  ;r.5ooo.)  so  to  45 
tm  = 0. 
to  = in. 
go  to  ;; 

4:  IF(PZEMl) . ST.  20000.  ) 30  TO  50 

TM  = 0. 

TO  = 22. 

GO  TO  55  ' . 

50  TM  = 0. 

TO  * 33. 

35  AM  = -20. 

AO  = 10. 

°t  (n+ 1)  = ri 

PTCM+2)  = TO 
°AHG (N*l ) s AM 
PAHG(N*2)  = A0 
PAVG(NH)  = AM 
°AVGOH2)  = AO 
PTS(N2H)  s TM 
■»TS(M2*2)  * TO 
°AVS(N?+1)  = AM 
PAV3 ( N2*’2)  * A0 
' CALL  PLOT (14  ,o.  ,-3) 

CALL  SPAXIStO. ,0,,15MA::EL  CMO  <G»3), 15, 4.0, 30. , AM, AO, -.60, 1.20, 
1 1190, 30. ,1,,-1,. 0305) 

CALL  S=>AXI3<0. ,0. ,10MTIME  ( SEC) ,-l 0 , 5. 0 , 0 . 0 , M, TO, 2. 50 , -. 5230 , 

1 1190,0.  ,1. ,-l , , 0305) 

CALL  LIMK »T,PAHG,M, 1,0,0) 

CALL  LIM£(?T,P4VG,M, 1,0,0) 

CALL  LINE(3T3,PAVS,M2,l,-3,i) 

. CALL  PLOT(0.  ,2.  ,3) 

CALL  PLOT ( 5 , , 2. , 2) 

CALL  P. OT(5.30,3.94,3) 

CALL  P. 0T  ( 5 . 35 , 3,  34,  2) 

CALL  STM90. (5. 70,3.31,  0 805, 3HAH5 , 0 . , 3) 

CALL  P.DT(5,30,3.31,3) 

CALL  P.3T(>, 35,3. 31,2) 

CALL  STM93i.(5,  20,3.72,  0 805, 3HAVG, 0 • , 3) 

CALL  STM9D. (5.  47,3.  807,,  0700, 1,0.,-1) 

CALL  P.OT(-l  ,-.8,3) 

CALL  P.OT(T.  0,-. 8,2) 

CALL  P.OT( 7 » 0 , 4. 3, 2) 

CALL  P.DT(-1  ,4.8,2) 

CALL  P.OTC-l.  ,-.8,2) 


D*AM  A_>MA  AND  BETA  VS. TIME  PLOT 
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PAL  (MU)  = AM 
PAL (MU)  =40 
®BE  (MU ) = 4M 
POE  < M ♦?  > = AO 

®B23(M?U)  = AH  . 

PB23(N2U>  = AO 
CALL  P.0T(1 A.  . 0. ,-3) 

CALL  S3 4X13(0. ,0.,13(4.,  BE  (0E3)  , 13, 4. 0, 90.  , 4M, 4 9 50 , 1.  30 , 
1 1190,90  ,1. ,-l,.030>) 

CALL  SP4XI5(0.  ,0.,10-fTIME  (SEC) ,-10 , 3. 0 , 0, 0 , T M ,TO , 2.50 , -. 5290, 
1 1190,0. ,1. ,-l,. 03051 
CALL  LrX-{=>r, PAL, 9,1, 0,0) 

CALL  LINE (PT, POE, X, 1,0,0) 

CALL  LrHE(PT3, PRES, N2,t, -3,1) 

CALL  P.OT(0, ,2. ,3) 

CALL  PLOH3.  ,2.  ,2) 

CALL  P.DT(5, 30,3.94, 3» 

CALL  PLOT (5, 55, 3. 94, 2) 

CALL  SYMRO. (3.70,3.90,  0805 ,5HAL®H4 , 0, ,5) 

CALL  PLOT (5. 30, 3. 81, 3) 

CALL  P.OTC5. 55,3.31, 2) 

CALL  SYMRD. (3. 70,3.77,  0 605 ,4HBET 4 , 0. , 4) 

CALL  SYMR3.(5.  47,3.307,.  0700,1,0. ,-l) 
call  plot ( - i . ,-.8,3) 

CALL  PLOT(7.0,-.8,2) 

CALL  PLOT (7. 0, 4. 8,2) 

CALL  PLOT ( - 1 . ,4.8,2) 

CALL  PLOT ( -1  ,-.8,2) 


C*******  ORAM  V-ODT  VS,  THE  PLOT 
C 

VM  = -100, 

VO  = 50. 

tF(PVOT(4) . 32. 0.)  PV ) T ( 4 ) = P VOT( 5) 

IF(PVOT(3) , GE. 0.)  3V)r(3)  = P VDT( 41 
IF(PVOT(2).3E.O.)  PV3T ( 2)  = PVOT(3) 

PVOT(l)  = 3VOT (2) 

PVOT(NU)  i VH 
°VOT (N»2)  = VO 
CALL  PL OT ( 1 4 . , 0.  ,-3) 

CALL  S3 4X13(0. ,0. ,15(V*OOT  (FT/SEC2) , 1 5, 4.  0 , 90 . ,VM, VO, -. 60, 1. 20, 
1 1190, 90. ,l.,-l,. 0305) 

CALL  S®  4X15(0.  ,Q.,10-(T(ME  (SEC  ) ,-10 , 3,  0 , 0. 0 , TM,  TO  , 2.50  , -.  5290, 

1 1190,0. ,1. 0305* 

CALL  LIN2(3r,®V0T,M, 1,9,0) 

CALL  PLOT(0, ,2. ,3) 

CALL  PLOT (5 , , 2. ,21 
CALL  PLOT ( - 1 ,-.8,3) 

CALL  PLOT( 7. 0 , 8, 2) 

CALL  PLOT( 7 . 0 , 4. 8, 2) 

CALL  PL3T( *1. ,4.8,2) 

CALL  PLOT(-l. ,-,8,2) 

C 

C**»***»  09AM  OELF  VS  TIME  P_OT 
C 

OH  a 0. 

00  « 4.  ' 

POF(N)  = 0. 


I*1 


* 


P 

P 


h 


"OF  (11*1)  = 3H 
"OF(M*»)  = DO 
CALL  PLOT ( 1 4.  .0. ,-3) 

CALL  SP»XI5(0. ,0.,1S)  OElF  (HZ)  ,15 ,4. 0 , 90 . , 3H ,3D» - . 60 , l. 20, 

I 1130,90. ,1.,-1,.080>) 

CALL  S’AXI3<0.  ,0.  ,10-UTMC  ( SEC)  ,-10  , 5.  0 , 0 , 0 , T 1 , TO , 2.  50  , 529  0 , 

1 1190,0. ,1. ,-l,. 08051 
CALL  Lrn:<*r,POF,H,t,OfO) 

CALL  P.OT(0, ,2. ,31 
CALL  P-OT(i.  ,2.  ,21 
CALL  PuOT(-l  ,-.8,3) 

CALL  P.OT( F. 0 , -«  8, 2) 

CALL  PLOT ( F . 0 , 4. 8,2) 

CALL  PLOTI-1.  ,4.8,2) 

CALL  PLOT(-l. ,-.8,2) 


C*******  ORAd  OY  ASQ  OZ  VS  RiNGE  PLOT 
C 

OH  * -50. 

00  = 25. 

RM  = 18000. 

RO  » -3000.  v 

«*OY(H3H>  » OH 

°OY(M3f2)  * DO 

°DZ(H3*1)  * DH 

POZOI3»2)  * DO 

PR31N3H)  * W 

PR3(N3f2)  » RO 

CALL  P.3T(14.  ,0. , - 3) 

CALL  S9AXI S ( 0. ,0.,Lld3Y,  OZ  ( FT), It , 4. 0 , 90. , 31,03, - . 50 , 1. 35, 

1 1190  , 90. , 1.,-1,.  0305) 

CALL  S»*XI5(0. ,0.,101*£MGE  fFT)»-10,5»0»0.0,Ri,R3,2.53*-«5290» 
1'  1130,0.  ,l.,-l,  , 0805) 

CALL  LISE(9R3,P0Y,S3, 1,0,0) 

CALL  LIME(9R3*P0Z,S3,1, 0,0) 

CALL  P.OT C 0 . ,2 . , 3) 

CALL  PLOT<5.  ,2.  ,2) 

CALL  P.3T1  5,  30,3. 94,  3) 

CALL  P.OT(5.S5,3.94,2) 

CALL  SYH33„<5. 70,3.93,  0805,2HOY,0. ,2) 

CALL  P. 3T(  5 , 30 , 3.  81,  3) 

CALL  P.OT(5,55,3.81,2) 

call  SYSOD. <5. 70,3.77,  0805,2HOZ,0.  ,2) 

CALL  SYSR3. <3. 47, 3. 837, .0700,  1,0., -11 
CALL  P.OT(-t.  ,-.8,3) 

CALL  P„OTCF.O,-.8,2) 

CALL  P. OT (F . 0 , 4. 8, 2) 

CALL  PiOTt-l  ,4. 8,2) 

CALL  PLOTl-t. ,-.8,2) 

CALL  P-OTE ( IOUMMYI 
C 

RETURN 


I o 
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BROOCH  0?TRAJ(lNPJT,  OUTPUT,  PLOT) 


OPTIMIZATION  DF  RETRAN  MISSILE  TRAJECTORY 


DEFINITION  OF  INPUT  VARIABLES 

X0,Y0,7n  - INITIAL  MISSILE  POSITION  IN  EARTH-FIXED 
COORDI NATES  (rT) 

VMO  - IMITIAL  MISSILE  TRJE  AIRSPEED  (FT/SEO) 

GAMO»°SIO  - INITIAL  MISSILE  FLIGHT  PATH  ANGLE  AND  HEAOING 
( RADIANS) 

AL3,9E0  - INITIAL  MISSILE  ANGLE  OF  ATTACK  ANO  SIDESLIP 
ANSLE  ( RADIANS) 

TO  - INITIAL  TIME  (SEO 

HO  - INTEGRATION  TIME  STEP  SIZE  (SEC) 

THAX  - MAXTMJM  FLIGHT  TIME  (SEC)  (SUBROUTINE  STATVEC 
TERMINATES  WHEN  TM AX  IS  EXCEEDED) 

XAO,YAO,?AO  - INITIA.  AIRCRAFT  POSITION  IN  EARTH-FIXED 
COORDINATES  (FT) 

V AO  - AIRCRAFT  TRUE  AIRSPEED  (FT/SEC) 

MAXRUN  - MAXIMJM  NUMBER  Or  ITERATIONS  OF  MAIN  »RD3RAM 
(PROGRAM  TERMINATES  IF  MIN  PAYOFF  NOT  REACHED 
IN  MAXRUN  ITERATIONS) 

TOL  - PROGRAM  TERMINATES  IF  J(N)-J(N-1)  IS  LESS  THAN  TOL 

CO  - INITIA.  VALJE  0s  Control  V ARIA9.E  A0JJ3TMENT  COEFFICIENT 

OELCO  - S T;*  SIZE  USED  IN  SEARCH  ROUTINE  TO  A0JJ3T  C 

KO  10» K020  - COEFFICIENTS  OF  2ND  ORDER  CJRVE  FIT  TO  VELOCITY 
ON  LOS  FJNCTTON 

KlOtKZOf <30  - PENALTY  TJNCTION  COEFFICIENTS 
TR  - ACCELERATION  COMMAND  REFERENCE  RANGE  (FT) 

GGdXAS  - 3IAS  COMMAND  IN  VERTICAL  PLANE  IN  3*3 

CT  - GAIN  COEFFICIENT  f0R  COMHAND  TO  HEAOING  GUIDANCE  ALGORITHM 
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3* I - 0AMPIM3  RATIO  IN  DOMMANO  TO  LDS  3JIOA.MC;  ALGORITHM 

IPT  - IMT-r.;?  SNITCH  3 DMT  ROLL  ING  P_DTTIMG.  PLOT  T INS  fAKES 
PL»Di  OMLY  IF  IPT  = 1 


COMMOM/Aei^:/XL 
OOMMOM/AIC/VAO, 
COMMOM'DNTR./AM 
i,TR,G3^iAs,;r,r. 
common/m^imd/vm 
commom'mic/  \/Mn 

CO**MON/pAYDFF/< 

common/ pi'j;d-/< 
C0MM0':/sTAr*/sx 

COMMOM/PLOATA/X 
1MPTSI,MPTS- 
COMMOM/PLDT2/PO 
1AVGIC  5?>,AM3S( 
■50MM0M/?LDM/oT 
1REPIC  5?>  , 4!_=»S  C 
PEAL  J1 , J2»  JT1 , 
PEAL  <110, <020, 


OSr»  Y ,DSr  »ZL  OSF  ,PSL“ » MLOS 
XAO , Y AO , 7 AO 

U13)  , A/  (410)  ,HO,IRJN,TMAX,OMAH(4lO>  .DMA  V (410) 

ZI 

F,3AMr»PSIF,XF  »YF»  Z”*'AL  F,  REF,  TF, *J,  I >LN 
.GAMO.PSI 0,X0,Y0,Z0, A_0 ,9E0,T0 
01,<02,<1 ,<2,<3,RAM3Er 
010, <023  , <10, <20, <30 

Ull)  ,5Y  (410)  ,SZ<419>  ,SVM<410>  , S®'?!  < 41 0 ) , S3  AM (410) 
I(  521, YI(  52) , ZI(  521, XS<  52>,YS(  52),ZS(  32), 

H I ( 32) , POZI ( 52) » PDHS ( 52),PDZS(  32>,AHGI(  52), 
52),  AMDS ( 52). 

It  32> » 3TS(  52 ) , RA  ‘131  ( 52),RAM3S(  52>,ALPI(  32)  , 
52),3£»3(  52),Y0rr(  52),VOTS(  52) 
JT2,<01,<02,<1,<2,<3 
<10, <20,  <10 


******  READ  INITIAL  CDNOITI 0 NS  ANO  TIME  STEP 

READ  11,XO,YO,ZO,\/MO,31MO,PSIO,ALO,DEO,TO,HO,TMAX,XAO,YAO,ZAO,  VAO, 
1NAXRUN,  TOL  ,30,  DEL39 
) FORMAT (8F10.  0/Z F10.  0 H T/3F1 0. 0) 

READ  15, <010, <020, <10, <20, <30  . 

> FORMAT (5E10. 0) 

READ  l3,Ti>,  333IAS,:T,DZI,IPT 
» FORMAT (4F10.Q.I2) 

°PIMT  21 

L FORMAT(ix//3x,“OPTIMAZITATION  OF  RETRAN  MISSILE  TRA JE3TORY“//3X 
l "PENALTY  FJNTIONS  3DEF- ICIENTS") 

PRINT  ?5, <010, <020, < 10, X 20, <30 

FORMA  T { 1 X/3  X,  •»  H<01  , *10.  4, 3X, 4HK02  , E10. 4,3X, 3H<1  , E10, 4,3X, 3HK2  , 
1E10.  4,3X,3H<3  ,E10.  4/10 

******  INITIALIZE  PLOTTIM3  ARRAYS  AND  RUM  COJNTER 


NPTSI  = 
NPTSF  = 
XI(t)  * 
VI(1)  = 
1IC1)  s 
XS(1>  = 
VS<1>  » 
?S(1)  * 
POHIC1) 
POMS { 1 ) 

porim 

POZS(l) 

AHGim 

AHGS(l) 

AVGI(l) 


S3RT ( (X0-X60) **2*(Y0-YA0)**2) 
PDMI(l) 

ZO-ZAO 

ZO-ZAO 

3. 

0, 

0, 
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Avnsd)  = o. 

°TI  (1 ) = TO 
PTSU)  = TO 

RANSKI)  * 3*>RT  <XO*X)4YO*Y04ZO*7Q> 


RANGS(l) 

* RANGI  H> 

ALPI (1) 

= A. 0*57. 2)5 

AL°S  < l ) 

= A’.  0*37.235 

REPIll) 

= 3*0*57.235 

RFPS(l) 

= 3 £0*57.235 

VOTI (1) 

- 0. 

VOTS ( t ) 

= 0. 

IRJN  = 0 
Cl  = 03 
OELCt  * 

a-',  co 

AH(1»  * 

0. 

AV(1)  s 

0. 

C 

(.»««*«  »•  INTE3RATE  STATE  E3  J®  TIOMS  FORWARD 
C 

30  IR!JN  = IRJNU 
PRINT  *3,I3UN 

33  FORMAT (1X//3X.**  ITERATION  NO.  “,I3/tX» 

IFdc'JN.-:a.  1)  PRINT  13 

19  ?ORHAT(3X,“  INITIAL  T Rft  JECT0RY"/1X> 

IFdRJN  1)  PRINT  ?O,TO,VM3,XO,n,?O,H0 

20  FORMATUX/3X.3HTO  , 1 , 3X,  3H\/H0,F3.2,  3X,  3HX0  ,F7, 0 ,3X, 3HY0  »F&,0, 
13X.3MZ0  »r4 . 0 . 3X  »3HH0  ,F3.2). 

CALL  STATY£C(0.  ) 

if(islx. >a.  i>  go  ra  so 

c 

c»****“*  calcjlat*  payof?  ano  penalties 

c 

CALL  PAYOFdJi) 

YT  S ViF-K01*RANGEr**2-K02*RANGEF 

RERR  = 1503  0, -RAN3*F 

*»SE?R  = »3Ir-°SLF 

CAERR  * G4MF-HL03 

PRINT  ?t» 

33  FORHAT  (lX//3X.”PAYDFr  SNO  ERRORS**) 

PRINT  35»Jl»YT»RER?,>31RR»3AERR 

3s  FORNATflX/  3X.3HJ1  313  7.3X.3HVT  , -10.  2, 3X,  5HR*  RR.  ,F10.2,3X, 
lGHPSERR  ,F3.&,3X,51GAERR  .F8.6/1X) 

C 

c* TEST  STOPPING  CONDITIONS 

c 

XFCIR'JN.  E5. 1)  GO  T3  40  • 

IFdRUN.  3t.  HAXRUN)  53  TO  80 
TF<A9S(J1-J?).LE.T3U  30  TD  80 
SO  J2  > Jt 
C 

{•***•*•  CALC JLAT * OMAN  AND  ONAV 
C 

CALL  OHJ 
C 

»-p»0RN  ONE  OINiNSIONAL  SEARCH  TO  NININIZE  J 
C 

PRINT  33 

33  PORNATCtXMX.-ONEOHSMSlONAL  ALPHA  SEARCH") 


i 

I 
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« * Cl 

DELC  = DEL31 
10  = 0 

do  50  r = i,ioo 

41  CALL  *TAT/iC(C> 

IFIISLM.-D.O)  GO  TO  ,5 
»RINT  43,3 

43  FO’MATdX/ JK.3HC  ,E10  4,3X,"N0  SOLUTION") 

IF(S. or. Oil  3D  TO  io 
4*  Cl  = Cl/50. 

0ELC1  * OE. 31/50. 

C s Cl 

OSLO  = DEL  3 1 
TO  = 10*1 

IF(TD  GT.iO)  60  TO  SO 
50  TO  41 

45  CALL  P»rOFf(JTl) 

3RINT  73,3, JT1,C1,DE.31 
IF(I- ED. 1)  3D  TO  43 
IFCJT1.5E. JT2)  GO  TO  4’ 

43  JT?  = JT1 
C * COEL3 
GO  TO  30 

4*  TFd.r.’E.  2)  3D  TO  3J 
31  = Cl/50. 

DEL Cl  = OE. 31/50. 

C = Cl 

OELC  = 0EL31 
53  COMTINJE 
60  C - C-DEL3 

°RINT  ’5,3, JT2,Cl,DE.3l 

75  FORMAT (IX/  3X.3HC  , E10, 4, 3X, 2HJ  , - 15.  4, 3X, 3R3  1 , E 10. 4, 3X«  ST0EL31 
1,E10.  •») 

DO  70  I = i,N 

AH(I>  = AH(  I)  -3*DHH(  II 

AV(I»  = 4/<I>-3*0HAV<I> 

70  CONTINUE 
C 

C*  ******  RETJRN  FDR  NeXT  ITERATION 
C 

GO  TO  ?0 
C 

C*******  PLOT  INITIAL  AND  -INAL  TRAJECTORIES 
C 

80  '’RINT  S3 

65  PORNAT(lX/3<t"DPTINI7rD  TRA  JECTORWIX) 

»RINT  ?0,TO,tfNO,XO, Y0,70,H0 
CALL  STAT/E3m.> 

PRINT  S6,N»TSI,NPTS«f 

85  rORNAT  I IX  / J X , 5HNPT5I  , * 3 , 3X  , eMPTS*  ,13) 

70TIC1)  = /DTI  12) 

V0TS11)  a /DTS12) 

IPlX^T.'O.l)  3 ALL  ’LTRiJ 
IFCIPON.GE. NAXRUN)  5D  T D 90 
IFCIO  3T. 15.)  30  TD  30 
CTO® 

C 

C**‘**»*  IF  JMIN  NOT  RCA3NE0,  PRINT  MESSAGE 
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LRDN  = I’JN-l 
PRINT  100,  .RUN,  Jl 

FORIATUX/JX,"  JMIN  NOT  RCACHEO  ITERATIONS, 

1/1X) 

ST0® 


“trl 2,5// 


C 

c*+»* 

c* 

c* 

c* 

c**  ** 
c 


•NO 

SU3R0UTINE  STATV£2( 3) 

**4444444*4444444444444, 44444444 44 44444444444444444444441 | 

THIS  S’mo-JTINS  INTEGRATES  THE  STATE  VECTOR  EQUATIONS 


c 

C*»»4 


C0H40M/  A - I NC/Xl OSr»  Y . 03  c »2L  OSF , PSLr  » THLOS 
COMMOM/AIO/VAO.XAO,  Y40,7A0 

COMMON' C'lTRL/ AH  (410)  , A ■/  < 4 10  ) , HO  ,1 2 JN,  T HAX,  OMA M ( ^1 0 ) , 0M  A V (4 10  > 
1,TR,GG3IAS,  0T.3ZI 

COMMON'  1CI N 3/ VMF,  3\  1 ”5  IF,  XF,YF, ZF,4LF,BEF,  Tr,N,  I3LM 

OOMMON'MIC'  VH0,GA'I0,??I0,X0,  Y0  ,Z  0,  A.O  , Sr.O , TO 

COMMON /ST4TE/SX( -13)  , S Y < 410  ) , SZ  (41!)  > , S VM  < 410)  , S^S I ( 413  > ,SS4N  (410  ) 
'0-‘MO‘KPlOAT4/XI(  52t,YI(  52>,ZI(  52),XS(  52)  , YS(  52),ZS(  52), 
1NPTST ,N°rSr 

C0HH0H/oL0T 2'PDHK  32),Pnzi(  52),?3HS(  52>,P37S(  52),4HGI(  >2), 
llV'K  52) , 4 MGS  ( 52)  , A (/••?<.  52) 

COf'MON'PLOT3/PTI  ( 32),^TS(  52),RANGI(  52),RANG3(  52>,4LPI(  52), 
15CPK  52) , 4 L°S  ( 5?)  , 3E3S  ( 52),VOTI(  52  ) , VDT3  C 52) 

PEAL  .M4CH,H,_I31,LI3?,-rG3,LlG4,LlG5,.IG:i,Ll37,LI3  3,LtG9, 
1LVI1,L*I?,LVI3,  LVIV, ‘.VIS  , LV  I& , LVI 7 , L VI 8 ,L  VI9,  . VG1 , WG2,  LVS3, 
lLVGVjLVGOjLVSojLVGX,  . V 3 S , L VG9 

***  INITIALIZE  THE  INTEGRATION  STt»  C3JNTER 


C 

C*«  '* 


1=1 

**•  INIrIALIZE  THE  MISSLE  PARAMETERS 

AL  = A.O 
RE  = B*  0 
VM  = VNO 
GAN  s 5AH0 
®SI  = »StO 
X * XO 
Y * VO 
Z « ZO 
T e TO 
G * .525 
IGJIO  * 0 
H = HO 
ISLN  = 0 
IPLOT  = 0* 

IF(C*  NE. T3, ) GO  TO  1) 

. C » 0. 

IPLOT  = 1 


FIN"  TINE  VARIANT  PARAMETERS 


| 


V 

c 

10  tau  = .01j;fcb*»*51». >9 

OHO  a *.  3305>15*rVJ**«i.26 
HASH  = . 120A*V1/Sm<riJ) 

0 a RM*I*  /H*  V' 1/2. 

M = 12.  *4? 

IF(T.  Lr.  <*•  5)  H = l + .39-.476<»*  (T-.6J) 

TH  = 0. 

IF(T.  IT.  4.  5)  TH  = 2000 
C 

C«***«*«  CAi:jLAri  COMMAXDiD  ALPHA  AND  3;TA 
C 

AVV  s AVU)  -3*OHAV< II 
IFCAVV. 57.20. > AVV  = 20. 

IF (AVV.  17,-20.  ) AVV  = -20.  . 

AVV  = AVV* 32 -2 

AHH  = AH  ( I ) -0*  DHAH (I) 

TF (ASH. OT.  20  ) AHH  = 20. 

IF(AHH.LT.-2n.»  AHH  = -20. 

AHH  = AHM*  32.2 

AL  * (»VV  *32.  2*CDS  ( S’-.  H ) ) *H/  ( 20,  1*S*Q) 

IFUL.ST. . 350)  AL  = .370 
IF( AL . L 7 • - . 350)  AL  = - 350 
9E  = -AHH  * \t  ( 20. 1*5*3) 

IF(3-  GT..350)  HE  = .3:0 
IFOZ.LT 350)  = -.350 

ALTOT  a S3RT1 AL**Z*3E*‘2> 

C 

C****,M  CALC JLA  T £ ATRODVNAHIC  FORCES 
C 

COO  a -51,  S53*HACH**7*?90.974*MACH**3-115l.  2«4  4* MA  CH**  >♦  18A’4.  332* 
lHACH“*t-t,!3  731»HACH**3+865.  362*  HACH*  * 2-243.  2 93*H  ACHi-28. 253 
OELCC  a -2,  851£-05*A.TDr**2-1.82d6L-03*ALror*3.371’4t-03 
OELCOF  = -7. 4060c-l2*Z‘7-1.0785E-07*(-Z)-6. 8773E-Q5 
FX  = -5*0*  (C.)0*DELCOO‘LCOF) 

FY  a -S*C*<23.  l*BE) 

FZ  * -2*1*(20. 1*AD 

r 

C*******  JNIiSRAI-i  HISSLi  -aJATIONS  OF  HOTION 

c 

CAL  = :OS(A.) 

SAL  a SI-HAD 
C9£  a ;3S(3i) 

3B£  = 5 1 M ( 3 £ ) 

VMOOT  a ((rX*TH)*CAL*C3E*FY*SBt*FZ*SAL*Cf)E)/i-32.2*SIV(CA'0 

ay  a<-(-x»rn) *cal*»9:*-v*C3E-fz*sa.*3DE>/h 
AZ  =(-(FX*rH)*SAL*rZ*CAL)/H*32.2*:3S(GAH) 

®SIDOT  a A V/ { VH*COS ( GAM ) ) 

GAXOOT  = -AZ/VM 
VH  * VHf H*  VMOOT 
CAM  = *,AH*H*r.AMOOr 

psi  * »3iH*»snor 

xoor  a vH*:D3(r.AM)*:Ds(Psi) 

VOOT  a VM*C3S<r,Art)»SI'HPSI) 

• ZOOT  a - VH*  SIN ( GAH) 

. X S X*H*XDDT 
Y a V+H* YODT 


Z a 7*H*Z3Dr 


I 

i 


I 


■ 


i 


c 

C*«**M*  ST07;  STATE  VA.JES 

c 

TF(C.N*.0.)  30  TO  12 
SVKH  a YN 
3GAHd>  = SAN 
SPSI(T)  = ?5I 
SXd)  = X 
SYd)  a r 

szm  * z 
c 

C*“**“*  T?sr  F3*  STOPPING  33M0IT ION 

c 

12  T = T*N 

IF(T  GT.THAX)  GO  T 0 53 
TFCZ.OT.I.)  SO  TO  ;o 

°SL0S  a ATAN(YA0/XA0> -(  T-T0>*VA0/S0RT(XA0**2l-fl0**2> 

THLOS  a ATAN< ZA0/3Q?r( X A0»*2+YA0**2> » 

RANG;  a 30RT(X»*2fY*»?l 
XLOS  =-R*N3;*SOS(:>SLOS) 

YLOS  a-RANGi’SINdS.OSt 
ZLOS  = RAN3E*TAN(TNL0S) 

OZ  = Z*  ZLOS 

OH  a S0»7(C  X-X’.OS)  **>♦<  Y-YLOS)  **21 
W * X-X'.OS 
OH  a SIVUOH.W) 

OVZ  = /N*SIN(T-IL0S-3*Hi 

OVH  a XOOT»:OS  (PSLDSt-3  1M5/2.  )+YOOT*SIN(P5l03*T.  1 V16/ 2 • > 
l-VA0*OANG;/3GRT  <XA3**2»-YA0**2) 

DERR2  = <.X-XlOS)**2«'(Y-YLOS>**2*(Z-Z-OS)**2 
IF(0'PZ2.Lii2500.)  GO  TO  15 
C 

£».»«»**  MAKE  ACSiLEPATION  30MMAN0  OECISIOM 
C 

T * T + l 

iF(i»Lor.  n.  i>  go  ro  23 
IFdRUN.EQ.l.  ANO.D.tO.O.)  GO  TO  20 
IFdPUN.NE.  1.  ANO.I.L;.  M>  GO  TO  10 
IF(C. Ns.  9t  * AMO. I.LE. M)  SO  TO- 10 
OELAH  a AHC  M) * AH  <M*1 3 
DELAY  a AV(N)-AV<N-1) 

OELOHH  a ONAH(M)-DNAN<N-i> 

OELOHV  = ONAY(N)-DNArf  < *4-  1) 

1H( I)  a AH ( I - 1 ) +Ot. A H 
AV(I)  a AYd-lM-OE.A/ 

•*HAH(I)  a OHAHd-lJ  OE.OHH  * 

JHAYUJ  = OHAV  <I_ll  OELOHV 
GO  TO  10 
C 

r-coro  final  coNoirioNS 

c 

15  TF  a T 
XF  » X 
XF  * Y 
ZF  » Z 
VHP  = YH 
3AMF  r SAN 
*$IF  a PSI 
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\ X- 


ALF  = U 
3EF  s 

if(c  n.n.i  n=i 

XLOSF  = XLD5 
YLOSF  » Y.OS 
SLOSF  r ?lD5 
PSLF  = PSLOS 

IF(0  EM.I  3RINT  ir,  r,  V'1iX,Y,Z,3-».')%AM(I-1I  .«V(r-l> 

17  FO»M4T(1X/3X,2HTF  , - 5. I , 3X, 3HVMF  , *9.  ? , 3X, 3MX-  , F7 , 0, 3X , 3MY-  , 

1-6.0,3X,3MZ-  ,'8.0,3<,4HDHF  , F8.2,  3Xt  «H0ZF  , ?J  . 2,  1 < , 3RAH  ,FS.Z,3X, 
13HAV  ,r5.2» 

C 

RETURN 

C 

£••«««•»  30m3JT;  T3JAmSF0RMUtou  FROM  GIMJA.  AXES  TO  EARTM  AXES 
C 


CPL  = 

cost? 

3L0  5 ) 

SPL  = 3 IN<  3 

3L0S) 

CTL  = 

cos  ( r 

HL0S) 

STL  : 

SIM(THLOS) 

LI51 

= COL* 

CTL 

L1G2 

3 -SP. 

L1G3 

= SPL* 

STL 

LIC4 

= 3PL* 

OIL 

LIG5 

3 C°L 

LIGS 

s S°L* 

STL 

LIS  7 

= -ST. 

LIG4 

= 0. 

LIG9 

= STL 

c 

COMPUTE  TRAfJSFORMUTOM  FROM  EARTH  AXES  TO  VELOCITY  AXES 
C 

ops  = :rtsc3si) 
sps  = 3 1 ’I  ( 3 s i ) 

CGA  s OOS(3AM) 

3GA  * 3IM<;AM> 

lvii  = :ps*ssa 

LVI2  = SOS»25A 
LVI3  = -SSA 
LVI4  » -SP3 
LVI3  * 203 
LVI5  s 0. 
lvi7  = cps*ssa 

LVIS  = SPS*3GA 
LVI9  = CGA 

C . 

C0M3UTE  TRANSFORMATION  FROM  GIM3A.  AXES  TO  VELOCIT 1 AXES 
C 

LVS1  * L/Il»LlGl*L/I?*LlG4»LVI3*LI37 
• LVG2  = LVri*LlC,2+LVI?»LlG5fLVI3*LICa  . 

LVG3  * L'/Il»Ll‘.3+L/I2*-.rG6»LVI3*L  TJ9 
LV3<*  s LYI«».I3l+.vr3*.IG4»LVI6*tr37 
LVG  > * LVI«*LIG?*L/I 5*'.IG5tLVI6*LI3S 
LVSS  * lVI«»LK.3*.VI5*.IG6*LVI6*Li;0 
LVG7  s L7IT*LIGlHm*.ir,4«>LVI9*Ln7 
LVG 9 3 LVI**lT3?*LVH*lIG*»»tV19*Li;a 
LVG9  s L7ir*LlG3«-Lm*'.ir.b*LVI9*LlC9 

n * LVSl*LVS»»LVG3»LV3»»LVG6*LVG7».v;j*LVG4*.V39-./G3*LVG5*.VG7- 


( 


i 

f 

1 
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o o o ooo  ooo 


c 


tLVG2*LV34»_/3:)-LV3l»-V*.5»LVG9 
IF<t»LDT.i5.1)  GO  TO  31 


••••*•“  COMMUTE  SUIOANCi  SO-MANOS  FOR  FIRST  RUN 

CH  = CT*T<*/?ANGE 
Cl  = .05*1MTR/RAn;;»»»2 
C2  * C*I“.  2^9-  CTR/RAN3D 

•*  IF  HEAD  OF  LOS#  OO-MANO  TURN  .TO  PAOAL.EL  . DS 

IFCIGUTO. SO. l>  GO  TO  23 
IFIOU.LT.  0.  ) GO  TO  25 
AH< I)  s CH* (3SLOS-3Sl) 

AXG  = -LV32/0*AH<I) 

AHS  = LVG5/D*A-*(I) 

AVG  * C1»0?*C2«0VZ>333IAS 
AV(I)  = -(uVr,7*AX3«-L»/5i*  AHG-L VG9* AVG) 

AHTSST  = -Ol'OH-CZ’DVH 
IF ( AHTSST. w T . 0 . ) 30  TO  31 
IFtT.GT.?.)  I3UI0  = l 

**  COMMAND  TO  LOS 

25  AHG  = -31*DH-C2*OVH 

30  AV3  = 3l*0Z*02*DVZ*,3»3t  AS 

AH( I)  » LVG5»AHG-L/G>»& VG 
AV( I)  * LV33*AVG-LVG9MHG 

C 

C»***»»»  PRINT  OJTPUT  ANO  U»OATE  PLOT  ARRAf S FOR  FIRST  A NO  rINAL  RUNS 
C 

31  L * 2 0/H0 
K = I/L-(I*t)/l 

PRINT  J«5,T,  VM,X,Y,Z,OM,OZ,AH(I-1),AV<I-1) 

3S  PORMAT(lXZIX,2-tT  , F 5.  1 , 3X,  3H  VM  # r 9. 2t3<t3H<  ,-7.0,3<»3HY 

lp6.  0 , 3< , 3HZ  »F«.0,3<»  HOM  , Ffl.2, 3X,  !»HOZ  , r 9 i 2,  3<#  3R  AH  ,FS. 
13HAV  ,FS.2) 

IF(K  SO. V 30  TO  1J 
IF <I*LDT .*0>1)  GO  TO  V5 
NPTSI  * MPT SI *1 
XKMPTSI)  * X 
YKMPTSI)  » r 
ZKMPTSI)  x -Z 
IFUHS.3T.20  > AHJ*23 • 

IFUHG. Lv.  -20.  ) AH3  = *20» 

AHGI ( NBTSI)  * AHG 
TFUVS.3T.20  ) AV3*21. 

IFUVG.LT.-20.  > AV3*-20. 

AVGI  ( fPTSIl  = AVG 
PO«I(fPTSI»  » OH 
POZHN^TSI)  * OZ 
°TI < NPTSI)  * T 

RAN3KNPT5P  x soRr(<*<»v*r«-z*z> 

ALPI(NPrsi)  = AL*57#  293 
•l£PI<N«*TSn  « Pc’57,23: 

VOTKN’TSII  s vmoot 

▼FCViITT  CPTSI)  . GT.tO).)  VOTUNPTSI)  = 100. 

IF(VOTMM9rSIl  . LT.-100  ) VOTUNPTSI)  = -100. 

\ ’• 
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r\>  • 


I I 


40 

43 


C 

50 

.53 


C 

C 


C 

C*J 

C* 

C* 

C* 

c* 

c*- 

c 


r, o ro  to 

•jptsf  = mptsf* 

X3C*l°T5e)  * X 
YSOJPTSF)  = Y 
ZS(N"T3F)  = -Z 
AMGS(HSTSFI  =( 
IFCAHG3  CIPT  3F) 
IFl  AHG?  (*IPTSF) 
Avgs(m»tsf>  = 
IF|AVG^C1?T3C) 
iFuvGscpr.sn 

POHS(N3T Sr)  = 

pozs(m»tsf>  = 

PTSCMPTSF)  r T 
RANGS(N°T3r)  s 

ALPS ( N*T3F)  = 
9EPS(Nns->  = 
VOT  3 ( NpTSF)  = 
TF(YOT«:{fl:»rSc) 
IF(Y0T>  {’l3r3F) 
GO  TO  10 


LVG9*AM(  r 
. GT. 20.  > 

. LT. -20. > 
(LVGJMV  ' 
. GT, 20.  ) 

. IT. -20.  i 

OH 

OZ 


♦L7G6*AV(I)  \n 
- jS  (fJPTSP)  =20, 
l ^GSOJPTSF) =-20* 
, ♦ ;.7G5*AV(I))/0 
«'G3  (NPTS")  =20, 
,A75S(NPTSF>=-20. 


sort  t x* / * 

AL*5Z, 29' 
3£*5Z. 29' 
VMOOT 
. GT, 10  3.) 

. LT. *13  0 i 


'*Y*Z*Z> 


/CTS  ( NPTSF)  = 
OTS(NPTSF)  = 


100. 

>100. 


IF(C  EO. 0.1  »RINT  55 
FORMAT!  1X/3X, "NO  SDL  JT T G *»*“> 
ISLN  = 1 


RETURN 


ENO 

SU3R0UTINE  DHU 


THIS  SUBROUTINE  I NTS'  >i“‘-S  THE  HMDA-03T  AMD  X-ODT  EQUATIONS 
9AC<WARDS  AND  STDREr  '*:u<D  ^ DHAV(I) 


COMMON/ AFI.MC/XLOSp,Y_Or 
COHiON/CMTRL/AH (413) , A / 

1,TR,GG3IA3,CT,ZI 
COMMON/MIC/  YMO.GAMO,®' 

COMMON/ MP I MC/VMF,51Mr»:* 

00NM0M'oAYD‘r:r/K0l,<0  2,< 

COMMON/ST  AT  E/SX  (*.10  ) , S* 

REAL  N104,M,LVMF,L°SS,-  --  

1<02,  <li  <2,<3,L0TV,.0T*,L  j*LDTX»LDTY,LOTZ 


,r.OSFfP5LF,rHLOS 

-I  I) ,HO,IRJN,TMAX,OHAM(4lO)  , DM A V (4 10) 


* • ' C.YO.ZO, ALO.OiO.TO 
I - . XP,YF,ZF,ALF,3EF,r-,M,I3LN 
, * T »i<3,PA,IJEf  ,/ 

--E)  ,SZ  0.10)  ,SVM(  410)  ,S°3I(  410)  , S3  AM  (410) 
- - • LXF.LYF.LZT.LYM.L’S.LSA.LX.LY.LZ.XOl, 


c 

£••«»»«•  CALC JLAT E LAMOA(Tr) 
C 


RANGE F = S3RI(XF*XJt^YF*  /F’-EF*ZF) 
IF(RAN»EF, ED  0.)  R1N3E?  * 


LWMF  s -1, 

LPSF  = <1*(PSIF-PSLF) 

LGAF  = <2*( 3AMF-TM.03) 

LXr  * *. •<31»Xe>K0,*<F/''1  *3Ec-K3*XF/RANGEr*(13003. -RAMGEF) 
LYF  = ».  *<0l*YF»K02*YF/>i  ‘ 3Er-K3*Y?/RANGEF* (15000. -RAMGEF) 
LZF  ■ ?.  •Ol*RFfK0?*  E*/**-’  ♦•-P“K3*ZF/RANG£F*  (13000. -RAMGEF) 
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r 


c 

£•»«*•«»  INITIALIZE  PARAHETERS 

C 

AL  * ALP 
IE  = REP 
VM  = VHP 
CAM  = SAMP 
PSI  = °5IP 
X * XF 
V = VF 
! * ZF 
T = TF  * 

S = . 753 
H * -HO 
LVM  = LVHF 
LPS  = '.°SP 
LGA  * LCAP 
LX  * L<r 
LY  = LVF 
LZ  * LZF 
I « M*1 
PRINT  24 

24  P0RMAT(1X/3X,"RAC<HAR3  IMTEGRATION'VIX) 

C 

cal:jlate  time  variant  parameters  and  derivatives 
c 

10  I * 1-1 

TA'J  = . 09356:,*  24515,  59 
RHO  * 5. 3303i-15*TAa»*:. 26 
HACH  = . 0234*VM/S3RT(TAU> 

R s RHO* VM*  VM/  2 • 

M * 12.  42 

IF(T.LT.4.5)  M = 14. 39-. 4764* (T-. 60) 

TH  * 0. 

IF(T.  LT  • 5 ) TH  * 2000- 

OTAUZ  = -033366 

0RHO7  * 2.71HIE-14*TAJ**3.26*0TA'JZ 
IMA CV  » . 020  4/SQRT (TAJ) 

OMACZ  * 0102*VM/TAJ*»1.5*OTAUZ 

RQV  * RHO* VM 
OQZ  * VM**2/2. *DRH1Z 
C 

c*******  cal:jlate  commanded  alpha  and  seta  and  derivatives 
c 

AW  = AV  (II 

IF( AVV, GT, 2 0- 1 AVV  = 21. 

IF(AVV.LT.-20.  > AVV  * -20, 

AVV  * AVV* 32. 2 
AHH  = AH  (II 

TF(AHH. 5*. 20. > AHH  *20. 

IF (AHH, LT.-20.)  AHH  * -20. 

AHH  * AHH*32  2 

AL  » (»VV  ♦■3?.2*C1S(GiMn*M/(20.1*S*ai 
IF(AL.5T.. 350)  AL  = ,350 
IF(AL..T.-. 350)  AL  = - 350 
HE  ■ -A MM  *M/(  20. 1*5*  Q ) 

IF(RE.3T..350)  OE  * .350 
IF(SF..T.-. 350)  US  = • 350 
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ALTOT  x SaU(Al**2*GE**2) 
ifcaltdt  £i.o. i Airor  » i.e-io 
DALV  = -AlO’OOV 

OALG  x -32. 2*97(20. 1*5'D)*SIN (GA9) 

oal7  * -u/o*oaz 
docv  = -3E/a*oov 
ooez  = -G£/D*OOZ 

OALTV  = AL741T0TO1LV*  ? t /ALTOT  *03  £/ 

OALTG  = AL/4LTOT*01LO 

0ALT7  * AL'ALTOT’OXLZ^E/ALTO^OGEZ- 

OAVVA  x 32. 2 

OAHHA  * 32.2 

0AL4V  * 9/(20. 1*S'0)*0?VVA 
DBEAH  x -9/(20. l*S*a« *3AHHA 
OALTAH  x 9E*57, 295'A_TDT*D9E4H 
GALTAV  = 1. *57. 295/4. TDT*OALAV 
C 

C»*»**»*  CALO'JLAf  £ 4ER0DT N49!  0 FORCES  AMD  DERIVATIVES 
C 

COO  = -51,  S53*rtAC9**'*?B0.974*MAOH**3-1151.244*9409**S*1844.  332* 
1 MACH* *4- 135 9 781*910 9* »?*835.  38 2*9 AC H* *2-24&,  2 3 3* 9 1C9* 28, 253 
OELCC  x -2.  3sl£-05*A.T-'T**2-1.8286i-13*ALTOTH.571'4E-03 ' 

DELC'DF  = .n30£-l2*Z* Z-i. 0785E- 0 7*  ( - Z) -3.  8779E-05 

FX  = -S*a*(ODO*OELOOO-  LCOF) 

FY  = -9»G*(20. 1*BE» 

FZ  x -S*0*(20. 1*AL» 

OCOOV  * (-331. S06MA0H**  6*  2235.  8-»4*M10H*  *3-5755. 223*910  H**4> 
f*37r.  325*9409**3-5039.  ?43*94CH**2*177l. 754*9A09-2'»3.293  ) *0913V 
OCDOZ  x (-331.  306*9409**6+2235.  844  *9409**5-575  3. 22  0 *9109**  41- 
17377.  323*9109**3-5033. ? 43*9 ACH**2*t 7*1 , 73V* 910 9-2 4 5* 2 9 8) *0915Z 
OOSLCV  x (-5  702E-05»A'.TOT-1.  3286E-03)  *OALTV 
OOtlOS  x (-5  702E-05*A.TOT-1.  82e3E-03> *01LrG 
DOELCZ  x C-5.  702£-J5*A-.TOT-l.d266£-03)*OALTZ 
3DELPZ  = -14.  3120£-l2*Z*1.0785c-07 
OFXV  x -cx/D‘DOV-S»1*(jOOOY*OOELOiO 
OFXG  = -5*D*mELC0 

OFXZ  x -FX/3*OOZ-S»3»(DOOOZ*OOELCZ*ODELFZ) 

OFYV  = -S*D*20. 1*D3E V-3 * 20. 1* BE*DQV 
DFZV  » -S*D*20.  1*D1LV-5*20.1*AL*OQV 
3FZ5  x -S*D*20.  1*0U; 

30ELAM  x (-5  702E- 05*4.7 OT- 1.  9286£-03>*01LTA9 
ODEIAV  x (-5.  702E- 05*4. TOT-1.  8286E- 03)  *0AlT4V 
3FXAH  x -5*D*D9EL49 
OFXAV  x -S*  D*  'DDELAV 
OFY4H  x -S*D*20. 1*D3:A9 
DFZAV  x -3*3*20.  i*D4;.4V 
C 

C**  *■♦**  CALO-JUATE  V9DOT,f»SI*OT  ANO  GA90DT  AND  DERIVATIVES 
C 

CAL  = OOS(AL) 

SAL  * «IMUL> 


COE  x OOSOE) 

SBE  * *IN(3E> 

VH30T  » ((r<*rM)*0AL*03E*FY*SBE*?Z»SAL*CDE>79-32.2*SI9(GA9» 
AY  x|-(rx»rH) * CAL* SB- f - Y*C3E-FZ*S4l*$9E> /9 
A7  x(-(?x»r9l *SAL*rZ*0»L)/H»32.2*0DS(GAM) 

PSI01T  x 4r/(VM*CD3  (349M 
GA900T  X -AZ7VH 


v 


uuu 


0V3TV  = (3- XV*CAL*39i-rX*SAl.*C0E*0Al/“FX,'CAL*:>9r-*O3SV*-0FYV*39E*rV* 
lCnE’OP*  VOr?V*SAL*:3E*'Z*CAl*CaC*3ALV-FZ*SAl*5  3t*33tV-TH*SA.*305* 
nALV-T-l*04.*>«V:*99£  VI  /f 

OVOTG  5 (0?XG*CAL»:3:*rX*SAL*CPE*0At.3*DF73*SAt*C9E*FZ*CAL*33E*0AL3 
1-TM*SA;.*C9£ *3ALG)/-I-3Z  ?*C3S<GAM) 

9V9TZ  * <3rX7*CAl*:9:-  - X*SAL*C3E*34LZ-FX*CA_*3  3E»  3 3EZ*eY*CB£»OOE7* 
1FZ*CAL*31£*3AL7-FZ»SU'  SOE*O  JEZ-TH*SAl.*C3i»DA.7-TH*r.A.*SnL*30EZ)/1 
0P3TV  * ( - < 3rX V • •:A_*S9E-(FX»H)  *5Al*53£*JA.V-(*X*M) *CA.*39E* 
10BE7O- VV*:3>FX*S3£*33EV-DFZV*SAL»S3E-FZ*CAL*59E*3ALV-FZ*34L*C3E* 
n9.-V>/<M‘:3S<GAM)*V1)  -^SIOOT/VM 

OP3TG  * t-O-XG  \ *;*.‘S9£-(FX*T-I>  *34L*S3E*3At-r.-3  = ZS*S  AL*S3E-FZ 
i*CAL*S3"*3AL0)  /<M*:03tr,4M)*VM)  +SIN(64M)*PSI00r/C0S{SA'O 
npor?  s <-(3-XZ  ) *;A_*S3£MFX+r-l)  *?AL*33"*34L7-(rXfM)  *CA.*CBE* 

108EZ-Fr»S3£*33-:z-F7*3A.*SO£*OALZ-Fr*SAL*COE*D3*T)  / ( M»  30 S (G A M)  * VI ) 
0G9TV  * - ( - ( 3FX  V ) »Sit-(FX*TH)*:4L*0ALY«-DFZY*36L-FZ*SAL*0UV)/ 

1 <M*V'M-‘34M3DT/VM 

OGOTG  = -(-(3FXG  ) * SAL- (FX+TH) *2 AL*D ALS+0FZ3* CAL-FZ*S Al» 04LG) / 

1 (H»VM)*3?.J*SIH(GM)  /VM 

0G3TZ  = -(-(OFXZ  ) »S1L-(FX  + TH)*C4L*9ALZ-FZ*'3M»3U7) /<M*V1) 
9V0TAH  = ( ( 3rX  AH  ) • :4L*C3E-  (FXHH)  *CAL*S9£*  3BE»rtfOFYAH*SBE*FY* 
1CBE*09£A4--Z*SAL*33£*03£AH)  /M 

OVOTAV  = ((3-XAV  ) *C5L*C9E-<FX«-n>  *SAL*CBE*94LAV*-0rZ4V*SA.*CBE* 

1FZ*CAL»C3E*34LAV)/M 

OPOTAH  = (-(1FXAH  ) *2AL*SBE-(FX«-rH)  *CAL*CBE*93EA-t«-DrYAH*C3E-FY* 
iSB£*Qfi*AH-r Z*SAL*23£*3SEAM)  /(  M*V1*30S  (GAM)  1 
0P-3T4V  - (•  ( 3FXAY  >*:AL*S9E+(FX*rH)*SA:.*S3E*34L4V-0rZAV*S4;.*S9E- 
1F7*CAL*53£»3ALAV)/('1*Y‘'*C0S  (GAM)) 

0G3TAH  = (3r<AH  )*3AL/(M*YM) 

3G0TAV  = -(-(9FXAY  ) ♦ S At-  (FX  + TH)  *C4L*D«lAVlOr7AV»CA.-FZ*S4L* 

10ALAV)/t‘1»YM) 

C 

C*»«**f  UPGUE  Yi.PSI  AMO  GAM  ANO  CAL3U. ATE. DERIVATIVES 
C 

VM  = SVM(I-l) 

PSI  = SPSKI-l) 

GAM  = SSAM(I-l) 

CGA  = 30S( 3AM) 

5G4  = SI’JOAM) 

CPS  = :0S(’SI) 

SPS  = 5IM(95I) 

XOOT  = V t*CG4*CPS 
YOOT  » VM* 2 34*  SPS 
*?00T  * -VVS3A 
9X0TV  * CG4*3PS 
3X0TP  = -VM*2GA»S?3 
QXOTG  » -VM*5GA*CPS 
OYOTV  » CG»*SPS 
OYOTP  * VM*C3A*CPS 
OVOTG  * -V1*3GA*S93 
OZOTV  * -S3 A 
OZOTG  * -VM*2GA 

**  U°OAr:  X,Y  AND  l 

X « SXCI-1) 

Y « SY(I-l) 

Z « SZ<!-tl 


( 
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CALCULATE  LAMDAODT 
C 

LOTV  = -(HM»0V0mL»3,DP0TV*LGA*0;0TY*L<*3OrY»LY*DYm*lZ*DZDTY> 
LOTP  s - UX*OXOTP».Y*  DYOTP) 

LOTG  * -Ci./'l’iJYDTDH’S*  DPOTG*LGA*DOOTS*LX*OOT3fLY*9YDTG*LZ* 
1DZOTG1 

LOTX  =0.  . 

LDTY  =0.  . 

LOT?  s -(LrfM»OVOTZH»5’OPOTZ*LGA*OGOTZ> 

C 

C»»***»»  INTEGRATE  LAHDA-ODT  TO  GET  LAMDA. 

C 

LVM  = L YiH’LOTV 
LPS  = .OSM’LOTP 
LGA  = .34M*LOTG 
LX  = LXm’.DTX 
LY  = LYfH*wDTY 
LZ  = L?fH*.DTZ  ‘ 

C 

&>*•••*•  CALCULATE  ANO  STD?*  OHAH  AND  DHA V 
C 

OHAH  ( II  = !.\M*OVOn*H_:»S*DPOTAH4-LG4*DGDT4H 
OHAV(I)  = . Yi*OVOrAYH.a3*OPOTAV*LGA*DGOTAY 
T * T + H 
L = 2. /HO 
K - I/L*(IHI/L 
IFTK  ED. 0)  30  TO  31 

PRINT  ?5,T,  </M,  OHAH(I) , DHAY(I)  , X ,Y,  Z , AH  ( I)  , A V(D 
?•?  FOPH4T(lX/3X,2HT  5.1,3X,3HVM  , F3. ? , 3X, 5HDHA H , E 15.  7 , 3X  , 5H0HAV, 
1"15.7,3<,2-K  ,F7.0,3X,2HY  ,F7.0,3X,ZAZ  ,c7. 0, 3<,3H4rt  ,F6.2,3X, 
13HAV  ,-3.2» 

C 

C»»»V”  CHED<  ?DR  S TOP^I N3  CONDITION 
C 

30  IF(I  GT . 1)  OD  TO  11 
. RETURN 
EN3 

SUBROUTINE  PAYOFFOl 
C 

C*  ' • 

C*  THIS  SUBROUTINE  lOH^JTES  THE  VALJE  OF  THE  PA YOcF  FUNCTION  J • 

C*  * 

c 

oomhon/afihd/xlosp,y.def,zlosf,pslf,hlos 

CO-1HON'*1FlND/VNF,3AMFf  "S  IF,  XF  , YF,  ZF,  SI  F,  9EF,  T"  ,N,  I SLN 
OOMHON/°AY5-:r/K01,KO?,<l,K2,K3,RA'i;£- 
COMHON/P-:h3DF/K010,<021,  K10,K20,K30 
real  J,<01.<02,K1,<2,<Y 
REAL  K011,<0.20,K10,<20,K30 
C 

C»****«*  ESTABLISH  PENALTY  FJMCTION  COEFFICIENTS 
C 

<01  3 <010 

K02  * <020 
<1  > K10 
K?  m K20 
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K3  = <30 


. 1 


C 

MODIFY  PENALTY  r U N 0 T I ON  COEFFICIENTS  IF  FINAL  VALUES  ARE 
C WITHIN  r O'.ERANCES 

c 

• RANGE F s ST»T(XF*X-»YF' YF*ZF*ZF) 

IF<s>?:  = -°5.  = ,LE..?iO.A-n.PSIF-PS--.SE.O.)  <1  = 0. 
IF(A03(3AMF-THL0S  > . .£  .175)  <2=0. 

IFUAN3EF.SE  16000. > <3  = 0. 

c 

C*******  CALCULATE  PAYOFF  =UNCTION 

C 

J = -VMF+Ol  *RAN3EF**?f  <02  *RAMGE-+<1  /2.MPSIF-PSLF)  **2+<2  /2.* 
KGAMF-TMLOS)  **2  + <3  / 2.  • ( 160  00  .-RAN3Er)  **2 
C 

RETURN 

ZNO 

SUBROUTINE  PL T RAJ 
C 

C* 

C*  TMIS  5!I3 ROUTINE  *L0T>  COMPARISONS  3ETWIEN  TME  INITIAL 

C*  PERFORMANCE  PARAMETERS  ANO  THE  OPTIMIZED  PERFORMANCE 

C*  PARAMETERS 

C* 

C 

COMMON/ AFIMC/XLOSe,  Y. OS F , ZLOSF, PSLF, THLOS 
COMMON/ A 10/ V A0 ,XA0, Y 40, ZAO 

COMMON/PL3AT A/XI { 52),YI(  52)  ,ZI(  52>,XS(  52>,Y<5(  S2),ZS(  52), 
INI, NF 

COMMON/PLDT 2/®OHI ( 52),POZI(  52),P0HS(  52>,P0ZS(  52),AHGI(  52), 
1AVGK  5?>,AH3S(  52)  , A VGS  ( 52) 

COMMON/PLOT 3/PTI ( 52),  = TS(  52),RAMGI(  52>,RAN3S(  5>),ALPI(  52), 
19EPIC  52>,1_°S(  52),3E:>S(  52),VDTI(  52),V0TS<  52) 

OIMENSIOM  R3U)  ,Z0(L) 

C 

C**»»»»*  DRAW  HORIZONTAL  TRAJECTORY  PLOT 
C 

IF (ZAO .ME. -35000.)  GO  70  5 
XH  a -50000. 

XO  = 15000. 

YH  = 60000. 

YO  * -15000. 

GO  TO  15 

5 IFtZAO. ME. -20000.)  GO  TO  10 
XH  = -50000. 

XO  a 10000. 

YH  a 50000. 

YO  a -10000. 

GO  TO  15 
10  CONTINJE 
15  XI(NIM)  = <N 
XI  (MI  •*•?)  = O 
• XS(MFM)  I XM 
XS(NF ♦?)  = XO 
Vt(MUl)  s YH 

VI (MI ♦?)  = Y0 
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YS(nr*n  * ym 
v soif*?)  * ro 

CALL  PLOTSOO) 

CALL  P.OT(5, ,-5.,-3) 

CALL  P.OT(l.5,2.0,-3l 

CALL  SPAXtSCO.  ,0.,15H:’3SSRANGE  (FT) ,15,4.0,90. ,Y*S  YD, -.65,1.20, 
1 1150,90. ,1. ,-l,.030>) 

CALL  S9AXI  > ( 0.  ,4.0,1  .HIOUNRANGE  ( ? T ) , 1 4, 7. 0 , 0 . 0 , XM  , XD,  2. 80,  ♦.50, 
1.1190,0. ,1.,-t,. 09)51 
CALL  L?N“(XI,YI,MI,1,0, 0) 

CALL  LIM,;(KS,YS»M?»l»3,l) 

CALL  5Y59DL (♦. 00,4.  00,  126,2,  0.,-1) 

CALL  PLOT (5. 30,3.4V, 3) 

CALL  P.3T(5.  55,3.  44,2) 

CALL  SYH90L  (5.  70,3.  <»P,  0 805 , 4HINTL,  0 . , 4) 

CALL  PLOTls. 30,3.31, 3) 

CALL  PLOT (5, 55, 3. 31,?) 

CALL  SYM33.  ( 5.  7.0 ,3.  27 , 0 005 , 3HOPT,  3 . , 3 ) 

CALL  SYM33. (5.47,3,  307, . 0700, 1,0., -1) 

CALL  PLOT { - 1. 9,3) 

CALL  PLOT (3, Q,-.0,?) 

CALL  PLOT( 3. 0 , 4. 8,2) 

CALL  PLOT(-l. ,4.0,21 
CALL  P.OT(-t. ,-.8,2) 

C 

C*******  ORAW  VERTICAL  TRA  JE) T ORY  PLOT 
C 

IF(ZA0. ME. -33000.)  3D  TO  95 
RH  s 120000. 

RO  = -20000. 

ZM  * 0. 

ZO  » 10000. 

GO  TO  10“ 

9?  IF (ZAO. ME, -20000.)  33  TD  100 
RM  * 60000. 

»0  ■ -10900. 

ZH  « 0. 

TO  * 5000. 

GO  TO  105 
100  CONTI NJ- 

1Q5  TKNI+l)  » T9 

ZItNI*?)  * 73 

TS(MF*1)  « ZM 

TS(MF*2>  « ZO 

RAN5I (MI*1)  « RN 
RANGI  (MI >2)  « RO 
RANGS(XFM)  > RM 
RAMOS (Mr»2)  « RO 
RO (1)  * RAMGI(l) 

T0(1)  » -ZAO 
R0(2)  » 0. 

Z0(2)  « 0. 

RO (3)  • RM 
R0(4)  * RO 
ZO ( 3)  » TM 
T0(4)  « TO 

CALL  P.3T(14.  ,0.,-I) 

CALL  S*AXIS(0.,0.»13MALTITUOC  (FT) , 13,  4.0,90.  , ZM,Z3,-,  65,1, 30, 
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i it 10,30. ,i. ,-i  ,.080>) 

CALL  S’AXISd.  ,0.,t0M1'MGu  (FT)  ,-10, 5.  0,0.  0,RM»20,  2.59  , - 
1 1139,1.  .1.  , -It.  09051 
CALL  Lr*Ji(*vr,I,7I,Hl,1,0,0) 

CALL  LT‘r(?»NCS,Z5,Mr,l,3,t) 

CALL  91GHLM(P9,Z0,2,l> 

CALL  SXMTO.r,.  » 0. » > 1 25 » 2 , 0. »-l) 

CALL  P.O?n,  >0,3.3=.,  J) 

CALL  P.OTd,  95.3.34, 2) 

CALL  SM1D.  (0.  30.3.  31,  0805,i»HINTL,  0.  , »»> 

CALL  P.  0T ( 1 . >0  » 3. 31 , ?) 

CALL  PL9TO. 85,3.91, 2) 

CALL  SYM33.  n.  90,3,  77,  0 305, 3HOPT .0.  .3) 

CALL  Sr^lO. (0.  67,3.  817,.  070  0,1,0.  «-l) 

CALL  P.OT(-l  ,-.8,3) 

CALL  P„ OT ( 7 , 0 , -•  8,  ? ) 

CALL  P.OT(7, 0,A. 3,»> 

CALL  P.Or(-t.  , •»,  8, 2) 

CALL  P.OT(-l. ,-.9,2) 

C 

c*+'****  OPAH  VM33T  VS.  riH:  PLOT 
C 

VM  = -100. 

VO  = 53. 

IF(ZA0. Mi. -35000.)  53  TO  25 
TM  * 0.' 

TO  = 15. 

GO  TO  35 

25  IF( ZAO. Ml. -20000.)  33  TO  30 
TM  = 0. 

TO  = 10. 

GO  TO  35 
31  COMTINJE 

35  °TI (MT*1)  * TM 

PTI(MI*2)  = TO 

PTS(MFH)  « TM 

°TS (MF*2)  = TO 

VOTI(NIH)  = VM 
VOTI ( MI *2)  = VO 
vots(n-4-d  = vm 
VOTS(N-f2)  = VO 
VOTI  (D  = VOTI  (3) 

VOTI (2)  = VOTI (3) 

VOTS ( 1)  = V DTS ( 3) 

VOTSC2)  = V0TS(3) 

• CALL  PLV CIV  , 0.  ,-J) 

CALL  S°AXIS(9. , 0 • , 15 MVMOOT  (FT/SEC 2) , 15, V. 0,30. ,VM,VO,-. 
1 1190,91. ,t.»-l.. 03051 

CALL  S?AXIS(0.  ,0«,t0MTTME  ( SiC) ,-10 , 5. 0 ,0. 0 , TM, TO, 2.5) , - 
1.1110,0. ,1. , -l, .0935) 

' CALL  LTM-C»Tt, VOTI, MI, 1,0,0) 

CALL  LTM* (PTS,Vr)TS,M" ,1,3,1) 

CALL  CLOT (0,  ,2.  ,3) 

CALL  PLOT (a. ,2. ,2) 

CALL  Pi  OT(0.50,3.3V,3) 

CALL  PL OT (1 , 35 , 3. 3V , 2) 

CALL  STM-IO;.  (0.90,3.  90,  0 805 , 4HINTL,  0.  , 4) 

CALL  PL OT(], 50,3.81, SI 
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CALL  PLOT(3.35,3.5l,?5 

CALL  tY'OO. <1. 10,3. 77,  0 805 , 3H0PT , 0 . , 3) 

CALL  SVMOO.  (0.  87,3.  83?.  . 07 0 0 , 1 , 0.  , - II 
CALL  P.0T<-1. ,-.8,3) 

CALL  PLOT (7.0, -.3, 2) 

CALL  P.0T('.  0,4. 8, 2) 

CALL  P.OTC-1. ,4.6,2) 

CALL  P_OT(-l  ,-.8,2) 

c 

C*,ut”  DTAH  ah;  VS.  TIME  9.0T 

c 

AM  = -20. 

A0'=  10. 

AHGI  ( MI  H)  s AH 
AHGI (HI »2)  = AO 
AHGS  (Mr  1)  = AM 
AHGS(NC«-?)  = AO 
CALL  P.OT ( l 4 . ,0.  ,-3) 

CALL  S»AXI3(0.  ,0.,3H*,Hr.  (G  » S)  , 9 ,4.  0 , 30  . , AM,  AD  , SO  , i.  40  , 

1 1190,30  030?) 

CALL  S°AXI3(0.,0.,10HrrME  ( St C) ,- 10 , S, 0 , 0. 0 , T M , TO, 2. 50 , -. 52 30 

1.1100.0.  .1. ,-l,.0305) 

CALL  LINK3TI,  AHGI,  NT,  1,0,0) 

CALL  LIM-(=TS, AHGS.M- ,1 ,3,1) 

CALL  P.OTtO, ,2. ,3) 

CALL  P.OT(S, ,2. ,2) 

CALL  P_OT(5. 30,3.94, 3) 

CALL  P'.OT(5.S5,3.34,2) 

CALL  SMTD.  (3.  70,3.  91,  0 805 .4HINTL, 0 . , 4) 

CALL  P_OT{5, 30,3.31, 3) 

CALL  P. OT(j.S5»3.81,2) 

CALL  STMOD. ( 5. 70,3. 77,  0305,3HOPT,0.  ,3) 

CALL  SfMOO. (5.  47,3.  8 3 7, , 070 0,1,0., -1) 

CALL  P.OT(-l  ,-.8,3) 

CALL  P-DT(7.0,-.8,2) 

CALL  PLOT  ( 7 . 0.,  4 . 3, 2) 

CALL  PL OT ( - 1 . ,4.8,2) 

CALL  P.OT(-t. ,-.8,2) 

**•*  ORA<  AV3  VS.  TIMS  P.OT 

AVGI(NTH)  .=  AM 
AVGI  (MTf  2)  = AO 
AVGS(N-H)  = AM 
AVGS  ( (lF*  2)  = AO 
CALL  P.OTU4.  ,0.  ,-3> 

CALL  PPAXIS  (0.  , 0..3MAV*.  (G  • S)  ,9,4.  0 , 30 . , AM,  AD, -.  SO  , 1,  40  , 

1 1190,30. ,l.,-l,. 0305) 

CALL  S3AXI3(0.  ,0.,104TrME  ( SEC) ,- 10 , 3. 0 , 0 , 0 , T M , TO, 2. 50 , -. 5290 

1.1130.0.  .1. ,-l,. 0305» 

CALL  LIM“(3TI» AVGI, NT , t , 0 , 0 ) 

CALL  LtM;(»T3, AVGS, N-, 1,3,1) 

call  plot ( o . ,2  ,3) 

CALL  P-OTJS. ,2. ,2) 

CALL  PLOT(5.30 ,3.94, 3) 

CALL  P. OT ( 5 , f>5 ,3.94,21 

CALL  SVM9D- (>.  70,3.  90,  0805.4HINTL, 0.  ,4) 

CALL  P.OT(5. 30,3,31, 3) 


o o o 
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CALL  PL0T(5.65,3.8l,?) 

CALL  SYMOD.  (3.  70,3.  77,  0 505  , 3 MOPT  , 3 . , 3 ) 

CALL  SYMBD.  (3.47,3.807,  . 07  0 0,  1,0.  , -1) 

CALL  PLOT < - 1 ,-.8,3) 

CALL  °LOT(7#0»-.8,?) 

CALL  P.OT(7.0,4.8»?) 

CALL  P. OT ( - 1 . ,4.8,?)  ' 

CALL  PLOT(-l. ,-.8,?) 

**•  ORAM  AL VS.  TIME  PLOT 

ALPitnr^i)  = am 
ALPI  ( MI  »-2)  = AO 

ALPS  (Nr*-l)  = AM  • 

ALPS(Nr*2)  = AO 
CALL  P.DT(U  ,0.,-3) 

CALL  SPlXIS(f).  ,0,,lHA.PHA  (OEG) , 11 , 4,  0 , 90.  , AM  , AO  , - . 60 , 1. 30, 

1 1190,00. ,1. ,-l,. 0303) 

CALL  SPAXIStO.  ,o.,10UfM£  (SEC)  ,-10,5.  0,0.0,  TM,  TO,  2.  SO, -.5290, 
1 1100,0. ,1. ,-l,. 08)5) 

CALL  LI*IC(»TI,  ALPI,  MI,  1,0,0) 

CALL  LIMKPTS,  ALPS,  M?, 1,3,1) 

CALL  P.OT(0. ,2. ,3) 

CALL  P.OT( 3 « ,2, ,2) 

CALL  P.0T(5, 30 ,3. 94, 3) 

CALL  P'.OT(5. 65,3.94,2) 

CALL  SYMBD. (5.70,3. 90,  0 805 , 4HINTL , 0 . , 4) 

CALL  P.OT(5, 30,3.31, 3) 

CALL  P. OT(5.o5,3.8l,2) 

CALL  SYMBD. (5. 70,3. 77,  0 805,3HOPT,0. ,3) 

CALL  SYMBD. (5.  47,3.  3 0 7,.  0700,1,0. ,-i) 

CALL  PLOT ( - 1 . ,-.8,3) 

CALL  PLOT (7,0, -.8, 2) 

CALL  PLOT(7, 0,4. 8,2) 

CALL  PLOT(-l. ,4. 8,2) 

CALL  pLOT ( - 1 . ,-.8,2) 

C 

C**»»*‘*  ORAM  BETA  VS.  TIM-  PLOT 
C 

4EPKNI  + 1)  = AM 
9£PI  ( NI  *2)  » AO 
9EPS(NrH)  * AM 
BEPS(Nr*2)  = AO 
CALL  P.OT ( L 4( , 0.  .*  3) 

. CALL  S°AXI 5 ( 0. ,0»,L1H3ETA  ( OEG)  ,11,4,  0,90., AM, AO, -.60, 1.30, 

1 1190, 90. , l. ,-l,. 0)05) 

CALL  SPAXIS(0. ,0.,10MT!ME  ( SEC) ,-10 , 5,  0 , 0. 0 , TM ,TO,  2.50  , 5290 , 
1.1190,0. ,1. ,-l ,. 0305) 

CALL  LIM-(»TI,REPI,*lli  1,  0,0) 

CALL  LrN€(PTS,OEPS,N-, 1,3,1) 

CALL  P'.OTO  » «?.  ,3) 

CALL  P.OT(b, ,2, ,2) 

CALL  PLOT(5. 39,3.94,3) 

CALL  P. OT(3.65,3.9*»  2) 

CALL  SYMBOL ( >• 70,3. 90,  0305,4HINTL, 0. , 4) 

CALL  P.  OT  ( 5 . 30 , 3. 31 , 3) 

CALL  P.OT (5. 65,3.91,?) 

CALL  SYMBOL (5.  7 0,3,  77,.  P 805 ,3H0PT , 0 . , 3) 


wr*r***  -v" 


CALL  SrMT)L(5.  47 , 3.  83 7 , , 070 G , 1 , 0. , -1> 

CALL  PLOT(-l.  ,-.8,3) 

CALL  P.OT(7,1, -.9,2) 

CALL  P.OT(7.0,4.8,2) 

CALL  PL  QT ( - 1 . ,<..8,2) 

.CALL  P.  OT  ( • 1.  ,-.9,2) 

' c 

OP  Ad  OH  VS.  RANGE  3.0T 
C 

IFtZAO.N", -35000.)  GD  TO  75 

OH  = -1000. 

00  = 1000. 

GO  TO  35 

75  IF(ZA0.  Hr. -20000.)  GD  TO  90 
OM  = -1000. 

00  = moo. 

GO  TO  35 
90  COMTIHJ? 

95  POHKNC4-U  = OH 

P0Hioir«-2)  = 5o 

®OHS(Hr*l>  = OM 

POHS(N-f2)  = 00 

CALL  P.  0 T (14.  , 0.  ,-3) 

CALL  S3 4X13(0.  ,ti.»7H0H  ( FT) ,7 ,4.0 , 9 0 . , DM, DO , 30 , 1, 50, 

1 1190,90. , l. , - 1, . 03  0 5) 

CALL  S34XIS(0. ,0.,10HR4NGE  (FT) ,-10,5. 0,0.0,RH,R0,  2.50, -.5290, 
1- 1190,0. ,1. ,-l ,. 09)5) 

CALL  LTN-<R4nGI,PDU, NT,  1,0,0) 

CALL  LIH"{RAMGS,POHS,HF, 1,3,1) 

CALL  PL  OT ( 0 . ,1. ,3) 

CALL  P. 0T ( 3 » ,1. ,2) 

CALL  P. 07(5*30,3.94,3) 

CALL  p„OT( 5 . 55 , 3 • 94, 2) 

CALL  SYM3D. (5. 70,3. 90,  0 805, 4HINTL , 0 . , 4) 

CALL  P.DT(5. 30,3.91, 3) 

CALL  P.0T{5.53,3.5t,2) 

CALL  SYMCO.  (5.  47,3.  917,  . 0700,1,0.  ,-D 
CALL  SYM3D- (5. 70,3. 77,  C805,3HOPT,0, ,3) 

CALL  P.0T(-1.  ,-.8,3) 

CALL  P.  07  ( 7 » 0 , -•  8*  2) 

CALL  PLOT(7.0,4.8,?) 

CALL  PLDT(-1. ,4.8,2) 

CALL  P.07 (-1. ,-.8,2) 

C 

£».*•••«  0PA(<  oz  ,/S.  RANGE  P.OT 
C 

IFCZAO. N\. -35000.)  GD  TO  110 
CM  = 1000. 
f)D  * -500. 
r.0  TO  120 

110  IF(740. ME. -20000.)  GD  TO  115 

CM  * 500. 

CO  * -250. 

GO  TO  120 
113  CONTINUE 
120  ®OZI ( NT  ♦ 1)  - OM 
®07I  ( NT  ♦ ?)  s 00 
pOZS( Nc»l)  = OM 
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onzs  (»)*«■?>  =oo 

CALL  PLOT (l 4.  ,0.  ,-3) 

CALL  S»4XI5(0.  ,0.,M3Z  C TT)  , 7 , 4.0 , 3 0 . , 0M,  3D, SO  * l . 50» 

1*.  1190  ,°9 . , 1 • *-l».(H05) 

CALL  ”4X13(0. ,0.,10-m‘JGE  (rT) 10 , i. 0 , 0. 0 , <«,  2. 5) , 5210 , 
1 ,1.  ,-i  0305) 

CALL  L!‘J2(*4‘I3I, P32I, ST  , 1,0,0) 

CALL  LTN2<?.MIGS,P9?S,N-,  1,3,1) 

CALL  PLOT(0. ,2.  ,3) 

CALL  PL  0 T ( a , , 2 . ,2) 

CALL  P.0T(5, 30 ,3.34,3) 

CALL  P'.CT  ( 5 . 35 , 3.94, 2) 

CALL  SYM13-_  (5.  70,3.  90,  P 605 ,4  HINTL » 0.  , 4) 

CALL  P.0T<;.3n,3.3l,!> 

CALL  P.OT(3.S5, 3.81,2) 

CALL  SYM-13.  (3.  70,3.  77,  0 805 , 3HOPT  , 0.  , 3 ) 

CALL  SYilD. (5. 47,3.  30  7,.  070  0,1,0.  t-1) 

CALL  P'.OM-l  ,-.8,3) 

CALL  PL07<7, o, -.8,2) 

CALL  PLOT ( 7 , 0 , 4. 8, 2 ) 

CALL  o,.0T(-l  ,4.8,2) 

CALL  P. O T ( - 1 . ,-.8,2) 


call  plot* i i oummy) 

$ETU*N 

2ND 
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Appendix  D 

Retran  Simulation  Parameter  Plots 


Figs.  58-64  are  parameter  plots  for  a 15°,  35,000  ft  altitude 
launch.  Figs.  65-71  are  plots  for  a 90°,  35,000  ft  altitude  launch. 
Figs.  72-78  are  plots  for  a 15°,  5,000  ft  altitude  launch.  Figs. 
79-85  are  plots  for  a 45°,  5,000  ft  altitude  launch.  Figs.  86-92 
are  plots  for  a 90°,  5,000  ft  altitude  launch. 


i 
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Trajectory  Planview  - 15°,  35,000  ft  Launch 


DOWNRRNGE  (FT) 


DOWNRRNGE  (FT) 


Fig.  72  Trajectory  Plan  View  - 15°,  5,000  ft  Launch 


Fig.  79  Trajectory  Plan  View  - 45 
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Fig.  86  Trajectory  Plan  View  - 90°,  5,000  ft  Launch 


Several  horizontal  guidance  algorithms  were  tested  In  the  Retran 
simulation.  Initially,  the  study  concentrated  on  launches  from  low 
initial  squint  angles  of  15°  to  30°  (squint  angle  is  the  angular 
displacement  of  the  line  of  sight  from  the  aircraft  velocity  vector) . 
Later,  emphasis  was  shifted  to  medium  to  high  initial  squint  angles 
of  30°  to  90°.  The  algorithms  discussed  in  section  II  were  found  to 
yield  the  largest  usable  launch  envelope  and  the  longest  maximum  launch 
ranges.  The  algorithms  presented  here  generally  resulted  in  poorer 
missile  performance. 


Straight  Ahead  Intercept 


For  launches  at  lew  initial  squint  angles,  the  missile  tends  to 


fall  behind  the  line  of  sight  in  the  early  portion  of  the  flight. 
Therefore,  the  first  algorithm  tested  coins anded  no  horizontal  accelera- 
tion until  line  of  sight  intercept,  then  commanded  a turn  to  place  the 
missile  on  the  line  of  sight.  A typical  trajectory  is  illustrated  in 
Fig.  93. 


This  algorithm  required  determination  of  the  point  at  which  to 
start  the  intercept  turn.  Determination  of  this  point  is  complicated 
by  the  deceleration  of  the  missile  during  the  turn;  also,  for  a randomly 
maneuvering  aircraft,  the  rotation  rate  of  the  line  of  sight  during  the 
time  of  the  turn  cannot  be  accurately  predicted  at  the  start  of  the 


turn.  To  make  a precise  intercept  turn  would  require  varying  the 
acceleration  caomands  throughout  the  turn,  and  the  necessary  algorithm 
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Fig.  93  Trajectory  Plan  View  - Straight  Ahead  Intercept 
Algorithm 


would  be  very  complex.  In  order  to  keep  the  algorithm  relatively 
simple,  a constant-g  intercept  turn  was  commanded.  A conservative 
estimation  of  the  start  turn  point  was  incorporated  which  would  allow 

f' 

completion  of  the  turn  with  the  missile  slightly  ahead  of  the  line  of 
sight.  The  rotation  of  the  line  of  sight  toward  the  missile  would 
then  result  in  a successful  Intercept. 

The  required  turning  angle  for  Intercept  is 

- <>»LOSst  + “z  *AtT  “ ’•'ST  - ’J'B  U°3) 

- $ * AtT  (104) 
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where 


( ) 


^losst 

wz 

AtT 

hr 


- heading  of  the  line  of  sight  at  start  turn 

- average  rotation  rate  of  the  line  of  sight 

- time  required  for  turn 

- missile  heading  at  start  turn 

- arbitrary  heading  bias  at  completion  of  turn 
(missile  heading  must  be  offset  slightly  from 
line  of  sight  in  order  for  missile  to  track 
the  rotating  line  of  sight) 

- average  missile  turn  rate 


The  missile  turning  rate  for  a constant  g turn  is 


£- 


Vv  3h 


(105) 


where  Vjj  is  missile  airspeed  in  the  turn  and  a^  is  the  commanded 
acceleration.  Assuming  constant  rotation  of  the  line  of  sight  and 
assuming  that  the  average  turn  rate  will  be  approximately  equal  to 
the  initial  turn  rate  (this  is  conservative,  since  the  turn  rate  will 
increase  as  airspeed  decreases  in  the  turn),  the  time  required  for 
the  turn  is 


Atr 


^cos 


ST 


- ♦ ST  - t 


B 


<P  - w. 


(106) 


In  order  to  start  the  turn  so  that  the  missile  will  complete 
the  turn  on  or  beyond  the  line  of  sight,  the  missile  must  have 
sufficient  overtake  speed  on  the  line  of  sight  to  close  the  distance 
to  the  line  of  sight  during  the  turn.  Here  it  was  estimated  that  the 
average  overtake  speed  would  equal  half  of  the  initial  missile  velocity 
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normal  to  the  line  of  sight.  The  approximate  arc  distance  normal  to 
the  line  of  sight  which  must  be  covered  during  the  turn  is 


L “ AyST  + * %T 


ST 


(107) 


where  AYgT  is  the  missile  distance  normal  to  the  line  of  sight  at 
start  turn  and  Rju  is  the  missile-to-target  range  at  start  turn 
(again,  this  is  conservative  because  will  decrease  during  the 
turn) . Therefore,  the  required  overtake  velocity  at  the  start  turn 
point  is 


Ay 


AtT  (AYST  + ^Z  ' RMTst) 


(108) 


• 

In  this  algorithm,  the  interception  turn  was  initiated  when  AY 
was  sufficient  to  insure  turn  completion  on  or  beyond  the  line  of 
sight.  A constant  4g  turn  was  commanded;  however,  if  during  the  turn 
the  current  A f was  insufficient  to  cover  the  distance  remaining  to 
turn  completion,  no  acceleration  was  commanded  and  the  missile 
continued  on  present  heading  until  overtake  was  again  sufficient  to 
continue  the  turn.  A command  to  track  algorithm  similar  to  the  one 
discussed  in  section  II  was  employed  for  line  of  sight  tracking  once 
tracking  heading  was  reached. 

This  algorithm  proved  to  have  only  limited  interception  capability. 
Successful  intercepts  were  made  from  an  initial  squint  angle  of  20°; 
however,  attempted  intercepts  from  larger  squint  angles  resulted  in 
the  missile  completing  the  turn  behind  the  line  of  sight,  with 
subsequent  hlgh-g  maneuvering  required  to  catch  the  line  of  sight. 

Even  when  Interception  was  successful,  tracking  time  on  the  line  of 
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Fig.  94  Trajectory  Plan  View  - Trun  to  Initial  Intercept 
Heading  Algorithm 


, ■ — ■ — ■ — ■ — — — ■ — 

At  ttg  < 10  seconds,  the  algorithm  reverted  to  command- to- track. 
This  algorithm  proved  Inferior  to  the  Initial  heading  algorithm  for 
low  squint  angle  launches,  as  the  reversion  to  command  to  track 
produced  a sharply  turning  trajectory  as  illustrated  in  Fig.  95. 
However,  the  algorithm  gave  superior  results  at  medium  initial  squint 
angles,  with  lower  terminal  airspeed  and  increased  tracking  time. 

A typical  trajectory  for  a medium  squint  angle  launch  is  shown  in 
Fig.  96.  It  was  observed  that  the  trajectories  produced  by  this 
algorithm  were  similar  to  those  which  could  be  obtained  by  utilizing 
the  command- to- track  algorithm  throughout  the  flight. 


Fig.  95  Trajectory  Plan  View  - Command  to  Fixed  Intercept 
Angle  Algorithm,  Low  Squint  Angle  Launch 
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airspeed  was  higher  than  obtained  with  the  previous  algorithms  and 
tracking  time  was  sufficient.  However,  at  initial  squint  angles 
greater  than  45*  the  missile  was  Initially  heading  away  from  the  line 
of  sight,  resulting  In  large  g turns  (6-8  g's)  to  bring  the  missile 
back  toward  the  line  of  sight.  Even  with  the  large  accelerations 
commanded,  sizable  distance  errors  developed  during  the  Initial 
turn;  consequently,  large  intercept  angles  coming  back  to  the  line 
of  sight  were  commanded.  This  resulted  in  high  g maneuvering  (10-15 
g's)  to  complete  the  interception  turn  onto  the  line  of  sight.  The 
trajectory  for  a 60°  angle  launch  is  shown  in  Fig.  99.  For  initial 
squint  angles  greater  than  60°,  the  required  maneuvering  sometimes 
exceeded  missile  capability  and  intercepts  were  not  accomplished. 
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Radar  - Retransmission  guidance  system  was  investigated  using  a digital  flight 
simulation.  A demonstration  flight  profile  was  assumed,  with  a minimum  of  20 
seconds  of  tracking  on  the  aircraft-to-target  line  of  sight  required.  A 
guidance  algorithm  was  developed  which  produced  satisfactory  trajectories.  A 
first  order  gradient  technique  was  employed  in  an  unsuccessful  attempt  to 
optimize  the  trajectories  for  maximum  ranop.  A useahle  launch  envplnne 
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for  this  missile  was  determined.  The  azimuthal  extent  of  the  envelope  was 
limited  by  radar  system  constraints  to  15°-90°  from  the  aircraft  velocity 
vector.  A maximum  slant  range  of  21  nautical  miles  was  obtained  from  a 
launch  altitude  of  35,000  ft.  Range  deteriorated  rapidly  with  decreasing 
launch  altitude,  with  5,000  ft  being  the  lowest  altitude  at  which  a useable 
launch  envelope  was  obtained.  Maximum  slant  range  at  5,000  ft  launch  altitude 
was  5.5  nautical  miles. 


